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Abstract 

Research  has  shown  that  perfonnance  of  liquid  hydrocarbon  fueled  pulse 
detonation  engines  is  limited  by  the  time  required  to  evaporate  liquid  fuel  droplets  within 
the  mixture.  Vaporization  of  liquid  fuels  prior  to  injection  has  been  shown  to  decrease 
ignition  times  and  also  increases  fuel  efficiency;  however,  the  size  and  efficiency  of  the 
vaporization  system  used  are  not  feasible  for  use  in  future  pulse  detonation  aircraft 
concepts.  The  purpose  of  this  research  is  to  harness  the  waste  heat  of  pulse  detonation 
engine  thrust  tubes  to  generate  a  steady-state,  self-sustained  flash  vaporization  and 
supercritical  heating  system  using  JP-8  as  the  working  fluid  and  fuel. 

Using  a  pulse  detonation  engine  thrust  tube  mounted  heat  exchanger,  the 
successful  flash  vaporization  of  JP-8  has  been  demonstrated.  Additional  testing 
demonstrated  the  successful  heating  of  JP-8  to  supercritical  conditions  with  fuel  injection 
temperatures  over  760  K.  All  JP-8  flash  vaporization  and  supercritical  heating  tests  were 
sustained  by  the  heated  fuel  and  run  to  steady-state  conditions.  Heat  addition  rates  to  the 
fuel  of  up  to  7.7  kW  were  achieved  during  superheated  testing.  A  method  for 
experimentally  determining  supercritical  fluid  density  is  presented  based  on  the  findings 
of  the  supercritical  heating  tests. 
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PULSE  DETONATION  ENGINE  THRUST  TUBE 
HEAT  EXCHANGER  FOR  FLASH  VAPORIZATION 
AND  SUPERCRITICAL  HEATING  OF  JP-8 

I.  Introduction 

Study  of  detonations  have  been  recorded  since  the  work  of  Hoffman  in  the  1940s 
(Hoffman,  1940).  However,  until  the  late  1980s  the  study  of  detonations  and  the  pulse 
detonation  engine  as  a  means  of  propulsion  had  seen  limited  interest.  Since  the  late 
1980’s  there  has  been  an  explosion  in  pulse  detonation  engine  research  rooted  in  the 
higher  thermal  efficiencies  of  the  constant  volume  process  which  detonations  closely 
emulate.  It  has  been  understood  for  sometime  that  the  constant  volume  process  has 
thennal  efficiencies  much  higher  than  that  of  constant  pressure  processes 
(Eidehnan,  1991)  used  in  most  of  today’s  current  aeronautical  propulsion  systems. 

In  addition  to  the  high  thermal  efficiency  of  pulse  detonation  engines,  benefits 
include  low  cost,  mechanical  simplicity,  few  moving  parts,  scalability,  and  a  wide  range 
of  operation.  Expected  applications  of  pulse  detonation  engines  include  cruise  missiles 
and  unmanned  aerial  vehicles.  Hybrid  concepts  using  pulsed  detonation  engines  as  an 
afterburner  in  turbojet  engines  or  as  an  additional  thrust  source  in  the  bypass  of  turbine 
engines  are  being  studied.  Other  research  efforts  include  combined  engine  concepts 
where  pulse  detonation  engines  are  used  up  to  hypersonic  velocities  at  which  time 
scramjets  engines  are  utilized  (Kailasanath,  2003).  Space  applications  include  pulse 
detonation  rocket  engines  which  are  currently  being  researched  and  tested  by  Air  Force 
Research  Laboratory  and  NASA  (Kailasanath,  2003) 
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Motivation 


While  the  prospective  applications  of  pulse  detonation  engines  are  extensive, 
there  are  numerous  technological  and  logistical  hurdles  that  must  be  overcome  before 
pulse  detonation  engines  may  transition  from  the  experimental  environment  to 
operational  use. 

The  majority  of  pulse  detonation  research  uses  hydrogen  or  gaseous  hydrocarbon 
fuels  (Glassman,  1996:224).  These  fuels  are  readily  available,  provide  excellent 
repeatability,  and  the  gaseous  state  of  the  fuel  contributes  to  excellent  detonability 
characteristics.  Conversely,  the  use  of  liquid  hydrocarbon  fuels  has  been  extremely 
limited  due  to  the  difficulty  in  obtaining  detonations.  While  there  has  been  limited 
success  in  using  liquid  hydrocarbons  and  various  aviation  fuels,  the  complexity  of  the 
systems  used  has  prevented  the  integration  of  liquid  fuels  as  the  standard  for 
experimental  pulsed  detonation  research.  As  the  maturity  of  pulse  detonation  engine 
technology  advances,  the  integration  of  liquid  hydrocarbon  fuels  is  paramount  to  the 
success  of  the  pulse  detonation  engine  as  a  viable  propulsion  system. 

Additionally,  the  United  States  Air  Force  and  Navy  have  invested  significant 
funding  and  research  into  pulse  detonation  research  in  hopes  of  high  military  payoff  for 
use  in  a  wide  array  of  aerospace  military  applications.  If  any  of  these  applications  are  to 
come  to  fruition,  the  use  of  military  grade  turbine  fuels  such  as  JP-8  and  JP-10  will  be 
essential.  The  processes  and  additives  meant  to  enhance  the  stability  of  these  fuels 
provide  additional  complications  by  further  decreasing  the  detonability  of  the  fuel. 
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Problem  Statement 


Recent  research  has  shown  that  the  difficulties  of  using  liquid  hydrocarbon  fuels 
and  even  military  grade  turbine  fuels  in  pulse  detonation  engines  are  surmountable  by  the 
use  of  complex  atomization  and  mixing  methods  or  by  the  use  of  flash  vaporization 
systems.  The  use  of  an  external  electrically  powered  flash  vaporization  system  was 
successfully  employed  to  flash  vaporize  JP-8  (Tucker,  2004)  While  these  systems 
demonstrated  the  use  of  the  liquid  hydrocarbon  fuels  of  interest,  the  feasibility  of 
incorporating  such  systems  on  aircraft  designs  are  not  practical.  The  focus  of  this 
research  is  to  further  the  development  of  a  practical  fuel  vaporization  and  supercritical 
heating  system  that  will  allow  the  use  of  military  grade  turbine  fuels  in  pulsed  detonation 
engines  without  the  use  of  complex  fuel  atomization  or  injection  methods  and  without  the 
use  of  an  external  power  source. 

Previous  Flash  Vaporization  Systems 

This  research  is  being  completed  as  a  direct  follow  on  to  works  completed  by  Dr. 
Colin  Tucker  and  sponsored  by  Air  Force  Research  Laboratories  Pulse  Detonation 
Research  Facility  (Tucker,  2004).  In  the  previous  flash  vaporization  system,  a  20  kW 
external  electric  heater  in  a  nitrogen  inert  furnace  was  used  to  statically  heat  pressurized 
JP-8  to  temperatures  well  above  the  auto-ignition  temperatures  of  the  fuel.  The  heated 
fuel  was  pressure  fed  through  fuel  nozzles  allowing  a  premixed  flash  vaporized  fuel/air 
mixture. 

Tucker’s  work  demonstrated  the  first  successful  detonation  of  JP-8  in  a  working 
pulse  detonation  engine.  The  research  was  also  successful  in  measuring  the  quantitative 
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benefits  of  flash  vaporizing  the  fuel,  in  addition  to  characterizing  the  required  parameters 
necessary  to  achieve  flash  vaporization  in  liquid  hydrocarbon  fuels,  specifically  JP-8.  To 
the  author’s  knowledge,  no  documented  flash  vaporization  system  has  been  used  in  pulse 
detonation  research  prior  to  Tucker’s  work. 

Pulse  Detonation  Heat  Transfer  Research 

There  have  been  a  number  of  papers  written  on  the  heat  losses  and  heat  loads  in 
pulse  detonations  engines  (Radulescu,  2004;  Hoke,  2003;  Eidehnan,  2000;  Paxson, 

2004).  The  focus  of  these  papers  has  been  to  determine  the  heat  loads  and  the  effect  on 
performance  due  to  heat  losses.  To  the  author’s  knowledge  none  of  these  papers  have 
used  liquid  hydrocarbons  as  the  detonating  fuel  or  emphasized  the  possibility  of  using  the 
heat  losses  as  a  potential  energy  source. 

Hoke  et  al.  (Hoke,  2003)  performed  a  calorimetrical  heat  load  analysis  to 
determine  how  heat  load  varied  as  a  function  of  multiple  operating  parameters.  These 
tests  provided  insight  into  heat  loads  generated  with  hydrogen  detonations  and  with 
shorter  tube  lengths  than  those  required  for  liquid  hydrocarbon  fuels.  Analytical  work  by 
Radulescu  et  al.  studied  the  effect  of  heat  losses  on  performance  and  how  the  heat  losses 
varied  as  a  function  of  the  length  of  tube  to  tube  diameter  (Radulescu,  2004).  Eidehnan 
et  al.  developed  simulations  to  model  wall  temperature  and  heat  transfer  as  a  function 
various  design  parameters  using  hydrogen  as  the  fuel  (Eidehnan,  2004).  Paxson  and 
Perkins  studied  the  cooling  required  for  a  hybrid  engine  with  pulsed  detonation  tubes 
acting  as  the  afterburner  in  gas  turbine  based  engine  (Paxson,  2004).  It  was  their 
conclusion  that  the  cooling  required  drastically  reduced  the  performance  benefit  of  the 


4 


hybrid  engine  concept.  The  cooling  methods  used  were  forced  or  free  air  convection  and 
varying  engine  parameters  rather  than  using  the  fuel  as  a  heat  sink. 

Research  Goals 

It  is  the  primary  goal  of  this  research  to  use  the  waste  heat  generated  by  the  pulse 
detonation  engine  cycle  to  develop  a  steady-state,  self-staining,  flash-vaporization  and 
supercritical  heating  system  for  military  grade  JP-8  turbine  fuel.  The  safety  risks  and 
experimental  uncertainties  in  such  a  venture  require  many  intermediate  goals  to  be  met  to 
ensure  the  feasibility  and  safety  of  meeting  the  primary  objective.  The  following  is  a  list 
of  the  intermediate  goals  met  in  order  to  achieve  the  primary  goal. 

1 .  Determine  if  sufficient  waste  heat  is  available  to  flash  vaporize  and  superheat 
military  grade  JP-8  turbine  fuel 

2.  Design  and  construct  a  PDE  tube  mounted  heat  exchanger  and  flash 
vaporization  system 

3.  Complete  water-cooled  tests  using  hydrogen  as  the  detonation  fuel  to 
determine  heat  transfer  and  prove  design  concept 

4.  Complete  aviation  gasoline  and  JP-8  detonation  tests  without  the  use  of  the 
heat  exchanger  to  determine  thrust  tube  heat  transfer  and  axial  temperature 
profdes 

5.  Complete  water-cooled  tests  using  aviation  gasoline  as  the  detonation  fuel  to 
determine  further  heat  transfer  estimates 

6.  Complete  safety  approval  process  for  using  a  JP-8  cooled  heat  exchanger 
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7.  Conduct  JP-8  cooled  heat  exchanger  tests  using  heated  JP-8  as  the  detonation 
fuel. 

8.  Determine  and  test  at  operating  condition  where  steady-state,  self-sustained, 
flash  vaporization  of  JP-8  occurs  using  a  JP-8  cooled  heat  exchanger  using 
heated  JP-8  as  the  detonation  fuel. 

9.  Determine  and  test  at  operating  condition  where  steady-state,  self-sustaining, 
supercritical  heating  of  JP-8  occurs  using  a  JP-8  cooled  heat  exchanger  using 
heated  JP-8  as  the  detonation  fuel. 

Chapter  Summary 

The  high  thermal  efficiencies  of  pulse  detonation  engine  propulsion  systems  have 
drawn  significant  attention  in  the  research  community.  One  of  the  largest  challenges  in 
the  advancement  of  pulse  detonations  engine  technologies  and  fielding  of  pulse 
detonation  engine  powered  vehicles  is  the  transition  from  the  use  of  gaseous  fuels  to 
liquid  hydrocarbons.  A  complex  flash  vaporization  system  has  been  demonstrated  as  a 
viable  method  for  using  liquid  hydrocarbon  fuels  It  is  the  purpose  of  this  work  to  use  a 
PDE  tube  mounted  heat  exchanger  harnessing  the  waste  heat  to  develop  a  flash 
vaporization  and  supercritical  JP-8  fuel  heating  system. 

Organization 

Chapter  I  served  as  a  brief  introduction  to  pulse  detonation  engine  technology.  In 
addition  the  motivation,  problem  statement,  and  goals  for  this  work  are  discussed. 
Chapter  II  provides  the  engineering  foundation  for  this  research  beginning  with  the 
explanation  of  the  detonations  and  the  concept  of  a  pulse  detonation  engine.  The 
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development  of  JP-8  property  data  and  the  heat  transfer  methods  used  are  also  presented. 
In  Chapter  III,  the  facility,  pulse  detonation  engine,  instrumentation,  and  test 
configurations  are  discussed.  Chapter  IV  provides  the  results  and  analysis  of 
experimental  data  and  an  experimental  method  for  determining  supercritical  density. 
Chapter  V  discusses  the  conclusions  from  the  previous  chapters  and  provides 
recommendations  for  further  research. 
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II.  Background 


Detonation  Overview 

When  a  tube  open  at  both  ends  is  filled  with  an  explosive  fuel-air  mixture  and 
ignited  at  either  end  the  flame  front  will  progress  into  the  flammable  mixture.  The 
products  of  the  reaction  behind  the  flame  front  have  a  higher  temperature  and  specific 
volume  in  comparison  to  the  unburned  mixture.  However,  because  the  tube  is  open  at 
both  ends  the  local  pressure  rise  created  by  the  increase  in  specific  volume  is  allowed 
equalize  with  the  ambient  environment  and  the  flame  travels  at  the  steady  rate  at  which  it 
can  burn  the  reactants.  The  flame  front  in  this  reaction  is  a  deflagration  wave. 

When  a  tube  closed  a  one  end  is  filled  with  an  explosive  fuel-air  mixture  and 
ignited  at  the  closed  end  the  flame  front  will  progress  into  the  flammable  mixture.  The 
products  of  the  reaction  trailing  the  flame  front  have  a  higher  temperature  and  specific 
volume  than  the  unburned  mixture.  The  increase  in  specific  volume  creates  a 
compression  wave  that  travels  at  the  speed  of  sound  until  it  reaches  the  flame  front 
causing  the  flame  to  accelerate.  As  the  flame  continues  to  bum  the  compression  waves 
are  continually  emanated  by  the  increase  in  specific  volume  behind  the  wave.  The 
compression  waves  also  create  an  increase  in  temperature  and  consequently  an  increase 
in  the  speed  of  sound  which  causes  each  compression  wave  to  catch  up  to  the  previous 
waves  and  coalesce  at  the  flame  front  further  increasing  the  speed  of  flame  front.  The 
increase  in  the  flame  speed  also  generates  turbulence  in  the  flame  which  causes  further 
acceleration  in  the  flame  front. 
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If  the  tube  is  long  enough  the  acceleration  of  the  flame  front  will  continue  until  a 
shock  is  formed  from  the  coalescing  compression  waves.  The  shock  wave  generated  is 
strong  enough  to  ignite  the  mixture  ahead  of  the  flame  front.  The  reaction  behind  the 
shock  front  continuously  sends  forward  compression  waves  preventing  the  shock  wave 
from  decaying.  This  self-sustaining  shock  wave  coupled  with  its  following  reaction  front 
is  termed  a  detonation  wave. 


Typical  ratios  of  detonations  and  deflagrations  are  presented  in  Table  1  for 
comparison,  with  state  one  and  state  two  represented  in  Figure  1.  Note  that  the 
detonation  pressure  ratio  is  an  order  of  magnitude  greater  than  deflagrations  waves  and 
also  that  the  detonation  waves  travel  at  Mach  numbers  several  orders  of  magnitudes 
higher  than  deflagration  waves. 


Table  1.  Typical  detonation  and  deflagration  property  ratios  (Glassman,  1996:223) 


Usual  magnitude  of  ratio 

Ratio 

Detonation 

Deflagration 

ul/al 

5-10 

0.0001  -0.03 

u2/u1 

0.4 -0.7 

4-16 

P2/P1 

13-55 

0.98  -0.976 

T2/T1 

8-21 

4-  16 

pl/p2 

1.4 -2.6 

0.06-0.25 

Detonation  Background 

Consider  a  one  dimensional  combustion  wave  that  travels  through  a  channel  from 


left  to  right.  Relative  to  the  wave,  the  reactants  moving  into  the  wave  at  a  velocity,  ui, 


and  the  products  are  leaving  the  wave  at  a  velocity,  m,  as  shown  in  Figure  1 . 
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Figure  1.  One  dimensional  combustion  wave  traveling  through  channel  with  velocities  relative  to  the 

wave  front 


The  one  dimensional  conservation  of  mass,  momentum,  and  energy  equations  are 


Plul  =  P2'u2 


Pi +p  flit  -P2+P2'u2 

1  2  1  2 
VT1+  9U1  +  *  =  cp'T2  +  T'u2 


If  the  ideal  gas  law  and  constant  gas  constant  can  be  assumed  for  both  the  reactants  and 


the  products  then 


pt =  p  i'R'Ti 


p2  -  P2'R'T2 


Additionally  the  speed  of  sound  and  Mach  number  relations  are 
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M 


m2  = 


U2 

a2 


(6) 

(7) 


al  =  ^-R'Tl  = 


Y’pl 

■J  pl 


a2  “  J?'  r'T2  - 

yj  P2 

By  substitution  and  re-arranging  the  above  equations  (Glassman,  1996:228)  it  can  be 
shown  that 


(8) 

(9) 


Y  -  1 


P2  _  )_1 

^P2  Pl  J  2 


(P2-Pl) 


+  ±1 

Pl  P2 )_ 


(10) 


Equation  (10)  is  known  as  the  Hugoniot  equation.  Assuming  values  for  density 
and  pressure  at  state  one  and  assigning  a  value  for  q,  the  possible  solutions  for  the 
pressure  and  density  of  state  two  can  be  detennined  by  setting  a  range  variable  for  one  of 
the  unknowns  and  solving  for  the  other.  In  this  manner  the  Hugoniot  curve  is  generated 
as  shown  in  Figure  2. 
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Sample  Hugoniot  Curve 


Figure  2.  Pressure  versus  inverse  density  for  initial  state  and  Hugoniot  curve  for  state  2 
By  drawing  parallel  lines  to  each  axis  through  state  one  and  also  drawing  lines 
from  state  one  to  the  tangents  of  the  Hugoniot  curve  for  state  two,  the  curve  can  be 
segmented  into  five  portions.  Only  certain  portions  of  the  Hugoniot  curve  represent 
solutions  that  are  physically  possible.  These  portions  can  be  further  differentiated  into 
steady-state  and  transient  solutions  represented  in  Figure  3.  Only  one  steady-state 
solution  occurs  at  pressures  greater  than  the  initial  state.  This  point  is  known  as  the 
Chapman-Jouguet  (CJ)  point  and  corresponds  to  a  self-sustained  detonation.  The  dashed 
region  above  the  CJ  point  refers  to  a  strong  detonation  which  is  an  unstable  condition  that 
will  equalize  to  the  CJ  point.  The  dashed  region  below  the  CJ  point  refers  to  a  weak 
detonation.  For  this  solution  to  exist  extremely  fast  chemical  kinetics  are  required.  The 
kinetics  of  liquid  hydrocarbon  fuels,  such  as  JP-8,  does  not  pennit  this  solution.  The 
solutions  to  the  right  of  the  state  one  refer  to  deflagrations.  (Glassman,  1996:232-235) 
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Sample  Hugoniot  Curve 
with  Physically  Possible  Solutions 

- -  Strong  Detonation 

(Transient  Condition) 

- Weak  Detonation 

(Not  possible  in  JP-8) 

-  Weak  Deflagration 

(Steady-State  Condition) 

Deflagration 

/ 

i/P 

Figure  3.  Pressure  versus  inverse  density  for  initial  state  and  Hugoniot  curve  with  physically 

possible  solutions  for  state  two 

For  propulsion  purposes  the  CJ  point  is  the  desired  steady-state  solution.  Of  the 
possible  solutions  the  CJ  point  corresponds  to  the  state  with  the  minimum  entropy 
(Glassman,  1996:236)  and  reaches  supersonic  detonation  wave  velocities  with  large 
pressure  ratios.  The  wave  speed  of  the  CJ  point  has  become  one  of  the  governing  metrics 
for  detennining  when  detonations  are  present.  The  CJ  wave  speed  for  most  hydrocarbon- 
air  mixtures  with  equivalence  ratios  from  1.0  -  2.0  are  between  1700  m/s  and  2000  m/s 
(Glassman,  1996;247) 

Pulse  Detonation  Engine  Cycle 

The  generalized  pulse  detonation  engine  (PDE)  cycle  is  composed  of  fill, 
detonation,  and  purge  processes  represented  in  Figure  4,  Figure  5,  and  Figure  6 
respectively.  In  the  fill  process  a  valve  is  opened  filling  the  PDE  tube  with  an  explosive 
fuel-air  mixture.  The  ratio  of  the  fuel-air  mixture  volume  to  the  tube  volume  is  referred 
to  as  the  fill  fraction.  A  fill  fraction  of  one  represents  the  entire  tube  filled  with  a  fuel-air 
mixture  while  a  fill  fraction  of  0.5  indicates  half  of  the  tube  is  filled.  For  the  purpose  of 


Chapman-Jouget 
Steady-State 
/  Detonation  Wave 
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this  work  the  fill  fraction  was  kept  at  one.  When  the  fill  process  is  complete  the  fill  valve 
is  closed. 


Figure  4.  Generalized  PDE  fill  process  with  valve  opening  to  fill  the  PDE  tube  with  fuel-air  mixture 


The  detonation  process  is  typically  thermally  initiated.  For  this  work  the  ignition 


energy  was  provided  by  an  automotive  spark  ignition  system.  The  ignition  begins  a 


deflagration  at  the  closed  end  of  the  tube  which  transitions  to  a  detonation  by  the  process 


mentioned  earlier  in  this  chapter.  The  transition  from  a  deflagration  to  a  detonation  is 


known  as  the  detonation  to  deflagration  transition  (DDT).  The  time  and  distance  required 
to  achieve  DDT  are  to  be  minimized  for  optimum  PDE  performance.  Thrust  for  the  PDE 
is  generated  during  blowdown  when  the  products  are  exhausted  by  the  finite  change  in 


momentum  and  also  due  to  the  pressure  ratio  at  the  exit  area. 
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Figure  5.  Spark  initiated  PDE  detonation  process  with  transition  from  deflagration  wave  to 

detonation  wave 


In  the  purge  process  the  purge  valve  opens  and  fills  the  tube  with  a  non-explosive 
air  mixture  to  purge  the  exhaust  products  from  the  tube.  The  purge  air  also  cools  the 
PDE  tube  and  remaining  products.  The  purpose  of  the  purge  process  is  to  prevent  a 
standing  flame  from  occurring  that  would  lead  to  a  backfire.  The  purge  fraction  is  the 
ratio  of  purge  air  volume  to  the  PDE  tube  volume.  The  purge  fraction  for  this  work  was 
0.5.  Once  the  tube  has  been  filled  with  the  designated  purge  fraction  the  purge  valve  is 
closed  and  the  PDE  cycle  is  repeated. 


Figure  6.  Generalized  PDE  purge  process  with  valve  opening  to  purge  the  PDE  tube  with  air 
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Flash  Vaporization 

Flash  vaporization  is  a  process  that  occurs  between  two  states.  The  initial  state 
required  for  flash  vaporization  is  characterized  by  a  fluid  at  an  elevated  temperature  and 
pressure.  The  second  state  is  characterized  by  a  significantly  lower  pressure  and  constant 
enthalpy  if  the  process  is  assumed  adiabatic.  If  the  pressure  drop  occurs  rapidly  the  fluid 
passes  through  the  vapor  dome  and  the  final  state  of  the  substance  will  be  a  vapor.  The 
state  flash  vaporization  temperature  is  defined  as  the  temperature  that  intersects  the  state 
pressure  with  the  saturated  vapor  line.  The  process  from  state  one  to  state  two  is 
represented  in  the  pressure-enthalpy  diagram  for  the  AFRL  SUPERTRAPP  JP-8 
surrogate  in  Figure  7.  The  development  of  the  data  for  this  figure  is  presented  later  in 
this  chapter. 


Pressure  vs.  Enthalpy 
AFRL  SUPERTRAPP  JP-8  Surrogate 


Enthalpy  (kJ/kg) 


Figure  7.  Pressure  versus  enthalpy  diagram  of  AFRL  SUPERTRAPP  JP-8  surrogate  with  vapor 

dome  and  flash  vaporization  process  path 

In  the  PDE  engine  the  fuel  is  injected  into  an  air  stream  through  the  use  of  fuel 
nozzles.  The  pressure  after  the  injection  point  may  be  assumed  to  be  at  ambient  pressure. 
Assuming  the  pressurized  heated  liquid  fuel  is  at  or  above  the  saturated  vapor 
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temperature  for  ambient  pressure;  the  pressure  drop  that  occurs  across  the  nozzles  will 
flash  vaporize  the  fuel  locally  in  the  fuel  stream.  If  the  air  stream  is  not  at  a  temperature 
at  or  above  the  dew  point  temperature  of  the  mixture,  some  of  the  fuel  in  the  mixture  will 
re-condense. 

The  concept  of  local  flash  vaporization  and  mixture  dew  point  temperature  is 
illustrated  in  a  generic  mixture  temperature  versus  fuel  temperature  diagram  shown  in 
Figure  8.  In  two  cases  the  air  is  heated  prior  to  fuel  injection  to  air  temperatures  Tairi  and 
Tair2  with  Tairi >Tair2.  For  the  fuel  temperature  range  the  mixture  temperature  for  the  case 
with  Tair2  does  not  exceed  dew  point  temperature  of  the  mixture  (Tdpmix)  and  does  not 
completely  vaporize  the  mixture,  but  it  does  locally  flash  vaporize  the  fuel  at  fuel 
temperatures  greater  than  or  equal  to  the  fuel  flash  vaporization  temperature  (Tfvfuei)- 
Over  the  same  range  of  fuel  temperature  for  the  case  with  Tairi  the  mixture  temperature 
surpasses  the  mixture  dew  point  temperature  and  the  fuel  mixture  becomes  completely 
vaporized. 
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Figure  8.  Mixture  temperature  versus  fuel  temperature  for  two  air  temperatures  at  constant  air 

pressure 

The  mixture  temperature  versus  JP-8  fuel  temperature  for  various  air  temperatures 
is  presented  in  Figure  9.  This  data  was  generated  analytically  and  experimentally  verified 


(Tucker,  2005).  For  JP-8  the  mixture  temperature  required  for  a  completely  vaporization 
is  400  K  (Tucker,  2005).  The  JP-8  flash  vaporization  temperature  shown  in  Figure  9  is 
530  K. 
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Figure  9.  Stoichiometric  JP-8  surrogate  air  mixture  liquid  vapor  equilibrium  in  the  intake  manifold 

for  4  air  temperatures  at  2  bar.  (Tucker,  2005) 

Supercritical  JP-8 

The  critical  point  is  the  point  where  the  slope  of  the  vapor  dome  is  zero  as  shown 
in  Figure  10.  The  critical  pressure  and  temperature  are  values  that  intersect  at  the  critical 
point.  For  pressures  and  temperatures  greater  than  or  equal  to  the  critical  pressure  and 
temperature  the  substance  can  no  longer  be  considered  a  liquid  or  vapor.  At  these  values 
the  substance  is  referred  to  as  supercritical.  The  pressure-temperature  diagram  used  for 
JP-8  is  presented  in  Figure  10.  The  development  of  the  figure  is  presented  in  the 
following  section.  The  critical  pressure  and  critical  temperature  of  JP-8  are  23.3  atm  and 
683  K  respectively  (Szetela,  1981). 
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Figure  10.  Pressure  versus  temperature  diagram  of  AFRL  SUPERTRAPP  JP-8  surrogate  with 

vapor  dome  and  critical  point 

JP-8  Thermodynamic  Properties  Data 

The  flash  vaporization  system  design  and  analysis  required  the  thennodynamic 
property  data  of  JP-8  from  room  temperature  to  supercritical  temperature  ranges  of  up  to 
800  K.  The  lower  temperature  property  data  can  be  obtained  from  the  CRC  Handbook  of 
Aviation  Fuel  Properties  (CRC,  2004).  Unfortunately,  the  property  data  is  only  provided 
for  temperatures  up  to  393  K  for  JP-8.  To  the  author’s  knowledge,  there  is  no  complete 
set  of  thermodynamic  data  for  JP-8  above  this  temperature  that  includes  density,  specific 
heat,  viscosity,  and  thermal  conductivity. 

The  lack  of  high  temperature  JP-8  property  data  is  attributed  to  two  reasons.  First 
and  foremost,  there  has  not  been  a  widespread  need  for  data  at  these  temperatures.  In 
most  engineering  applications  it  is  not  necessary,  or  desirable,  to  operate  with  fuel  at 
supercritical  temperatures  and  pressures.  Second  there  are  no  current  methods  to 
determine  the  fuel  properties  at  these  temperatures  and  pressures.  The  enthalpy  and 
specific  heat  of  high  temperature  fuels  may  be  determined  by  heat  transfer  tests  but 
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without  the  other  key  thermodynamic  properties  there  is  limited  benefit.  As  a  result  of 
this  work  an  experimental  method  for  detennining  high  temperature  fuel  density  is 
presented. 

Since  JP-8  is  not  a  pure  substance  it  also  difficult  to  model  the  fuel  based  on  its 
chemical  makeup.  To  model  characteristics  or  properties  of  JP-8  a  surrogate  was  used. 

A  surrogate  is  a  mixture  of  hydrocarbons  that  are  pure  substances  for  which  the  chemistry 
is  known.  The  surrogate  is  developed  to  mimic  certain  characteristics  or  properties  of  the 
fuels. 

PPDS2  and  SUPERTRAPP  are  two  known  software  packages  that  predict  the 
thennophysical  property  data  based  on  a  surrogate  mixture.  PPDS2  is  a  commercially 
available  thennophysical  properties  calculation  engine  and  associated  databanks  for 
obtaining  the  thermodynamic,  transport  and  phase  equilibria  properties  for  fluids  and 
fluid  mixtures  (PPDS,  2005).  SUPERTRAPP  is  a  product  of  the  National  Institute  of 
Standards  and  Technology  (NIST).  SUPERTRAPP  is  an  interactive  computer  database 
for  the  prediction  of  thennodynamic  and  transport  properties  of  fluid  mixtures  (NIST, 
2003). 

For  this  work  SUPERTRAPP  was  used  to  develop  the  required  thermodynamic 
property  data  for  a  surrogate  JP-8  mixture.  The  surrogate  used  in  the  development  of  the 
data  was  provided  by  Air  Force  Research  Laboratory  Propulsion  Directorate  (AFRL/PR). 
The  surrogate  composition,  hereto  referred  to  as  the  AFRL  SUPERTRAPP  JP-8 
surrogate,  is  presented  in  Table  2. 
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Table  2.  AFRL  SUPERTRAPP  JP-8  surrogate  composition  (Spadaccini,  1998) 


Component 

SUPERTRAPP  Synonym 

Mole  Fraction 

Mass  Fraction 

methylcyclohexane 

MCC6 

0.065 

0.042 

meta-xylene 

MXYL 

0.07 

0.049 

ethylcyclohexane 

ECC6 

0.067 

0.050 

n-decane 

CIO 

0.157 

0.147 

butyl  benzene 

C4BNZ 

0.056 

0.050 

isobutylbenzene 

IC4BNZ 

0.056 

0.050 

nampthalene 

NAPH 

0.058 

0.049 

n-dodecane 

C12 

0.175 

0.207 

1-methylnapthalene 

1MNAPH 

0.052 

0.049 

n-tetradecane 

C14 

0.113 

0.156 

n-hexadecane 

C16 

0.066 

0.104 

2,5-dimethylhexane 

25DMH 

0.065 

0.049 

The  determination  to  use  the  AFRL  surrogate  was  based  on  agreement  with  the 
CRC  Handbook  of  Aviation  Fuel  Properties  in  the  available  temperature  range  and  with 
other  surrogate  model  data  available  at  the  limited  pressures  and  temperatures  available. 
The  comparison  surrogate  data  was  also  provided  by  AFRL/PR.  The  first  comparison 
surrogate  data  was  compiled  in  SUPERTRAPP  based  on  a  surrogate  developed  by  Schulz 
(Henegan  and  Schulz,  1993).  The  second  surrogate  data  was  compiled  in  PPDS2  based 
on  a  proprietary  surrogate.  Both  of  these  data  sets  were  generated  for  a  pressure  of  6.895 
MPa. 


The  AFRL  SUPERTRAPP  data  generated  compared  well  with  low  temperature 


CRC  data  and  with  the  two  comparison  surrogate  data  sets  at  the  prescribed  pressure  of 


6.895  MPa  as  shown  in  Figure  11. 
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Figure  11.  Modelled  JP-8  surrogate  density  at  6.895  MPa  and  CRC  JP-8  comparison 
The  AFRL  SUPERTRAPP  surrogate  presented  in  Table  2  met  the  measures  of 
success  prescribed  and  is  used  for  all  design  calculations  and  analysis.  The  surrogate 
shall  be  referred  to  as  the  AFRL  surrogate  for  the  remainder  of  this  document.  Density, 
specific  heat,  viscosity,  and  thermal  conductivity  tables  were  generated  for  temperatures 
ranging  from  273-998  K  and  pressures  ranging  from  1-85  atm  using  the  AFRL  surrogate 
in  SUPERTRAPP.  Using  the  state  information  provided  by  SUPERTRAPP  the  vapor 
pressure  curve  was  also  generated.  All  AFRL  surrogate  data  generated  by 
SUPERTRAPP  is  provided  in  Appendix  A. 


Power  Required 

For  steady-state  conditions  with  no  work  and  negligible  contributions  from 
velocity  and  elevation,  the  specific  energy  required  to  heat  JP-8  over  a  temperature 
differential  can  be  determined  from  the  reduced  form  of  the  energy  equation  (Incropera 
and  DeWitt,  1996:399) 
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q  = 


'^outlet 

cp(T)dT 

"Millet 


(11) 


Graphically  this  equation  represents  the  area  under  the  curve  between  any  two 
temperatures  on  a  specific  heat  versus  temperature  curve.  With  the  specific  heat  data 
generated  for  the  AFRL  SUPERTRAPP  JP-8  surrogate  the  finite  difference  method  was 
used  to  determine  the  specific  energy  required  between  any  two  temperatures.  For  the 
purpose  of  this  work  the  cold  temperature  was  conservatively  set  to  290  K.  The  hot 
temperature  chosen  for  the  initial  heat  exchanger  design  was  530  K.  The  heat  of 
combustion  of  the  JP-8  used  in  the  tests  was  43.3  MJ/kg.  From  Eq.  1 1  the  specific 
energy  required  is  554  kJ/kg  or  1.3%  of  the  energy  available  in  the  fuel. 


To  determine  the  power  required,  the  mass  flow  must  be  calculated.  The  mass 
flow  calculation  is  determined  from  the  engine  operating  parameters  and  the  chemistry  of 
the  fuel-air  mixture.  The  engine  operating  parameters  provide  the  volumetric  flow  rate 
and  mixture  chemistry  provides  the  correlation  to  mass  flow  rate. 


The  volumetric  flow  of  the  fuel-air  mixture  is  a  dependent  on  engine  frequency, 
fill  fraction,  number  of  PDE  tubes,  and  tube  volume.  Engine  frequency  is  a  user  input 
parameter  and  for  design  purposes  is  constant  at  15  Hz.  The  tube  volume  for  one  1.829 
m  tube  with  a  5.25  cm  diameter  is  0.004  m  .  Mechanical  limitations  of  the  current  PDE 
engine  require  a  two-tube  configuration  for  liquid  fuel  operation.  Since  only  one  heat 
exchanger  is  to  be  installed  during  the  tests  for  this  work  the  mass  flow  requirements  will 
be  based  on  a  two-tube  configuration.  Fill  fraction  was  unity  for  this  work.  The  total 
volumetric  flow  of  the  engine  is  then  determined  by 
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- -  =  (Frequency  )-(Tube_Volume)-(Fill_Fraction)-(Number_of_Tubes  ) 

dt  (12) 

The  equivalence  ratio  is  the  ratio  of  the  design  or  test  fuel-air  mixture  divided  by 
the  ratio  of  fuel-air  mixture  for  stoichiometric  conditions  (Glassman,  1996:21).  An 
equivalence  ratio  of  unity  is  tenned  stoichiometric.  A  stoichiometric  mixture  is  the  ideal 
ratio  of  fuel  and  air  so  that  all  of  the  air  and  all  of  the  fuel  are  consumed  in  the  reaction. 


f  mfuel  ^ 
-  act 

f  “fuel  ^ 
-  sto 

^  “air  J 


(13) 


Equivalence  ratios  of  greater  than  one  indicate  there  is  not  sufficient  air  to  allow 
for  complete  combustion  of  the  fuel  and  is  referred  to  as  rich.  Equivalence  ratios  of  less 
than  one  indicate  there  is  excess  air  for  the  fuel  present  and  is  referred  to  as  lean. 


Since  JP-8  is  not  a  pure  substance  the  chemical  formula  for  its  composition  is 
unknown  and  in  most  cases  proprietary  to  the  manufacturer.  The  chemical  formula  is 
approximated  as  C10.9H20.9.  To  determine  the  stoichiometric  ratio  of  JP-8  and  air  the 
following  chemical  equation  must  be  balanced  for  Y,  A,  B,  and  C. 

1-(C10.9'H20.9)  +  Y' 

Solving  for  A,  B,  C,  and  Y  yeilds 


Jn2]  =  a-(c-o2)  +  b.(h2.o)  +  c-(n2) 


A  =  10.900  B  =  10.450  Y=  16.125  C  =  60.66 
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The  mass  of  the  fuel  and  air  are  obtained  by  multiplying  the  molar  values  above  by  the 
molecular  weight.  The  stoichiometric  ratios  of  air-fuel  and  fuel-air  mixture  for  JP-8  are 
obtained 


“air  ^ 

-  sto  =  14.575 

“fuel ) 

(15) 

“fuel  ^ 

sto  =  0.069 

“air  ) 

(16) 

With  the  desired  equivalence  ratio,  stoichiometric  fuel-air  mixture,  volumetric 
flow  rate  and  assuming  steady-state  conditions  with  ideal  gas  law  for  air,  the  volumetric 
and  mass  flow  for  the  air  and  fuel  can  be  detennined  from 


Xdot_mix  \iot_air 

mdot_fuel  f 

- =  ' 

mdot  air  1 


+  ^dot_fuel 


“fuel  ^ 
-  sto 

“air  ) 


“dot  air 


V  .p 
dot_air  amb 

*\iniv 


MW  air 


“dotfuel  Vdot_fuel '  P  sur 


(17) 

(18) 


(19) 

(20) 


The  fuel  mass  flow  multiplied  by  the  specific  energy  yields  the  power  required  to 
heat  the  fuel  through  the  desire  temperature  range. 


Power_Required  =  m(jot  ^ej  -Specific_Energy 


(21) 


The  specified  and  calculated  heat  exchanger  design  parameters  are  presented  in  Table  3. 
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Table  3.  Heat  Exchanger  Design  Parameters 


Frequency 

15  Hz 

Single  Tube  Volume 

0.004  mA3 

Number  of  Tubes 

2 

Fill  Fraction 

1 

Equivalence  Ratio 

1.05 

Inlet  Temperature 

290  K 

Outlet  Temperature 

530  K 

Specific  Energy  Required 

554  kJ/kg 

Fuel  Mass  Flow 

0.443  kg/min 

Heat  Exchanger  Power  Required 

4.082  kW 

PDE  Heat  Transfer  Coefficient  and  Inner  Tube  Temperature 

Approximate  values  for  the  average  heat  transfer  coefficient  and  the  average 
internal  tube  temperature  were  required  to  begin  the  heat  exchanger  design.  Due  to  the 
high  temperatures,  vibrations,  and  impulses  generated  in  the  PDE  tube  these  values  could 
not  be  experimentally  determined.  These  values  were  extrapolated  from  experimental 
data  from  previous  hydrogen-air  tests  for  both  plain  tube  steady-state  external  tube 
temperatures  and  calorimetric  heat  transfer  testing  with  a  water-cooled  heat  exchanger 
(Hoke,  2003).  Using  the  experimental  data  the  heat  transfer  rates  were  calculated.  Using 
conservation  of  energy  the  heat  transfer  rates  are  used  to  calculate  the  inner  wall 
temperatures.  With  the  heat  transfer  and  inner  wall  temperature  for  both  tests  the  heat 
transfer  coefficient  and  inner  tube  temperature  may  be  solved  for  as  the  two  unknowns. 

For  the  referenced  plain  tube  tests  the  steady-state  external  wall  temperature  at  the 
hottest  section  of  the  tube  was  1005  K.  The  heat  transfer  per  unit  length  was  calculated 
based  on  free  convection  and  radiation  losses.  The  complete  calculations  can  be  found  in 
Appendix  B.  All  material  and  air  properties  were  calculated  based  on  table  values  and 
linearly  interpolated  to  the  film  temperature  (Incropera  and  DeWitt,  1996:326)  defined  as 
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Where 


(22) 


(tq  +  Tamb) 
Tfilm= - I - 


T 

0  =  External  surface  temperature  (K) 
T 

amb=  Ambient  temperature 


Radiation  losses  were  calculated  based  on  heat  transfer  to  a  black  body  (Incropera  and 
Dewitt,  1996:10) 


Where 


chad 


£-a-d0-ji- 


(23) 


8  =  Emissivity  of  the  tube  material  from  property  table 
(7  —  5.67 x  10”  8 
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s  K  (Stefan-Boltzmann  constant) 
cl'a°  =  Outer  tube  diameter  (m) 


Convective  losses  due  to  free  convection  (Incropera  and  DeWitt,  1996:8)  are  determined 
by 


Where 


*Ifc  hfc'^'^o'jTo  Tamb) 


(24) 


qfc  =  Heat  transfer  due  to  free  convection  (W/m) 

hfc  =  Free  convection  heat  transfer  coefficient  [W/(nr-K)] 
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The  free  convection  heat  transfer  coefficient  (Incropera  and  DeWitt,  1996:307)  is 
calculated  by 


Where 


(25) 


kamb  =  Thennal  conductivity  of  the  air  [W/(m-K)] 
Nufc  =  Nussclt  number  for  free  convection 


The  Nussult  number  for  free  convection  is  detennined  by  the  Churchill  and  Chu 
correlation  (Incropera  and  DeWitt,  1996:465)  for  a  long  horizontal  cylinder 

-,2 

1 

0.387-Ra6 

_8_ 

'  0.559s)  16 

v  Pr&  J  (26) 

The  Raleigh  number  (Incropera  and  DeWitt,  1996:456)  is  detennined  from 

g'P-(To-Tamb)'diao3 

Ra  = - 

v-a  (27) 

Where 

film  Expansion  coefficient  for  ideal  gas  (1/K) 
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Acceleration  due  to  gravity 


g  = 


9 

v  =  Kinematic  viscosity  (m  /2) 
a  =  Thennal  diffusivity  (nr/2) 

The  summation  of  the  radiation  and  free  convection  losses  equals  the  total  heat  loss 

^plairi-tube  -  ‘trad  +  <tfc 

For  steady-state  conditions  the  heat  transfer  to  the  ambient  environment  must  the  same  as 
the  heat  transfer  through  the  wall.  The  tube  is  considered  a  radial  system  and  the 
conduction  through  the  tube  wall  (Incropera  and  DeWitt,  1996:91)  may  be  represented  as 

“'7t''Hube'(^o  ~~  ^i_plain_tube ) 

Vainjube  =  /  dia.  ^ 

In  — 

VAiao)  (28) 

Where 

ktube  =  Thermal  conductivity  of  the  tube  wall  [W/(m-K)] 

T- 

1  =  Average  temperature  on  the  inside  of  the  PDE  tube  (  K) 
clut|  =  Inner  PDE  tube  diameter  (m) 


The  inner  surface  temperature  may  be  detennined  knowing  all  other  terms. 

For  the  calorimetric  water  cooled  heat  exchanger  tests  the  temperature  increases 
were  on  the  order  of  13  K.  Over  this  temperature  range  the  specific  heat  for  water  may 
be  assumed  constant  and  the  total  heat  transfer  (Incropera  and  DeWitt,  1996:399)  can  be 
determined  from 
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Qcal  ^Motwater '  cp_water  *  ( ^outlet 


Where 


Qeiil  =  Heat  transfer  from  calorimetric  heat  exchanger  tests  (W) 

cp_water  =  Specific  heat  of  the  water  (J/[kg-K]) 

T- 

inlet  =  Heat  exchanger  inlet  temperature  (  K) 

T 

outlet  =  Heat  exchanger  outlet  temperature  (  K) 


By  using  a  modified  fonn  of  Equation  28,  the  inner  surface  temperature  for  the 
calorimetric  water-cooled  heat  exchanger  tests. 


Qcal 


•^'^tube '  L-(T0 


In 
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Where 


L  =  Length  of  the  heat  exchanger  (m) 


Knowing  the  heat  transfer  and  inner  tube  temperature  for  the  free 
convection/radiation  test  and  the  water  cooled  test  there  are  now  two  equations  and  two 
unknowns  allowing  for  the  average  heat  transfer  coefficient  (hflame)  and  inner  tube 
temperature  (Tflame)  to  be  solved  for 


II  l«  linn  II  M3  II  !l!  IMIS  II  IN  lilt!  Ill  ft  mi; 
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The  calculated  heat  transfer  coefficient  and  inside  tube  temperature  were  calculated  to  be 

Aflame-  ^5  ^  Aflame-  H^IK 

n’"K  (33 

Heat  Exchanger  Design 

The  heat  exchanger  designed  for  this  work  was  a  concentric  tube  counter-flow 
heat  exchanger.  This  type  of  heat  exchanger  was  chosen  based  on  ease  of  incorporating 
into  the  PDE  tube,  size,  simplicity,  ease  of  manufacture,  and  the  high  heat  transfer 
efficiency  possible  by  having  the  fuel  in  contact  with  the  heat  transfer  source.  The  heat 
exchanger  was  designed  based  on  calculations  for  both  water  and  JP-8  as  the  working 
fluid.  The  general  methodology  and  assumptions  used  to  design  the  heat  exchanger  are 
presented  in  the  following  paragraphs.  The  actual  calculations  for  both  working  fluids 
are  presented  in  Appendix  B. 

The  heat  exchanger  design  is  based  on  a  one-dimensional,  steady-state,  uniform 
flow,  radial  system  using  a  finite  difference  method.  A  representation  of  a  heat 
exchanger  segment  and  cutaway  view  is  presented  in  Figure  12. 


Figure  12.  Representative  concentric  tube  heat  exchanger  segment  with  cutaway  view  and  finite  slice 

for  finite  difference  method 


32 


Looking  at  a  finite  slice  of  the  segmented  heat  exchanger,  indicated  as  dx,  a 
graphical  representation  of  the  finite  difference  method  is  presented  in  Figure  13. 


Figure  13.  Finite  difference  method  representation  of  finite  slice  dx  of  concentric  tube  heat 

exchanger 


Heat  is  transferred  to  the  working  fluid  by  convection  from  the  inside  of  the  PDE 
tube  to  the  inner  tube  wall.  Heat  is  then  conducted  through  the  inner  tube  wall.  The  heat 
is  transferred  to  the  working  fluid  by  convection  from  the  tube  wall.  The  heat  is  then 
rejected  to  the  ambient  environment  or  transmitted  in  the  form  of  a  temperature  increase 
to  the  next  finite  slice  fluid  element.  The  rejected  heat  is  transferred  by  convection  to  the 
outer  tube  wall,  conducted  through  the  outer  tube  wall  and  rejected  to  the  ambient 
environment  by  free  convection  and  radiation. 

The  heat  transfer  equations  between  each  temperature  node  in  Figure  13  may  be 
represented  as 
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If  the  initial  inlet  temperature  is  used  for  the  first  finite  slice  then  there  are  ten  unknowns 
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Equations  (17)-(21)  represent  ten  equations  for  the  ten  unknowns  allowing  the 
determination  of  all  temperatures  and  heat  transfers  for  the  finite  slice.  The  temperature 
of  fluid  in  the  following  finite  slice  is  also  calculated  permitting  the  same  set  of  equations 
to  be  solved  for  as  each  subsequent  finite  slice. 


The  finite  difference  method  begins  with  the  inlet  temperature  as  the  fluid 
temperature.  The  next  fluid  element  temperature  is  calculated  as  an  unknown  and  the 
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cycle  begins  over  again  for  the  next  finite  slice.  This  continues  until  the  desired  outlet 
temperature  is  achieved.  The  number  of  cycles  required  multiplied  by  the  slice  size 
represents  the  length  of  the  heat  exchanger  required. 

Turbulent  flow  effects  on  the  heat  transfer  coefficient  of  the  working  fluid  are 
accounted  for  in  the  model,  as  are  material  and  fluid  property  variations  as  a  function  of 
temperature.  While  too  extensive  to  discuss  in  this  chapter  the  interested  reader  is 
referred  to  Appendix  B  for  the  final  calculations  for  the  method  used. 

Other  Design  Considerations 

To  match  the  current  PDE  tubes,  the  inner  tube  of  the  heat  exchanger  was  limited 
to  2”  stainless  steel  pipe.  Using  the  existing  2”  pipe  geometry  and  threading  on  the  heat 
exchanger  inner  tube  ensured  that  the  heat  exchanger  could  easily  be  moved  along  the 
axial  length  to  adjust  the  heat  transfer  to  the  fluid. 

The  annular  region  thickness  was  the  primary  consideration  in  determining  the 
outside  tube  diameter.  It  was  desired  to  have  the  smallest  annular  region  possible  to 
reduce  residence  time  and  minimize  the  amount  of  hot  fuel  in  the  heat  exchanger  while 
maximizing  the  surface  contact  area  with  the  heated  inner  tube.  The  outer  tube  selection 
was  limited  to  the  selection  of  the  pipe  stock  available  and  manufacturing  capabilities. 
With  these  criteria  in  mind  2.5”  stainless  steel  schedule  10  tube  was  selected  for  the  outer 
tube  which  provided  a  3.302  mm  annular  thickness. 

A  concern  in  the  early  design  process  was  that  boiling  may  occur  in  regions  prone 
to  stagnation.  To  prevent  stagnation,  a  helical  rod  was  incorporated  into  the  design.  It 
was  believed  that  installing  the  helical  rod  would  force  the  flow  in  a  helical  path  from  the 
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inlet  to  the  exit  preventing  any  hot  spots  and  cold  spots.  It  is  shown  in  the  chapter  four 
that  the  helical  rod  did  not  work  as  intended  due  to  buoyancy  in  the  working  fluid. 
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ITT.  Facilities  and  Instrumentation 


Pulsed  Detonation  Research  Facility 

This  research  was  conducted  at  the  Pulse  Detonation  Research  Facility  located  at 
Building  7 1  A,  D  Bay,  Wright-Patterson  Air  Force  Base,  Ohio.  The  Pulse  Detonation 
Research  Facility  is  a  component  of  the  Air  Force  Research  Laboratory  Propulsion 
Directorate  (AFRL/PR). 

This  facility  was  originally  designed  and  constructed  for  turbojet  testing,  and  has 
since  been  modified  to  support  pulse  detonation  research.  The  major  components  of  the 
facility  include  the  test  cell,  control  room,  and  fuel  room.  The  21,200  m  test  cell 
features  a  static  thrust  stand  capable  of  measuring  thrust  of  up  to  267,000  N  (Schauer, 
2001:3).  For  pulse  detonation  research,  the  thrust  values  and  the  size  of  the  research 
equipment  does  not  necessitate  the  use  of  the  high  capacity  static  thrust  stand.  A  smaller 
damped  thrust  stand  has  been  installed  above  the  static  thrust  stand.  It  is  on  the  small 
damped  thrust  stand  that  all  pulse  detonation  research  is  completed.  During  testing  the 
products  are  ejected  into  an  exhaust  tunnel  and  vented  to  the  atmosphere.  The  test  cell  is 
explosion  proof  and  is  separated  from  the  control  room  and  the  fuel  room  by  0.61  m  thick 
concrete  walls. 

The  control  room  is  located  adjacent  to  the  test  cell.  Personnel  are  not  permitted 
in  the  test  cell  during  testing.  To  monitor  the  test  cell,  a  number  of  closed  circuit  cameras 
are  placed  throughout  the  test  cell  and  fuel  room.  All  testing  is  controlled  by  use  of  a 
facility  specific  control  panel  and  LabView  control  software  run  from  a  control  computer. 
Real-time  data  is  recorded  by  the  LabView  control  software  run  from  the  control 
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computer,  while  high-speed  data  is  recorded  by  a  second  computer  also  running 
additional  LabView  software. 

The  fuel  room,  located  adjacent  to  the  control  room  and  the  test  cell,  provides  a 
low  point  ventilation  system,  liquid  fuel  storage,  and  fuel  conditioning  equipment.  All 
liquid  fuels  are  processed,  stored,  conditioned,  and  pressure  fed  from  the  fuel  room  to  the 
test  cell. 

Air  Supply  System 

Compressed  air  for  the  PDE  is  provided  by  one  of  two  Ingersoll-Rand  Pac  Air 
Compressors  (Model  No.  PA  300V).  Each  261  kW  compressor  provides  40  m  / min  at  a 
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rated  operating  pressure  of  689  kPa.  Air  from  the  compressors  is  stored  in  a  4.5  m 
collection  tank  fabricated  by  Buckeye  Fabrication  Co.  Inc.,  Dayton  Ohio  (Serial  Number 
10894).  From  the  collection  tank  the  air  is  piped  into  the  test  cell. 

The  main  supply  line  is  split  into  two  separate  lines  under  the  test  stand  for  the  fill 
and  purge  supply  lines,  indicated  in  Figure  14.  A  large  range  of  air  mass  flow  rates  are 
possible  by  the  use  of  fixed  diameter  orifice  plates  in  line  with  the  fill  and  purge  supply 
lines.  For  the  purpose  of  this  work  an  8.99  mm  diameter  orifice  plate  was  used  in  the 
purge  supply  line  and  a  12.7  mm  diameter  orifice  plate  was  used  in  the  fill  supply  line. 
During  testing,  air  mass  flow  rate  is  controlled  from  the  control  computer  by  Tescom 
Electropneumatic  PID  Controllers  (Model  No.  ER  1200)  which  actuate  dome  loader  type 
pressure  regulators  upstream  of  the  orifice  plates.  Surge  tanks  are  located  down  stream 
of  the  orifice  plates  to  prevent  compression  waves  generated  in  the  valve  system  from 
disrupting  the  flow  at  the  orifice  plates. 


38 


Under  Test  Stand 


Tescom  Pressure 

To  Purge  To  Air  Controller 


Figure  14.  Diagram  of  PDE  main  air  supply  split  to  fill  and  purge  air  supply  lines  with  associated 

hardware 

The  purge  air  enters  the  test  stand  and  is  piped  into  the  purge  manifold  for  use  in 
the  PDE.  The  fill  air  enters  the  test  stand  where  it  flows  through  a  15  kW  Chromalox 
Circulation  Heater  (P/N  053-500870-187).  The  air  temperature  can  be  adjusted  by  setting 
the  upper  temperature  limit  on  the  Chromalox  temperature  controller  (Model  No.  2104) 
on  the  facility  control  panel  or  by  adjusting  a  setting  on  the  control  computer.  After 
exiting  the  heater,  the  fill  air  is  mixed  with  the  fuel  and  enters  the  fill  manifold  for  use  in 
the  PDE. 

Hydrogen  Fuel  Supply  System 

The  hydrogen  fuel  supply  is  provided  by  a  hydrogen  tuber  trailer  located  outside 
of  the  research  facility.  The  trailer  holds  38  individual  hydrogen  bottles  each  with  a 
starting  pressure  of  16.55  MPa.  The  total  water  volume  of  the  trailer  is  15.69  m  .  During 
testing,  individuals  bottles  are  used  until  the  bottle  pressure  falls  below  the  pressure 
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required  for  a  given  test  condition,  at  which  time  a  new  bottle  is  opened.  This  continues 
until  the  entire  hydrogen  trailer  has  been  emptied  at  which  time  the  trailer  is  replaced. 

The  hydrogen  enters  the  research  facility  through  the  test  cell.  The  hydrogen  is 
routed  under  the  test  stand  through  a  series  of  manual  and  pneumatic  ball  valves.  The 
pneumatic  ball  valves  provide  the  ability  to  turn  the  hydrogen  flow  on  or  off  during 
testing  from  the  facility  control  computer.  The  hydrogen  mass  flow  rate  is  determined  by 
a  Flow-Dyne  critical  flow  nozzle  (P/N  N08006-SA).  During  testing,  the  hydrogen  mass 
flow  rate  is  regulated  from  the  control  computer  by  a  Tescom  Electronic  Pressure 
Controller  (Model  No.  ER3000)  and  TESCOM  Pressure  Regulator  (Model  No.  44-1313- 
2 122A200)  upstream  of  the  critical  flow  nozzle.  A  surge  tank  is  located  downstream  of 
the  critical  flow  nozzle  to  prevent  shock  waves  generated  in  the  valve  system  from 
disrupting  the  flow  at  the  critical  flow  nozzle. 

For  this  research  the  hydrogen  was  injected  prior  to  the  inlet  manifold  by  means 
of  a  T-section  in  the  air  supply  line  located  at  the  exit  of  the  air  heater  2.5  m  upstream  of 
the  injection  valve.  Due  to  the  gaseous  state  of  the  hydrogen,  nozzles  or  mixing  devices 
were  not  needed  in  the  referenced  tests.  Instead  the  turbulent  nature  of  the  airflow  and 
fuel  injection,  in  addition  to  the  flow  path  are  sufficient  to  ensure  a  near  homogenous 
mixture. 

Liquid  Fuel  Supply  System 

The  two  liquid  hydrocarbon  fuels  used  for  this  work  were  aviation  gasoline  and 
JP-8  military  turbine  fuel,  hereto  referred  to  as  avgas  and  JP-8  respectively.  Avgas  was 
used  due  to  its  availability,  detonability,  and  numerous  similarities  to  JP-8,  which  is  the 
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primary  fuel  of  interest  for  this  work.  The  similarity  in  the  fuel  system  for  these  fuels  is 
near  identical  and  hence  will  be  presented  in  such  a  fashion  with  any  pertinent  differences 
noted.  For  JP-8  flash  vaporization  and  supercritical  heating  tests  the  JP-8  requires 
conditioning  which  is  discussed  in  the  following  section.  This  section  will  explain  the 
fuel  system  as  if  the  fuel  has  already  been  conditioned  or  does  not  require  conditioning. 

Avgas  is  directly  available  from  the  fuel  farm  adjacent  to  the  facility  and  is 
pumped  directly  from  the  fuel  fann  to  the  fuel  room.  JP-8  is  obtained  locally  from 
AFRL/PR  Fuels  Branch  and  stored  in  18.93  L  fuel  containers  until  use. 

For  testing,  the  fuel  was  placed  in  a  41.64  L  stainless  steel  general  purpose 
pressure  vessel  (S/N  28108-007).  The  fuel  is  pressure  fed  by  a  standard  nitrogen  bottle 
to  two  9.46  L  Greer  hydraulic  accumulators  (Model  No.30A-2/4A)  rated  to  20.86  MPa. 
The  accumulators  are  pressurized  by  a  standard  nitrogen  bottle  and  high  pressure 
regulator.  The  accumulators  use  a  rubber  diaphragm  that  to  separate  the  nitrogen  source 
from  the  fuels  and  is  compatible  with  the  fuels  used  in  this  work  (Tucker,  2005).  The 
fuel  fill  and  pressurization  schematic  are  present  in  Figure  15. 
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Figure  15.  Diagram  of  fuel  room  configuration  and  process  of  filling  the  accumulators  and  also 
pressuring  the  accumulators  providing  pressurized  fuel  for  testing 

The  fuel  line  out  of  the  fuel  room  is  opened  by  manual  ball  valves  allowing  a  flow 
path  to  a  point  underneath  the  test  stand.  At  this  point  two  manual  ball  valves  control  one 
of  two  flow  paths.  The  first  is  a  flow  meter  bypass  used  for  system  pressurization.  The 
second  path  is  through  a  Flow  Technology  turbine  volumetric  flow  meter  (Model  No. 
FT4-8AEU2-LEAT5).  A  thermocouple  immediately  after  the  flow  meter  permits 
temperature  compensation  in  fuel  density  for  fuel  mass  flow  to  be  calculated  by  the 
control  computer.  From  the  flow  meter,  the  fuel  line  traverses  the  underside  of  the  test 
stand  to  a  pneumatic  valve  actuated  from  the  control  computer.  This  valve  is  referred  to 
as  the  last  chance  valve  and  controls  the  flow  of  fuel  during  test.  After  the  last  chance 
valve  the  fuel  line  enters  the  test  stand.  The  diagram  of  the  flow  meter  and  flow  meter 
bypass  are  presented  in  Figure  16. 
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Figure  16.  Diagram  of  fuel  flow  meter  and  flow  meter  bypass  system  with  last  chance  valve 
For  testing  without  fuel  heating,  the  fuel  line  was  routed  to  the  fuel  inlet  manifold. 
For  fuel  heating  tests  the  fuel  entered  the  flash  vaporization  system,  discussed  in  a 
following  section,  and  terminated  in  the  fuel  inlet  manifold.  The  fuel  inlet  manifold 
consists  of  two  spray  bars  that  have  been  welded  into  the  fill  air  manifold  as  shown  in 
Figure  17.  Each  spray  bar  is  drilled  and  tapped  with  five  threaded  holes  that  accept  the 
Delavan  fuel  flow  nozzles.  The  number  of  nozzles  and  flow  number  of  the  nozzles  used 
varied  from  test  to  test  depending  on  the  fuel,  fuel  pressure,  and  mass  flow  requirements. 


Figure  17.  Upstream  view  into  the  fill  air  manifold  of  the  fuel  inlet  manifold  with  spray  bars 

Delavan  fuel  flow  nozzles 
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Fuel  Conditioning 

For  heated  JP-8  tests,  fuel  conditioning  is  required  to  minimize  carbon  deposits  in 
the  flash  vaporization  system  due  to  the  dissolved  oxygen  reacting  with  the  fuel.  Fuel 
condition  begins  with  placing  the  JP-8  in  a  41.64  L  general  purpose  pressure  vessel.  A 
coiled  piece  of  stainless  steel  tubing  in  the  bottom  of  the  vessel  with  numerous  small 
diameter  holes  drilled  in  it  is  attached  to  a  piece  of  tubing  that  exits  the  pressure  vessel 
and  is  connected  to  a  standard  nitrogen  bottle  and  regulator.  A  ball  valve  on  the  top  of 
the  pressure  vessel  is  vents  the  tank  to  the  outside  atmosphere.  The  nitrogen  pressure  is 
raised  until  bubbling  through  the  coils  is  audibly  detected.  After  the  proper  volume  of 
nitrogen  (Tucker,  2005)  has  been  bubbled  through  the  JP-8,  the  vent  to  atmosphere  is 
closed  and  the  holding  tank  is  pressurized.  At  this  point  the  JP-8  may  be  used  to  fill  the 
fuel  system  as  mentioned  in  the  previous  section. 


Figure  18.  Top  view  of  fuel  conditioning  holding  tank  with  nitrogen  bubbling  coiled  tube  at  the  tank 

bottom 

Ignition  System 

The  ignition  system  timing  is  determined  by  the  angular  position  of  the  camshaft 
in  the  PDE.  The  camshaft  position  is  read  by  a  BEI  optical  encoder  (Model  H25,  S/N 
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YOO 13039)  and  sent  to  the  control  computer.  An  ignition  delay  may  be  set  in  the  control 
computer  by  the  user  and  is  often  required  to  prevent  backfiring.  With  the  encoder  input 
and  the  user  defined  ignition  delay  the  control  computer  will  detennines  the  ignition 
timing  and  sends  a  signal  to  the  ignition  relay  box.  The  relay  box  sends  the  signal  to  12 
VDC  MSD  Digital  DIS-4  Ignition  system.  The  ignition  system  releases  105-1 15  mJ  per 
spark.  For  each  fire  cycle  the  control  computer  commands  four  sparks  signals  per  tube 
for  total  ignition  energy  of  420-460  mJ  per  tube  per  cycle. 

The  spark  plugs  used  are  a  modified  NGK  plug  that  has  the  grounding  electrode 
removed  and  a  small  piece  of  tube  welded  to  the  end  of  the  spark  plug.  The  cylinder 
provides  the  grounding  capability  for  the  spark  and  a  small  non-turbulent  environment  for 
explosive  reaction  to  develop. 

Pulse  Detonation  Engine 

The  PDE  used  for  this  research  uses  the  head  from  a  General  Motors  Quad  4  head 
with  dual  overhead  camshafts  shown  in  Figure  19.  The  camshafts  are  belt  driven  by  a 
variable  speed  14.9  kW  Baldor  electrical  motor  (Model  No.  M4102T,  S/N 
C0303240142).  Motor  control  and  frequency  are  provided  by  the  control  computer. 
Instead  of  the  intake  and  exhaust  cycles  typical  in  automotive  engines,  the  valve  system 
in  the  PDE  corresponds  to  fill  and  purge  cycles  respectively.  On  the  fill  cycle,  the  intake 
valves  open  allowing  a  fuel-air  mixture  to  flow  into  the  PDE.  Likewise,  on  the  purge 
cycle  the  exhaust  valves  open  allowing  purge  air  to  flow  into  the  PDE  engine.  The  valve 
train  is  lubricated  by  automotive  oil  which  is  pumped  into  the  head  from  an  external 
reservoir  and  electric  oil  pump.  The  cooling  water  is  provided  by  a  closed  loop  radiator 
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cooling  system,  which  is  pumped  through  the  existing  head  cooling  water  ports  by  use  of 
an  external  electric  water  pump. 


Figure  19.  GM  Quad  4  head  being  used  as  PDE  valve  train  for  fill  air  manifold  (top)  and  purge 

manifold  (bottom) 

PDE  tube  mounting  plates  are  attached  to  the  head  with  the  existing  head  bolts. 
The  mounting  plates  are  1 .27  cm  stainless  steel  plate  that  mount  flush  with  the  head  and 
are  sealed  by  the  use  of  a  stock  head  gasket.  In  general,  the  mounting  plates  are  threaded 
to  accept  male  2”  national  pipe  thread  (NPT). 

When  using  liquid  hydrocarbon  fuels  it  is  necessary  to  use  detonation  initiating 
devices  to  achieve  detonations  in  reasonable  lengths  of  tube.  Numerous  types  of 
detonation  tripping  devices  exist,  however  for  this  research  a  structurally  enhanced 
shchelkin  spiral  (Shchelkin,  1940)  was  and  is  shown  in  Figure  20.  The  spiral  is  installed 
prior  to  the  mounting  plates  and  held  in  place  once  by  the  mounting  plates.  The  PDE 
tube  is  then  slid  over  the  spiral  and  screwed  into  the  mounting  plate. 
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Figure  20.  Shelkin  Like  Spiral  with  Structural  Support 

The  PDE  tubes  used  for  this  research  were  1.829  m  long  and  fabricated  from  2” 
stainless  steel  schedule  40  pipe  sections  joined  with  2”  NPT  collars.  During  tests  with  a 
heat  exchanger  installed,  various  pipe  sections  were  used  to  allow  axial  positioning  of  the 
heat  exchanger  while  maintaining  the  tube  length  of  1.829  m. 

Heat  Exchanger  Configuration 

For  this  research  two  heat  exchangers  were  developed.  The  first  heat  exchanger 
used  was  designed  based  on  assumptions  and  heat  transfer  calculations  discussed  in  the 
previous  chapter  and  presented  in  Appendix  B.  After  multiple  tests  it  was  realized  that  a 
smaller  second  heat  exchanger  was  required  due  to  the  high  heat  transfers  and  to  validate 
theories  developed  using  the  first  heat  exchanger.  The  first  heat  exchanger  was  76.2  cm 
long  and  the  second  heat  exchanger  was  38.1  cm  long;  and  therefore  the  heat  exchangers 
shall  be  referred  to  as  the  long  heat  exchanger  and  the  short  heat  exchanger  for  the 
remainder  of  this  document. 

The  long  heat  exchanger  is  a  76.2  cm  long  concentric  tube  configuration  with  the 
inner  tube  constructed  from  2”  stainless  steel  schedule  40  pipe  and  the  outer  tube 
constructed  from  2  1/2”  stainless  steel  schedule  10  pipe.  It  would  seem  counter-intuitive 
that  a  thinner  walled  pipe  would  be  used  on  the  outside  given  the  high  pressures  in  the 
heat  exchanger,  however  the  annular  region  in  between  the  two  tubes  was  the  critical 
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concern  in  the  design  as  long  as  the  heat  exchanger  could  pass  structural  design 
considerations. 

The  two  pipes  were  welded  to  6.35  mm  stainless  steel  end  plates  allowing  equal 
annular  spacing  of  3.302  mm  around  the  perimeter.  In  an  attempt  to  prevent  stagnation 
of  the  fluids  within  the  heat  exchanger  a  3.175  mm  diameter  stainless  steel  rod  was  spot 
welded  to  the  inner  tube  in  a  helical  path  with  six  full  turns.  The  helical  rod  is  shown 
welded  to  the  inner  tube  in  Figure  2 1 .  The  helical  rod  did  not  ensure  a  fluid  tight  flow 
barrier,  but  it  was  hypothesized  that  the  helix  would  force  a  helical  flow  path  and  a  more 
even  temperature  distribution  throughout  the  length  of  the  heat  exchanger  flow  path 
would  be  developed.  The  outer  tube  had  two  1/4”  Swagelok  unions  welded  at  opposite 
ends  and  180  degree  radially  offset.  These  fittings  corresponded  to  the  fluid  inlet  and 
outlet  positions. 


Figure  21.  Construction  of  the  long  heat  exchanger  with  helical  rod  welded  in  place 

To  mount  the  heat  exchanger  to  existing  PDE  tubes,  two  15.24  cm  extensions 
were  fabricated  from  2”  stainless  steel  schedule  40  pipe.  At  one  end  of  the  extension  a 
6.35  mm  stainless  steel  endplate  for  mounting  to  the  heat  exchanger  is  welded  in  place. 
The  other  end  of  the  extension  is  threaded  with  male  2”  NPT  mate  for  mounting  to  the 
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PDE  head  or  to  other  PDE  tube  sections  with  female  2”  NPT  pipe  collars  as  shown  in 
Figure  22.  In  each  extension  two  3/8”-24  stainless  steel  bolts  are  welded  three  inches 
apart  from  one  another  for  wave  speed  measurement. 


Figure  22.  Heat  exchanger  connecting  extension  with  end  plate  for  heat  exchanger  installation, 
instrumentation  ports,  and  male  2”  NPT  connected  to  female  2”  pipe  collar 

The  short  heat  exchanger  was  also  a  concentric  tube  configuration.  The  inner 
tube  is  fabricated  of  2”  stainless  steel  schedule  10  pipe  and  the  outer  tube  is  2.5”  stainless 
steel  schedule  10  pipe.  The  short  heat  exchanger  was  38.1  cm  long.  The  short  heat 
exchanger  did  not  have  the  helical  flow  coil  installed.  The  same  size  inlet  and  outlet 
fittings  were  used  and  installed  at  opposite  ends  of  the  tube  and  180  degrees  radially 
offset.  1/8”  swagelok  fittings  were  installed  at  the  12.7  cm  and  25.4  cm  axial  positions 
and  radially  aligned  with  the  outlet.  These  fittings  were  used  to  insert  thermocouples  into 
the  flow.  The  short  heat  exchanger  utilized  the  same  endplate  construction  and 
connecting  extensions  used  in  the  long  heat  exchanger.  The  short  heat  exchanger  is 
presented  in  Figure  23. 
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Figure  23.  Profile  view  of  short  heat  exchanger  with  inlet  and  outlet  ports  at  opposing  ends  and  two 
spaced  thermocouple  flow  ports  on  the  outlet  side 

Both  heat  exchangers  were  hydrostatically  pressure  tested  in  accordance  with 
ASME  B3 1 .3.  The  rated  working  pressure  for  both  heat  exchangers  is  5.5 16  MPa  at 
588  K. 

Water  Flash  Vaporization  System 

Safety  requirements  necessitated  that  initial  tests  be  completed  using  water  for 
flash  vaporization.  This  served  as  a  proof-of-concept  for  a  JP-8  flash  vaporization  system 
(FVS).  These  tests  also  provided  insight  into  the  heat  transfer  rates  and  flow 
characteristics  of  the  FVS  and  heat  exchangers. 

The  water  FVS  consists  of  a  standard  nitrogen  bottle,  water  reservoir,  one  of  the 
two  previously  mentioned  heat  exchangers,  high  temperature/pressure  ball  valve, 
thermocouples,  pressure  transducer,  and  all  associated  tubing  and  fittings  required  to 
properly  connect  all  components.  All  components  in  the  FVS  are  connected  by  %” 
stainless  steel  tubing  connected  with  varying  types  of  14”  stainless  steel  Swagelok 
compression  fittings. 

Reverse  osmosis  (RO)  water  was  used  in  all  tests  using  water  as  the  heated  fluid. 
RO  water  was  obtained  locally  and  transported  by  18.93  L  container.  The  water  reservoir 
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is  a  locally  manufactured  2.134  m  long,  2”  stainless  steel  schedule  80  pipe  with  end  caps 
welded  in  place.  The  reservoir  provided  3.785  L  capacity  and  is  rated  to  a  working 
pressure  of  10.342  MPa.  The  filled  reservoir  is  pressurized  by  means  of  a  standard 
nitrogen  bottle  and  high  pressure  regulator.  The  water  is  pressurized  up  to  6.895  MPa 
and  pressure  fed  to  the  heat  exchanger  inlet. 

An  Omega  13.79  MPa  pressure  transducer  (Model  No.  PX02K1-2KG5T)  at  the 
heat  exchanger  outlet  monitors  static  pressure.  The  water  flows  from  the  heat  exchanger 
to  a  pneumatically  operated  Bonetti  stainless  steel  ball  valve  (Model  No.  941-NPT-600- 
3 1)  rated  to  4.72  MPa  at  823  K.  The  pneumatic  valve  provided  water  flow  control  during 
test  from  the  control  computer.  A  diagram  of  the  water  FVS  is  presented  in  Figure  24. 


Figure  24.  Diagram  of  PDE  engine  with  water  FVS  and  instrumentation  installed 

After  passing  through  the  hot  valve  the  heated  water  is  exhausted  to  the  test  cell 
atmosphere  by  use  of  a  spray  bar  using  the  same  type  of  Delavan  spray  nozzles  used  in 
the  fuel  system.  A  picture  of  the  spray  bar  with  nozzles  installed  is  present  in  Figure  25 


51 


The  spray  bar  was  manufactured  out  of  19.05  mm  stainless  steel  hexagon  bar  stock.  The 
spray  bar  was  drilled  and  ported  for  five  spray  nozzles.  A  1/4”  Swagelok  union  was 
welded  to  provide  an  inlet  for  the  hot  water  supply  from  the  FVS.  The  exhausted  heated 
water  is  vaporized  after  crossing  the  nozzles  and  experiencing  the  pressure  drop  to 
ambient  atmosphere. 


Figure  25.  Water  spray  bar  with  Delavan  spray  nozzles  installed 
Water  mass  flow  was  measured  by  exhausting  the  water  into  an  18.927  L  stainless 
steel  general  purpose  pressure  vessel  (S/N  109702-024).  The  water  vapor  was  ejected 
into  the  vessel  which  was  partially  filled  with  water  to  a  level  above  the  spray  nozzles. 
The  vessel  had  a  small  fitting  less  than  one  centimeter  that  vented  the  vessel  to 
atmosphere.  The  mass  losses  to  the  environment  due  to  steam  leaving  the  vessel  were 
assumed  to  be  negligible.  The  vessel  hung  from  a  1334  N  capacity  BLH  load  cell  (Type 
U3XXA,  S/N  36557).  Load  cell  data  is  sent  to  the  control  computer. 


JP-8  Flash  Vaporization  System 

The  JP-8  FVS  consists  of  a  pneumatic  fuel  control,  a  nitrogen  bottle  for  purging 
the  hot  fuel  from  the  system,  a  pneumatic  control  for  nitrogen  purge,  a  check  valve  for 
the  nitrogen  purge,  one  of  the  two  previously  mentioned  heat  exchangers,  thermocouples, 
pressure  transducer,  and  all  associated  tubing  and  fittings  required  to  properly  connect  the 
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components.  All  components  in  the  JP-8  FVS  are  connected  by  stainless  steel  tubing 
fastened  with  varying  types  of  stainless  steel  Swagelok  compression  fittings,  unless 
otherwise  noted.  The  purpose  of  the  nitrogen  purge  is  to  maintain  FVS  pressurization 
and  purge  the  remaining  JP-8  in  the  FVS  after  the  fuel  supply  has  been  shut  off. 

The  JP-8  FVS  begins  at  the  last  chance  valve  under  the  test  stand.  Immediately 
after  the  last  chance  valve  a  tee  connection  introduces  the  nitrogen  purge  into  the  system. 
The  nitrogen  purge  is  controlled  by  a  pneumatic  ball  valve  and  inline  check  valve.  The 
nitrogen  purge  pneumatic  valve  controls  the  flow  of  the  nitrogen  and  the  check  valve 
ensures  that  the  nitrogen  purge  will  initiate  if  the  fuel  pressure  drops  below  the  nitrogen 
pressure  set  point.  For  the  purpose  of  this  work  the  nitrogen  purge  was  set  values  above 
the  critical  pressure  of  JP-8. 

The  flow  path  enters  the  test  stand  through  a  manual  ball  valve  where  the  flow  is 
routed  to  the  heat  exchanger  inlet.  An  Omega  pressure  transducer  rated  to  13.79  MPa  is 
placed  at  the  heat  exchanger  outlet  to  monitor  static  pressure.  From  the  heat  exchanger 
outlet  the  flow  traveled  through  an  insulated  fuel  line  which  is  terminated  in  the  fuel  inlet 
manifold.  The  JP-8  FVS  diagram  is  presented  in  Figure  26. 
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Figure  26.  Diagram  of  PDE  engine  with  JP-8  FVS  and  instrumentation  installed 


Temperature  Instrumentation 

Inlet  and  outlet  heat  exchanger  flow  temperatures  were  measured  using  J-Type 
thermocouples  placed  at  the  centerline  supply  lines.  The  measurements  were  taken  7.62 
cm  upstream  and  downstream  of  the  heat  exchanger  entrance  and  exhaust  respectively. 
Internal  heat  exchanger  flow  temperatures  were  measured  by  1/8”  J-Type  thermocouples 
placed  in  the  middle  of  the  annular  region  of  the  heat  exchanger. 


During  JP-8  heated  fuel  tests  the  fuel  manifold  inlet  temperatures  were  read  by  a 
J-Type  thermocouple  placed  at  the  centerline  of  the  supply  line  7.62  cm  upstream  of  the 
fuel  manifold  inlet.  The  fuel  line  was  insulated  during  these  tests  and  the  temperature 
losses  from  the  thennocouple  reading  to  fuel  nozzles  was  assumed  to  be  negligible. 
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Air  flow  temperature  measurements  were  taken  by  J-type  thermocouples  0.5  m 
upstream  of  the  fuel  inlet  manifold.  Fuel-air  mixture  temperature  was  recorded  by  a  J- 
type  thermocouple  0.25  m  downstream  of  the  fuel  inlet  manifold. 

External  PDE  tube  wall  temperatures  were  measured  by  J-type  thermocouples 
mounted  externally  by  compression  clamps  to  the  PDE  tube.  Figure  27  shows  the  use  of 
external  thermocouples  in  addition  to  flow,  inlet,  and  outlet  thennocouples  installed  on 
the  short  heat  exchanger. 


Figure  27.  Short  heat  exchanger  installed  with  surface,  flow,  inlet,  and  outlet  thermocouples 
Facility  Control  Software 

The  Labview  PDE  Facility  Control  and  Data  Acquisition  software  provides 
realtime  control  and  data  acquisition  of  all  PDE  operating  parameters.  The  software 
controls  engine  frequency,  purge  and  fill  air  mass  flow  rates,  purge  and  fill  fraction, 
ignition  delay,  spark  control,  and  fuel  flow  control,  just  to  name  a  few.  The  software 
provides  real-time  monitoring  of  all  of  the  control  parameters  in  addition  a  multitude  of 
air/fuel  pressure  and  temperatures. 
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Test  Configuration  for  PDE  Tube  Tests  without  Heat  Exchanger 


Plain  tube  tests  were  completed  with  avgas-air  and  JP-8-air  detonations  to 
determine  steady-state  external  wall  temperature,  heat  transfer  rates,  and  wavespeed  data 
as  functions  of  axial  position.  Two  1.829  m  PDE  tubes  with  1.219  m  shelkin  like  spirals 
were  used  for  these  tests.  For  the  avgas  test,  fill  air  was  heated  to  32 1  K  and  the 
equivalence  ratio  was  varied  between  0.95-1.26.  For  the  JP-8  test,  fill  air  was  heated  to 
394  K  and  the  equivalence  ratio  was  varied  between  1.05-1.26.  Engine  operating 
frequency  was  maintained  at  15  Hz  and  an  8  ms  ignition  delay  was  used.  Ion  probes  and 
external  thermocouples  locations  are  shown  in  Table  4. 


Table  4.  Plain  tube  test  ion  probe  and  thermocouple  locations 


Ion  Probe  and  External 

Axial  PDE  Tube  Location 

Therm  couple 

(cm) 

1 

19.69 

2 

46.99 

3 

73.82 

4 

105.89 

5 

121.76 

6 

137.01 

7 

154.62 

8 

169.86 

Water  FVS  with  Hydrogen-Air  Detonation  Configuration 

Water  FVS  tests  with  hydrogen-air  detonations  were  completed  with  long  heat 
exchanger.  For  hydrogen  fueled  tests  one  1.829  m  PDE  tube,  with  heat  exchanger,  was 
used.  The  spiral  for  the  hydrogen  fueled  tests  was  .305  m  long.  The  hydrogen-air 
mixture  was  maintained  at  an  equivalence  ratio  of  1.0  with  an  engine  frequency  of  10  Hz, 
and  ignition  delay  of  6  ms. 
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Water  FVS  with  Avgas- Air  Detonation  Configuration 

Water  FVS  tests  with  avgas-air  detonations  were  completed  with  long  and  short 
heat  exchanger.  For  avgas  fueled  tests  two  1.829  m  PDE  tubes  were  used,  each  with  a 
1.219  m  spiral.  Only  one  heat  exchanger  was  incorporated  for  all  tests. 

JP-8  FVS  Configuration 

JP-8  FVS  tests  were  completed  with  long  and  short  heat  exchanger.  Two  1.829  m 
PDE  tubes  were  used,  each  with  a  1 .2 19  m  spiral.  Only  one  heat  exchanger  was 
incorporated  for  all  tests.  A  photograph  representing  a  two-tube  test  configuration  with 
the  long  heat  exchanger  installed  at  the  hot  section  is  shown  in  Figure  28.  Inlet  and  outlet 
thermocouples  can  be  seen  in  addition  to  PDE  tube  ion  probes. 


Figure  28.  Generic  two-tube  configuration  with  instrumented  long  heat  exchanger  installed  on 
closest  PDE  tube  with  the  inlet  at  the  end  of  the  tube  and  the  outlet  toward  the  front  of  the  tube 


57 


IV.  Results  and  Analysis 


Determination  of  Water  Mass  Flow  Rate 

Water  mass  flow  rate  was  determined  by  ejecting  the  heated  water  into  a  holding 
tank  that  was  hung  from  a  load  cell.  The  slope  of  the  load  cell  recordings  with  time 
provided  the  mass  flow  rate  for  water  tests.  A  sample  image  of  the  load  cell  recordings 
and  how  mass  flow  rate  is  calculated  is  shown  in  Figure  29. 


Load  Cell  vs.  Time 
(Representative  Sample) 


Figure  29.  Sample  mass  flow  calculation  based  on  slope  of  load  cell  versus  time 


Heat  Transfer  Calculation 

For  heat  exchanger  tests,  the  heat  transfer  rate  to  the  working  fluid  was 
determined  from  the  mass  flow,  the  specific  heat  of  the  working  fluid,  and  the  inlet/outlet 
temperatures.  Mass  flow  was  determined  by  the  method  in  the  previous  paragraph  for 
water  FVS  tests  and  using  flow  meter  data  for  JP-8  FVS  tests.  If  specific  heat  is  assumed 
to  be  constant  the  determination  of  heat  transfer  rate  is  governed  by 

C1  •=  n\lot'cp' (^outlet  ~  ^ inlet)  (44) 
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Due  to  the  elevated  temperatures  achieved  in  the  heat  exchanger,  specific  heat 
may  no  longer  be  assumed  constant  and  must  be  accounted  for  as  a  function  of 
temperature.  An  approximation  of  specific  heat  was  obtained  by  using  the  specific  heat 
value  for  the  average  temperature  between  the  inlet  and  outlet  temperature.  This  specific 
heat  value  is  obtained  by  linearly  interpolating  the  average  temperature  with  the  fluid 
property  values.  The  calculations  for  the  heat  exchanger  tests  are  presented  in 
Appendix  E. 

During  the  plain  tube  tests  the  measured  external  tube  wall  temperature  and 
ambient  temperature  were  used  to  calculate  the  heat  transfer  rate  from  the  tube  to  ambient 
environment  by  free  convection  and  radiation  calculation.  The  heat  transfer  rate  changes 
with  axial  tube  length  due  to  the  varying  temperature  profile;  therefore,  heat  transfer  rate 
is  calculated  as  a  function  of  length.  The  method  used  to  calculate  the  heat  transfer  rate  is 
presented  in  Appendix  C. 

Wave  Speed  Calculation 

Wave  speed  is  determined  by  ion  probe  and  time  data  collected  from  the  high 
speed  data  computer.  The  ion  probes  are  continuously  supplied  with  4.5  volts.  When  a 
wave  travels  past  the  ion  probe,  a  circuit  is  completed  causing  a  steep  decrease  in  the 
recorded  voltage  as  shown  in  Figure  30  and  Figure  3 1 .  The  energizing  system  quickly 
recharges  the  electrode  before  the  next  ignition  event  occurs. 

The  time  the  wave  passes  a  probe  is  determined  by  the  threshold  value  below  the 
source  voltage  when  the  circuit  is  completed.  A  threshold  value  of  4.45  volts  was  used. 
This  value  was  chosen  because  it  is  outside  of  the  noise  level  of  the  voltage  reading  but 
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within  the  range  to  obtain  weak  wave  events  (See  Figure  31).  The  voltage  drop  for  a 
measurable  wave  event  is  a  nearly  linear  function  allowing  for  linear  interpolation  of  data 
points.  The  first  data  sample  after  the  threshold  value  and  the  point  prior  to  the  threshold 
value  are  used  to  interpolate  the  time  of  the  wave  event.  The  analyzed  time  data  has  a 
data  resolution  of  50  microseconds.  The  code  used  to  determine  the  wave  speed  is 
presented  in  Appendix  D. 


High  Speed  Data  vs.  Time 
(Representative  Sample) 


142.50  147.50  152.50  157.50  162.50 

Time  (milliseconds) 

Figure  30.  Sample  high  speed  data  with  spark  trace  shown  as  the  square  wave  and  the  ion  probe 

drop  due  to  the  wave  passing  the  sensors 


- Spark  Trace 

- 1  st  Ion  Probe 

- 2nd  Ion  Probe 

- 3rd  Ion  Probe 


High  Speed  Data  vs.  Time 
(Representative  Sample) 


- 1st  Ion  Probe 

- 2nd  Ion  Probe 

- 3rd  Ion  Probe 

- Ion  Wave  Speed  (WS) 

Threshold 


Figure  31.  Detail  of  ion  probe  voltage  as  wave  passes  including  the  wave  speed  threshold  about 

which  the  time  is  interpolated 
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PDE  Tube  Tests  without  Heat  Exchanger 

The  steady-state  surface  temperature  results  for  the  plain  tube  avgas-air 
detonation  test  are  shown  in  Figure  32.  Maximum  steady-state  temperatures  of  850  K  to 
875  K  were  seen  at  the  end  of  the  spiral  in  the  120  cm  -160  cm  axial  positions.  There  is  a 
linear  temperature  rise  from  the  head  of  the  PDE  to  the  end  of  the  spiral.  The  noticeable 
temperature  drop  at  the  end  of  the  tube  is  due  to  the  entrainment  of  cool  air  into  the  end 
of  tube  when  an  expansion  wave  is  created  after  the  wave  and  exhaust  products  have 
been  ejected  from  the  tube. 

Equivalence  ratios  of  near  1 . 1  provided  the  highest  steady-state  temperatures  in 
the  hot  section.  The  richer  mixtures  showed  over  a  50  K  drop  at  the  hottest  portion  of  the 
tube.  Near  stoichiometric  mixtures  provided  highest  temperatures  forward  of  the  hot 
section  and  temperatures  near  the  maximum  values  in  the  hot  section. 


Figure  32.  Plain  tube  surface  temperature  vs.  axial  distance  w/  varying  equivalence  ratio 
avgas  (298  K)  -  air  (322  K),  15  Hz,  8  ms  ignition  delay, 

1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
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Wave  speed  profiles  for  the  avgas  test,  shown  in  Figure  33,  closely  matched  the 
temperature  profile.  The  maximum  wave  speeds  were  obtained  at  the  110  cm  -  130  cm 
region  for  equivalence  ratios  just  above  stoichiometric.  It  is  at  these  locations  that  the 
wave  speeds  exceed  the  CJ  detonation  velocities.  Both  lean  and  rich  mixtures 
experienced  velocity  increases  with  length  but  not  to  the  magnitude  of  equivalence  ratios 
near  1.1.  The  richest  mixture  shows  an  early  velocity  spike  followed  by  a  rapid  decline 
to  deflagration  wave  speeds.  All  other  mixtures  achieve  CJ  velocities  by  the  end  of  the 
PDE  tube. 


Figure  33.  Plain  tube  wave  speed  vs.  axial  distance  w/  varying  equivalence  ratio 
avgas  (298  K)  -  air  (322  K),  15  Hz,  8  ms  ignition  delay, 

1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 

The  steady-state  surface  temperature  results  for  the  plain  tube  JP-8  -  air 
detonation  test  are  shown  in  Figure  34.  JP-8  temperature  and  velocity  profiles  closely 
mimic  those  of  the  avgas  tests  for  the  equivalence  ratios  tested.  The  JP-8  did  not  provide 
reliable  data  for  equivalence  ratios  that  were  lean  or  stoichiometric.  JP-8  temperature 
profiles  showed  tighter  grouping  throughout  the  range  of  equivalence  ratios.  As  was  seen 
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in  the  avgas  tests  the  maximum  temperatures  are  seen  in  the  120  cm  -160  cm.  The  JP-8 
achieved  from  25-50  K  higher  temperatures  in  the  hot  section  compared  to  the  avgas  test. 


Figure  34.  Plain  tube  surface  temperature  vs.  axial  distance  w/  varying  equivalence  ratio 
JP-8  (298  K)  -  air  (395  K),  15  Hz,  8  ms  ignition  delay, 

1.829  m,  2"  SS  Sch  40  Tube  w/  1.219  m  spiral 

The  JP-8  wave  speed  data,  in  Figure  35,  showed  a  tighter  wave  speed  grouping 
than  avgas.  Unlike  avgas,  JP-8  had  maximum  wave  speeds  at  the  richest  equivalence 
ratios  in  the  1 .2  range  for  unheated  fuel.  The  transient  overdriven  wave  speeds  can  be 
witnessed  at  the  end  of  the  spiral  with  a  decrease  to  an  equilibrium  CJ  point  towards  the 
end  of  the  tube  for  all  mixtures. 
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Figure  35.  Plain  tube  wave  speed  vs.  axial  distance  w/ varying  equivalence  ratio 
JP-8  (298  K)  -  air  (395  K),  15  Hz,  8  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 

The  heat  transfer  rate  to  the  surroundings  for  the  avgas  test  at  the  hot  section  was 
1.8-2. 5  kW/m;  while  JP-8  had  values  of  2. 4-2. 9  kW/rn,  shown  in  Figure  36  and  Figure  37 
respectively.  The  difference  between  the  two  tests  is  attributed  to  the  higher  maximum 
tube  temperatures.  The  heat  transfer  rates  are  below  the  required  heat  transfer  rate 
required  for  flash  vaporization  of  JP-8  for  mass  flows  at  the  same  operating  conditions. 
The  high  thermal  resistance  of  free  convection  and  radiation  heat  transfer  limit  the  heat 
transfer  possible.  As  will  be  shown,  the  heat  exchanger  tests  show  significantly  lower 
thermal  resistance  and  collect  much  higher  heat  loads  required  for  flash  vaporization. 
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Figure  36.  Plain  tube  heat  transfer  to  surroundings  vs.  axial  distance  w/  varying  equivalence  ratio 
Avgas  (298  K)  -air  (322  K),  15  Hz,  8  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 


Figure  37.  Plain  tube  heat  transfer  to  surroundings  vs.  axial  distance  w/  varying  equivalence  ratio 
JP-8  (298  K)  -  air  (395  K),  15  Hz,  8  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 

Water  FVS  with  Hydrogen-Air  Detonation 

Hydrogen  fueled  tests  with  water  FVS  were  required  to  prove  the  system  design 
and  detennine  initial  heat  transfer  rates.  For  the  first  test,  the  long  heat  exchanger  was 
used  and  spanned  from  0.876-1.638  m  measured  from  the  head.  Water  pressure  was 
maintained  at  6.205-6.895  MPa.  Water  mass  flow  was  varied  by  nozzle  size  until  a 
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steady-state  temperature  was  obtained  in  the  500  K-530  K  range.  At  a  water  mass  flow 
of  0.557  kg/min,  the  steady-state  water  outlet  temperature  of  5 12  K  was  achieved  as 
shown  in  Figure  38. 


Hydrogen  -  air  detonation,  phi  1.0, 10  Hz,  ignition  delay  6  ms 
1.829  m,  2"  SS  Sch  40  Tube  w/  0.305  m  spiral 
Water  mass  flow  0.557  kg/min,  heat  exchanger  location  (0.876-1.638  m) 

Steady-state  heat  transfer  rates  of  over  8  kW  calculated  as  shown  in  Figure  38. 

For  a  JP-8  fueled  test  at  the  same  operating  conditions  and  a  two-tube  configuration, 
using  equations  (1 7)-(2 1),  the  mass  flow  would  have  been  0.401  kg/min  and  the  heat 
transfer  required  to  obtain  the  same  outlet  temperatures  would  have  been  3.2 1 1  kW.  The 
high  heat  transfer  results  of  this  test  indicated  over  twice  the  required  heat  transfer. 

An  additional  hydrogen  fueled  test  was  completed  at  the  same  engine  operating 
conditions  with  lower  mass  flow  rates  and  lower  water  pressures.  This  test  was  originally 
intended  to  obtain  a  mass  flow  value  closer  to  the  JP-8  mass  flow  of  0.401  kg/min.  The 
reduction  in  water  pressure  resulted  in  boiling  within  the  heat  exchanger.  The  results  of 
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the  boiling  are  presented  so  that  the  unique  characteristics  are  known  and  can  be  avoided 
in  JP-8  FVS  testing. 

For  boiling  test  six  thennocouples  were  compression  clamped  to  the  heat 
exchanger,  three  at  the  top  and  three  at  the  bottom  of  the  heat  exchanger.  The  axial 
positions  were  determined  by  the  helical  coil  geometry  within  the  heat  exchanger  to 
ensure  that  the  thermocouples  were  placed  in  the  center  of  the  flow  channel.  See  Figure 
39  for  thermocouple  locations  and  terminology.  Measuring  from  the  head,  the  heat 
exchanger  spanned  from  0.152-0.914  m.  The  water  mass  flow  was  reduced  to  0.333 
kg/min.  The  water  pressure  recorded  at  the  outlet  was  4.413  MPa  during  the  test  which 
correlates  to  boiling  temperatures  of  530  K. 


Figure  39.  External  thermocouple  locations  for  hydrogen-air  detonation  with  boiling 
The  data  from  the  boiling  test  is  shown  in  Figure  40  with  temperatures  plotted 
versus  time.  Approximately  160  seconds  into  the  test  there  is  a  large  fluctuation  in  the 
thermocouple  readings.  Boiling  begins  in  this  time  frame  and  continues  for  the  duration 
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of  the  test.  The  fluctuations  in  the  various  temperatures  are  due  to  the  changing  flow 
characteristics  of  caused  by  boiling.  Looking  at  T-2  in  the  transient  portion,  the 
temperature  plateaus  momentarily  at  around  520  K  then  steadily  rises  to  temperatures  of 
561  K.  T-3  achieves  steady-state  temperatures  near  the  boiling  temperature  of  534  K. 
When  the  drop  through  the  wall  which  is  approximately  10-15  K  is  accounted  for  the 
boiling  may  extend  to  this  location.  T-l  and  all  three  lower  thennocouples  (not  plotted) 
remained  below  the  boiling  temperature  for  the  entire  test. 


Figure  40.  Temperature  vs.  time  for  boiling  test 
Hydrogen  -  air  detonation,  phi  1, 10  Hz,  ignition  delay  6  ms 
1.829  m,  2"  SS  Sch  40  Tube  w/  0.305  m  spiral 
Water  mass  flow  0.333  kg/min,  heat  exchanger  location  (0.152-0.914  m) 

The  temperature  readings  corresponding  to  the  time  of  the  highest  external 
thermocouple  reading  (632  seconds  into  the  test)  are  shown  versus  axial  location  in 
Figure  41.  The  boiling  temperature  is  also  plotted.  Note  that  the  external  temperatures 
do  not  account  for  the  10-15  K  loss  through  the  wall.  Boiling  occurs  locally  throughout 
the  top  portion  of  the  heat  exchanger  from  the  middle  to  the  outlet.  The  bottom  of  the 
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tube  experiences  a  linear  temperature  increase  throughout  the  length  of  the  heat 
exchanger. 


Figure  41.  Steady-state  external  surface  temperatures  632  seconds  into  the  boiling  test 
Hydrogen  -  air  detonation,  phi  1, 10  Hz,  ignition  delay  6  ms 
1.829  m,  2"  SS  Sch  40  Tube  w/  0.305  m  spiral 
Water  mass  flow  0.333  kg/min,  heat  exchanger  location  (0.152-0.914  m) 

From  the  top  and  bottom  temperature  profile  it  can  be  concluded  that  the  helical 
coil  which  was  meant  to  direct  flow  in  a  helical  path  was  not  functioning  as  designed.  If 
the  coil  had  performed  as  designed  the  temperature  profile  would  show  a  near  linear 
temperature  rise  from  inlet  to  the  outlet  on  both  the  top  and  bottom  as  the  flow  travel  in  a 
helical  path  to  the  outlet. 

This  temperature  data,  which  was  not  available  for  first  test,  provided  an 
additional  tool  for  monitoring  boiling  conditions.  For  the  remainder  of  the  tests  either 
radial  or  spanwise  temperature  distributions  were  monitored. 
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Water  FVS  with  Avgas- Air  Detonation 


The  first  hydrogen  fueled  water  FVS  test  proved  the  FVS  design  and  that  there 
was  more  than  enough  heat  transfer  to  the  fluid.  To  determine  the  heat  transfer  of  a 
liquid  hydrocarbon  fueled  detonation,  water  FVS  tests  were  completed  with  avgas  as  the 
detonation  fuel. 

One  steady-state  test  was  completed  with  the  long  heat  exchanger  using  avgas-air 
detonations.  The  equivalence  ratios  varied  from  1.04-1.10  throughout  the  duration  of  the 
test.  The  engine  frequency  was  15  Hz  with  a  4  ms  ignition  delay.  Water  mass  flow  was 
0.837  kg/min.  Water  pressure  was  set  to  6.619  MPa  corresponding  to  a  boiling 
temperature  of  556  K.  The  temperature  and  heat  transfer  results  are  shown  in  Figure  42. 


Figure  42.  Temperature  and  heat  transfer  vs.  time 
Avgas  (298  K)  -  air  (322  K)  detonation,  phi  1.04-1.10, 15  Hz,  ignition  delay  4  ms 
Water  mass  flow  0.837  kg/min,  heat  exchanger  location  (0.876-1.638  m) 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 


The  steady-state  temperature  of  540  K  was  achieved  with  calculated  heat  transfer 
rates  of  over  14  kW.  For  a  JP-8  fueled  test  at  the  same  operating  conditions  and  a  two- 
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tube  configuration,  using  equations  (17)-(21),  the  mass  flow  would  have  been  0.582 
kg/min  and  the  heat  transfer  required  to  obtain  the  same  outlet  temperatures  would  have 
been  5.693  kW.  Again  the  heat  generated  was  over  twice  the  required  amount  to  flash 
vaporize  the  JP-8  mass  flow  for  the  operating  conditions.  Based  on  the  results  of  this  test 
and  the  previous  water  cooled  tests  a  shorter  heat  exchanger  was  required  to  reduce  the 
amount  of  heat  transfer  to  the  working  fluid. 

Noting  the  high  thermal  stratification  in  the  previous  water  FVS  tests,  the  radial 
temperature  distribution  was  monitored  by  six  thermocouples  spaced  at  36°  increments 
from  0°  to  180°  as  shown  in  Figure  43  for  all  short  heat  exchanger  tests.  The  temperature 
profile  was  taken  12.7  cm  prior  to  the  heat  exchanger  outlet  unless  otherwise  noted.  An 
internal  thermocouple  at  the  same  axial  location  recorded  the  internal  flow  temperature  at 
the  top  of  the  heat  exchanger. 


Figure  43.  Thermocouple  locations  for  radial  temperature  profile 

For  the  first  avgas-air  test  with  the  short  heat  exchanger,  the  engine  frequency  was 
1 5  Hz  with  an  ignition  delay  of  6  ms.  Equivalence  ratio  varied  from  1 .05  to  1.11.  The 
heat  exchanger  spanned  from  1.130-1.511  m  measured  from  the  head. 
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With  a  water  mass  flow  of  0.523  kg/min  and  steady-state  temperature  of  480  K 
the  calculated  heat  transfer  rate  to  the  water  was  7  kW.  Using  equations  ( 1 7)-(2 1)  and 
the  same  engine  operating  parameters  and  two-tube  configuration  the  JP-8  mass  flow  and 
heat  transfer  would  be  0.572  kg/min  and  4. 14  kW. 

The  final  avgas-air  test  with  the  short  heat  exchanger  had  a  measured  water  mass 
flow  of  0.364  kg/min  as  compared  to  the  JP-8  flow  rate  of  0.572  kg/min  for  the  engine 
operating  parameters  and  two-tube  configuration,  per  equations  (17)-(21).  The  steady- 
state  water  outlet  temperature  was  540  K;  very  closely  matching  the  desired  JP-8  flash 
vaporization  temperature.  The  heat  transfer  rate  of  6.4  kW  was  generated  and  is  shown 
in  Figure  44.  For  the  JP-8  mass  flow  and  the  same  inlet  and  outlet  temperatures  the  JP-8 
requires  5.48  kW  of  power.  This  test  marked  the  first  time  that  the  heat  transfer  rate 
closely  matched  the  heat  transfer  rate  required  for  JP-8  with  the  same  inlet  and  outlet 
temperatures. 


"Water  Inlet 
"Water  Outlet 


Figure  44.  Inlet/Outlet  and  radial  temperature  vs.  time 
Avgas  (298  K)  -  air  (322  K)  detonation,  phi  1.06, 15  Hz,  ignition  delay  6  ms 
Water  mass  flow  0.364  kg/min,  heat  exchanger  location  (1.130-1.511  m) 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
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From  Figure  44  the  radial  temperature  differential  between  consecutive  radially 
spaced  thermocouples  appears  to  be  nearly  even.  The  radial  temperature  distribution  also 
follows  the  same  trends  seen  in  the  outlet  temperature.  The  radial  temperature  profile  at 
steady-state  conditions  is  shown  in  Figure  45.  The  steady-state  temperature  profile  had  a 
230  K  thennal  stratification  ranging  from  3 10  K  at  the  bottom  of  the  tube  to  the  540  K  at 
the  top  of  the  tube.  The  temperature  profile  showed  a  near  linear  growth  with  radial 
position. 

As  will  be  shown  later,  the  high  thennal  stratification  is  attributed  to  buoyancy 
forces.  For  the  final  outlet  temperature  to  so  closely  reflect  the  temperature  of  the  upper 
surface  there  must  be  an  asymmetric  velocity  profile  within  the  annular  region  with  the 
hottest  fluids  at  the  top  of  the  heat  exchanger  traveling  at  the  highest  velocities.  This  is 
an  unexpected  beneficial  feature  of  the  flow  which  reduces  the  residence  time  of  the 
hottest  fluids.  The  higher  velocities  in  the  hot  fluid  also  contribute  to  higher  local  heat 
transfer  coefficients  which  will  increase  heat  transfer  rates. 
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Figure  45.  Radial  Surface  Temperature  Profile  (K)  25.4  cm  Downstream  of  Inlet 
Avgas  (298  K)  -  air  (322  K)  detonation,  phi  1.06, 15  Hz,  ignition  delay  6  ms 
Water  mass  flow  0.364  kg/min,  heat  exchanger  location  (1.130-1.511  m) 

1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 

JP-8  Flash  Vaporization  System  Tests 

JP-8  FVS  tests  began  with  the  short  heat  exchanger  spanning  from  0.470-0.85 1  m. 
The  first  test  was  completed  with  equivalence  ratios  from  1 .0- 1 . 1  with  the  air  heated  to 
394  K.  The  engine  frequency  to  maintain  the  desired  equivalence  ratio  was  1 1.42  Hz  and 
fuel  mass  flow  of  0.355  kg/min.  To  prevent  boiling  of  the  JP-8  the  fuel  pressure  was 
always  a  values  above  the  critical  pressure. 

Steady-state  outlet  temperatures  of  625  K  were  achieved,  demonstrating  the  first 
successful  flash  vaporization  of  a  JP-8-air  mixture  using  PDE  waste  heat.  Upon 
achieving  steady-state  heat  exchanger  outlet  temperatures  the  air  heater  was  shut  off  to 
determine  if  heated  air  is  required  after  the  fuel  has  been  heated. 

Temperature  data  and  heat  transfer  rates  are  shown  in  Figure  46.  Steady-state 
heat  exchanger  outlet  temperatures  of  over  630  K  were  achieved  after  the  air  heater  was 
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shut  off  with  fuel  injection  temperatures  of  607  K  due  to  losses  from  the  heat  exchanger 
to  the  fuel  inlet  manifold.  The  rise  in  the  steady-state  temperature  is  attributed  to  a  lower 
equivalence  ratio  which  provided  better  engine  performance  as  will  be  shown. 

Steady-state  heat  transfer  rates  of  5  kW  were  observed.  Note  that  steady-state 
conditions  were  maintained  for  over  10  minutes  and  the  duration  of  the  test  was  only 
limited  by  the  fuel  storage  capacity. 


Figure  46.  Temperature  and  heat  transfer  vs.  time 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.470-0.851  m) 


It  can  be  seen  in  Figure  47,  that  fuel-air  mixture  temperatures  exceeded  the  400  K 


threshold  required  for  the  fuel  to  remain  in  the  vapor  state  after  mixing  (Tucker,  2005). 


In  this  test  the  mixture  temperature  did  not  achieve  steady-state  due  to  heating  of  the  inlet 


manifold.  Immediately  prior  to  the  air  heater  being  shut  off  the  mixture  temperature  was 


seven  degrees  Kelvin  hotter  than  the  upstream  air  temperature.  After  the  air  heater  was 


shut  off  the  transient  temperature  difference  between  mixture  temperature  and  upstream 
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air  temperature  increased  to  45  degree  Kelvin.  While  flash  vaporization  was  still 
occurring  locally  at  the  nozzles;  the  mixture  was  dual  phase  at  temperatures  below  the 
mixture  dew  point.  The  fuel  ran  out  prior  to  reaching  steady-state  mixture  temperature 
for  ambient  upstream  air  temperatures. 
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Figure  47.  Mixture  and  Upstream  Air  Temperature 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.470-0.851  m) 

A  twenty  percent  drop  in  mass  flow  occurred  over  the  duration  of  the  test  as 
shown  in  Figure  48.  For  the  fuel  nozzles  used,  the  mass  flow  is  proportional  to  the 
square  root  of  density  and  as  a  result  of  heating  the  fuel  to  near  supercritical  temperatures 
the  reduction  in  density  resulted  in  a  large  drop  in  the  mass  flow. 
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Figure  48.  Fuel  mass  flow  vs.  time 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w / 1.219  m  spiral 
Fleat  exchanger  location  (0.470-0.851  m) 


Looking  at  normalized  mass  flow  and  the  normalized  square  root  of  the  AFRL 
SUPERTRAPP  JP-8  surrogate  density  for  the  fuel  outlet  temperature,  shown  in  Figure 
49,  it  is  clear  the  correlation  between  the  two  properties.  A  drawback  of  the  fuel  system 
used  was  that  the  fuel  pressure  could  not  be  increased  during  the  test  to  account  for  the 
change  in  density  to  maintain  constant  mass  flow.  To  vary  equivalence  ratio  the 
frequency  was  adjusted,  effectively  increased  or  decreased  air  mass  flow  to  achieve  the 
desired  equivalence  ratios. 
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Figure  49.  Normalized  Mass  Flow  and  Density12 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.470-0.851  m) 

The  radial  temperature  profile  12.7  cm  from  the  outlet,  shown  Figure  50, 
reiterates  the  strong  thermal  gradients  present  in  the  heat  exchanger.  A  350  K  difference 


between  the  top  and  bottom  of  the  heat  exchanger  was  present  at  steady-state  conditions 


for  this  JP-8  FVS  test. 
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Figure  50.  Surface  Temperature  Profile  (K)  25.4  cm  Downstream  of  Inlet 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.470-0.851  m) 

Wave  speed  data  was  developed  for  this  test  as  a  function  of  time  and  is  presented 
in  Figure  5 1 .  From  the  initial  conditions  to  the  flash  vaporization  state  there  were 
negligible  changes  in  wave  speed.  When  the  air  heater  was  shut  off  there  was  a  drastic 
increase  in  wave  speed.  Also  plotted  with  time  is  the  equivalence  ratio  which  decreases 
due  to  the  decreased  mass  flow.  The  wave  speed  peaks  at  overdriven  velocities  at 
stoichiometric  mixtures  of  flash  vaporized  JP-8  and  then  decreases  but  is  still  near  the  CJ 
velocity  as  equivalence  ratios  become  lean.  Remember  from  the  plain  tube  tests  it  was 
not  possible  to  obtain  C  J  velocities  at  stoichiometric  and  lean  mixtures  with  unheated  JP- 
8. 
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Figure  51.  Wave  speed  and  equivalence  ratio  vs.  time 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.470-0.851  m) 

Supercritical  JP-8  Tests 

Having  achieved  flash  vaporization  temperatures,  the  scope  of  this  work  was 
expanded  to  see  what  JP-8  fluid  temperatures  could  be  achieved  with  the  FVS  and  what 
complications  or  benefits  arose  from  the  elevated  temperatures. 

To  obtain  the  highest  heat  transfer  to  the  JP-8  the  short  heat  exchanger  tests 
moved  to  the  hot  section  spanning  from  1 .130-1.5 11m  measured  from  the  head.  Engine 
frequency  started  at  above  16  Hz  but  was  adjusted  throughout  the  test  until  thermal 
equilibrium  was  reached  at  which  time  the  frequency  was  held  constant  at  13.88  Hz  at  an 
equivalence  ratio  of  unity.  The  air  was  heated  to  394K  and  ignition  delay  was  constant  at 
6  ms. 

Steady-state  heat  exchanger  outlet  temperatures  of  over  700  K,  with  fuel  injection 
temperatures  of  695  K  and  heat  transfer  rates  of  8.85  kW  were  achieved  as  shown  in 
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Figure  52.  Note  that  the  steady-state  fuel  injection  temperature  was  well  above  the  JP-8 
critical  temperature  of  680  K.  This  was  the  first  recorded  use  of  supercritical  JP-8  in  a 
PDE  engine. 


- Fuel  Manifold  Inlet 

- Fuel  HX  Inlet 

- 12.7  cm  Flow  TC 

- 25.4  cm  Flow  TC 

- Fuel  HX  Outlet 

111  “  Super  Critical 
Temperature 
—  Heat  T  ransfer 


Figure  52.  Fuel  flow  temperature  vs.  time 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w / 1.219  m  spiral 
Heat  exchanger  location  (1.130-1.511  m) 


The  fuel-air  mixture  temperature  reached  a  steady-state  temperature  of  428  K.  A 


radial  temperature  profile  was  obtained  at  the  outlet  with  temperature  stratification 


varying  from  367  K  at  the  bottom  of  the  tube  to  678  K  at  the  top  of  the  tube. 
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Figure  53.  Radial  temperature  (K)  profile 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w / 1.219  m  spiral 
Heat  exchanger  location  (1.130-1.511  m) 


Mixture  temperature  and  wave  speed  are  plotted  versus  time  in  Figure  54.  It  can 
be  seen  that  as  the  mixture  temperature  approaches  the  fully  vaporized  mixture 
temperature  of  400  K  that  thrust  tube  wave  speed  at  the  end  of  the  spiral  increase  to 
detonation  wave  speeds  of  greater  than  1700  m/s  (Glassman,  1996:  247). 
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Figure  54.  Wave  speed  and  mixture  temperature  vs.  time 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (1.130-1.511  m) 


The  drop  in  superheated  fuel  density  resulted  in  one  third  drop  in  fuel  mass  flow 
across  the  nozzles  due  to  the  mass  flow  proportionality  to  the  square  root  of  fuel  density. 
Figure  55  shows  normalized  values  of  heat  exchanger  outlet  temperature,  the  square  root 
of  AFRL  SUPERTRAPP  JP-8  surrogate  density  for  the  outlet  temperature,  and  mass 
flow.  Again  the  square  root  of  the  density  tracks  very  well  mass  flow  as  a  function  of 
time. 


83 


Figure  55.  Normalized  fuel  manifold  inlet  temperature, 
mass  flow,  and  square  root  of  density 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (1.130-1.511  m) 

The  final  JP-8  FVS  test  was  completed  with  the  long  heat  exchanger  installed  at 
the  hot  section  spanning  from  0.797-1.559  m  measure  from  the  head.  Engine  frequency 
began  at  20  Hz  but  was  adjusted  throughout  the  test  until  thermal  equilibrium  was 
reached  at  which  time  the  frequency  was  held  constant  at  12  Hz.  The  air  was  heated  to 
394K  and  ignition  delay  was  constant  at  6  ms.  Equivalence  ratios  varied  through  the  test 
with  frequency  and  temperature  changes  with  steady-state  equivalence  ratios  of 
approximately  1.0. 

Supercritical  steady-state  temperatures  of  760  K  were  achieved  with  a  calculated 
heat  transfer  rate  of  7.9  kW.  The  fuel-air  mixture  reached  temperatures  of  448  K  which 
is  near  the  auto-ignition  temperature  of  483  K. 

The  supercritical  mass  flow  dropped  significantly  throughout  the  duration  of  the 
test  with  final  steady-state  mass  flows  approximately  40%  of  the  starting  value.  The 
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square  root  of  the  AFRL  SUPERTRAPP  JP-8  surrogate  density  interpolated  for  the  fuel 
manifold  inlet  temperate  exhibits  the  same  trends  as  the  mass  flow,  but  a  divergence 
between  the  two  data  sets  is  apparent  in  the  data  sets  shown  in  Figure  56. 


Figure  56.  Normalized  fuel  manifold  inlet  temperature, 
mass  flow,  and  square  root  of  density 
JP-8  -  air  (394  K)  detonation,  6  ms  ignition  delay 
1.829  m,  2"  SS  Sch  40  Tube  w/ 1.219  m  spiral 
Heat  exchanger  location  (0.797-1.559  m) 


Experimental  Supercritical  Density  Calculation 

For  the  flow  nozzles  used,  the  mass  flow  is  proportional  to  the  square  root  of  the 
pressure  differential  and  the  square  root  of  the  fluid  density.  The  proportionality  is 
multiplied  by  a  constant  to  obtain  the  mass  flow.  This  constant  is  termed  the  flow 
number  (FN).  The  FN  varies  with  the  nozzle  flow  area  and  fluid  but  may  be  determined 
experimentally  with  known  mass  flows,  pressures,  and  densities.  The  three  forms  of  the 
nozzle  mass  flow  equation  (Bartok  and  Sarofim,  1991:553)  are 
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(45) 


“dot 

FN  = - 

4^p-4p 

“dot  =  FN-^P.^ 


“dot 

AP-FN2 


(46) 


(47) 


For  low  temperature  JP-8  testing  the  mass  flow,  pressure  differential,  and  density 
are  known  allowing  for  the  calculation  of  FN  for  a  set  of  fuel  nozzles.  This  is  termed 
calibration  FN  which  is  assumed  to  remain  constant.  For  all  JP-8  tests,  the  fuel  pressure 
was  constant  and  the  mass  flow  was  obtained  upstream  of  the  heat  exchanger  allowing 
for  the  density  to  be  calculated. 

For  the  above  method  to  be  valid  during  the  transient  portion  of  the  test  it  must  be 
proven  that  the  change  in  mass  flow  due  to  the  density  change  may  be  neglected.  This 
analysis  was  not  completed  under  this  work.  Therefore,  this  method  is  presented  as  an 
experimental  supercritical  density  calculation  only  at  steady-state  conditions  with  the 
assumption  that  the  above  equations  are  remain  valid  in  the  supercritical  temperature 
regime  tested. 

Free  Convection  versus  Forced  Convection 

The  initial  heat  exchanger  design  analysis  assumed  heat  transfer  due  to  uniform 
flow.  After  the  first  water-cooled  tests  it  was  apparent  from  the  thermal  gradients  that 
uniform  flow  could  not  be  assumed.  An  analysis  was  completed  to  detennine  whether 
the  buoyant  or  inertial  forces  dominated  flow  characteristics.  A  temperature  range  from 
290  K  and  650  K  was  used  for  the  analysis.  The  mass  flow  rate  varied  from  0.3  kg/min 
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to  0.6  kg/min.  The  property  data  used  for  the  analysis  was  the  SUPERTRAPP  JP-8 
surrogate  data  presented  in  Appendix  A. 

For  the  analysis  the  Reynolds  (Re)  number  and  Grashof  (Gr)  number  were 
required.  The  Reynolds  number  is  the  ratio  of  inertia  forces  to  the  viscous  forces,  while 
the  Grashof  number  is  the  ratio  of  buoyant  to  inertia  forces. 

It  can  be  shown  that  the  Reynolds  number  for  the  heat  exchanger  is 

4mdot 

Re=  - r- 

71  woi  +  dio)'P  (4; 

Where 


dio  60.3  5 111111  [nner  tube  outer  diameter 


doi  66.9_0nm  qu1c[.  jnner  diameter 


F  =  Dynamic  viscosity  (N-s/m  ) 


The  Grashof  number  is 


&'P'(Toi  Tfuel)'(doi  dio) 


Where 

T  • 

01  =  Inner  tube  outer  diameter  surface  temperature  (  K) 

Tfuel  =  pue[  temperature  at  Vi  thickness  of  annular  region  (  K) 

[>--(£) 

P  V  ai  J  Expansion  coefficient  ( 1/  K) 


(49) 
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Acceleration  due  to  gravity 


m 

g  =  9.807— 

2 
s 

2 

v  =  Kinematic  viscosity  (m  /2) 

The  temperature  data  for  the  inner  tube  outer  diameter  surface  temperature  was  estimated 
to  be  5  K  hotter  than  the  fuel  temperature. 

The  ratio  of  Grashof  number  to  Reynolds  number  squared  provides  insight  as  to 
whether  free  convection  due  to  buoyancy  or  forced  convection  dominate.  For  ratios  of 
much  less  than  1  the  flow  is  dominated  by  forced  convection.  For  ratios  much  greater 
than  1  the  flow  is  dominated  by  free  convection.  For  ratios  of  approximately  one  the 
combined  effects  of  free  and  force  convection  must  be  considered  (Incropera,  1996:451). 
The  ratio  of  Grashof  number  to  Reynolds  number  squared  is  presented  for  the 
temperature  range  in  Figure  57. 


Figure  57.  Ratio  of  Grashof  number  to  Reynolds  number  squared  for  simulated  mass  flows  and 

temperature 
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For  the  simulated  mass  flows,  which  span  the  range  for  all  tests  completed,  the 
ratio  indicates  that  buoyant  forces  must  be  considered.  This  analysis  explains  the  high 
thermal  gradients  witnessed. 

Coking/Deposits 

Clean  fuel  nozzles  were  photographed  prior  to  installation  of  J-8  FVS  testing. 
Two  supercritical  JP-8  tests  for  a  total  test  time  of  less  than  thirty  minutes  were 
completed.  The  nozzles  were  then  removed  and  photographed.  The  photographs  are 
presented  in  Figures  58  and  59.  It  is  apparent  from  the  blackening  of  the  nozzles  in  the 
post  test  condition  that  some  deposits  have  been  collected  on  the  fuel  nozzles. 
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Figure  58.  JP-8  fuel  nozzles  pretest  condition 
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Figure  59.  JP-8  fuel  nozzles  posttest  condition 
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The  deposits  visible  on  the  nozzles  prompted  interest  into  what  was  occurring  in 
the  heat  exchanger.  The  long  heat  exchanger  was  no  longer  required  after  the  completion 
of  testing  and  was  disassembled  to  visually  determine  if  fuel  deposits  were  a  significant 
issue.  The  disassembled  heat  exchanger  is  pictured  in  Figure  60.  It  can  be  seen  that 
surface  deposits  occurred  throughout  the  heat  exchanger  with  the  exception  of  the  local 
inlet  area  which  was  being  cooled  by  the  inlet  flow.  It  can  also  be  seen  that  the  helical 
coil  did  not  restrict  the  flow  to  the  helical  path  as  was  hypothesized  due  to  the  high 
thermal  stratification. 


Heat  Exchanger  Outlet 


Heat  Exchanger  Inlet 

Figure  60.  Disassembled  long  heat  exchanger  with  carbon  deposits 
Fuel  coking  and  deposits  did  not  appear  to  affect  any  of  the  experimental  data,  but 
it  could  be  seen  in  the  fuel  nozzles  and  long  heat  exchanger  that  deposits  did  occur  during 
JP-8  FVS  testing. 


90 


V.  Conclusions  and  Recommendations 


This  work  demonstrated  the  first  successful  use  of  the  waste  heat  generated  by  a 
PDE  engine  for  the  self-sustained  operation  and  flash  vaporization  of  a  JP-8-air  mixture. 
The  duration  of  the  steady-state  tests  exceeded  the  operating  time  of  any  previous  JP-8 
fueled  PDE’s  with  run  times  in  excess  of  twenty  minutes  and  limited  only  by  fuel  storage 
capacity. 

The  flash  vaporization  system  used  permits  research  to  be  completed  using  liquid 
hydrocarbon  fuels-air  mixtures  at  steady-state  conditions.  The  practicality  of  the  flash 
vaporization  system  used  also  serves  as  a  proof-of-concept  for  future  PDE  aerial  vehicles 
as  an  efficient,  simple  method  for  the  use  of  liquid  hydrocarbon  fuels. 

This  work  also  demonstrated  the  first  successful  use  of  supercritical  JP-8  in  a 
PDE  engine.  The  ability  to  heat  the  fuel  to  supercritical  temperatures  without  prohibitive 
amounts  of  fuel  deposits  or  fuel  thermal  stability  issues  are  attributed  to  the  steady-state 
flow  and  the  low  residence  time  of  the  fuel  in  the  heat  exchanger  and  feed  lines.  Heating 
the  fuel  to  supercritical  temperatures  and  detonating  the  superheated  fuel  provides  an 
important  step  to  using  endothermic  fuels  and/or  experiencing  endothermic  reactions. 

For  JP-8  FVS  tests  fuel  injection  temperatures  of  630  -760  K  were  achieved  with 
heat  transfer  rates  from  the  heat  exchanger  ranging  from  5. 0-7. 9  kW.  The  time  to  achieve 
steady-state  conditions  for  either  heat  exchanger  was  less  than  five  minutes. 

The  heat  exchanger  design  did  not  account  for  buoyancy  effects  which  were 
shown  to  be  a  contributing  factor  in  the  test  configuration.  The  asymmetric  temperature 
and  velocity  profile  of  the  fuel  in  the  heat  exchanger,  due  to  the  contributions  on  forced 
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and  free  convection,  led  to  significantly  higher  heat  transfer  rates  allowing  for 
supercritical  temperatures  to  be  reached. 

Another  milestone  of  this  work  was  the  development  of  a  method  to 
experimentally  detennine  the  density  of  supercritical  fuels  by  the  use  of  mass  flow  rate 
data,  temperature  data,  and  pressure  data.  The  method  presented  is  only  valid  in  during 
steady-state  operation  and  assumes  the  governing  equations  are  valid  in  the  supercritical 
temperature  range. 

Using  the  AFRL  SUPERTRAPP  JP-8  surrogate,  a  complete  set  of  thennodynamic 
property  data  including  density,  viscosity,  specific  heat,  enthalpy,  and  thermal 
conductivity  were  generated  for  temperatures  ranging  from  273  -  998  K  and  pressures 
from  1-85  atmospheres.  The  tables  provide  a  valuable  tool  in  analysis  and  compare 
well  with  existing  thermodynamic  properties  for  JP-8  at  low  temperatures  and  with 
existing  thermodynamic  property  prediction  models  at  elevated  temperatures. 

To  the  author’s  knowledge  the  AFRL  SUPERTRAPP  JP-8  surrogate  tables 
developed  and  presented  in  Appendix  A  encompass  the  largest  surrogate  data  set 
available  for  JP-8.  In  lieu  of  actual  JP-8  thermodynamic  data  this  data  provides 
reasonable  approximation  for  analytical  research  and  experimental  predictions. 

Recommendations 

The  fuel  pressurization  system  used  for  this  work  was  limited  to  a  constant 
pressure  during  testing.  For  all  of  the  JP-8  tests  the  reduction  in  the  density  due  to  the 
temperature  dependence  prevented  the  mass  flow  from  remaining  constant  until  steady- 
state  conditions  were  reached.  With  this  restriction  the  PDE  frequency  was  varied  to 
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obtain  the  desired  equivalence  ratio.  The  change  in  engine  frequency  and  equivalence 
ratios  throughout  the  transient  portion  of  the  test  was  prohibitive  in  detennining  the 
transient  PDE  characteristics.  It  is  recommended  that  a  variable  fuel  pressure  system  be 
incorporated  into  the  existing  PDE  configuration  to  allow  for  constant  mass  flow  and 
equivalence  ratios  in  further  high  temperature  liquid  hydrocarbon  testing. 

Supercritcal  temperatures  were  achieved  in  this  work  with  the  use  of  only  one 
tube  mounted  heat  exchanger  feeding  mass  flow  for  a  two-tube  PDE  configuration. 
Significantly  higher  temperatures  are  possible  by  the  use  of  one  heat  exchanger  per  tube 
running  the  mass  flow  only  for  that  tube.  It  may  be  possible  with  such  a  configuration  to 
achieve  the  temperatures  required  for  endothermic  reactions.  Just  as  performance 
benefits  were  shown  by  Tucker  (Tucker,  2005)  to  occur  when  using  a  flash  vaporized 
mixture  the  same  may  be  shown  by  achieving  endothermic  reactions  in  the  fuel. 

The  heat  transfer  to  the  fuel  was  significantly  higher  than  the  finite  difference 
method  model  predicted.  The  model  was  based  on  a  uniform  flow  force  convection  flow 
calculation  when  in  actuality  the  flow  had  over  a  300  K  circumferential  temperature 
gradient  indicating  a  strong  asymmetric  flow  profile  within  the  heat  exchanger.  The 
strong  thermal  gradient  and  Grashof  number  calculations  also  indicated  that  internal  heat 
transfer  was  dominated  by  free  convection  rather  than  forced  convection.  It  is 
recommended  that  the  fluid  velocity  profile  and  free  convective  heat  transfer  be  modeled 
to  determine  if  a  more  efficient  heat  exchanger  design  is  possible.  With  a  more  efficient 
heat  exchanger,  higher  temperatures  and  smaller  lighter  designs  are  possible. 
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It  is  recommended  that  further  research  be  completed  to  validate  the  method  and 
assumptions  used  to  experimentally  detennine  the  liquid  hydrocarbon  density  at 
supercritical  temperatures.  If  the  method  is  validated,  it  is  recommended  that  a  complete 
set  of  supercritical  density  tables  with  functions  of  temperature  and  pressure  be 
experimentally  generated.  This  method  is  not  limited  to  JP-8  and  may  be  incorporated  in 
determining  the  supercritical  density  properties  of  any  number  of  liquids  that  abide  by  the 
flow  nozzle  characteristics  presented. 
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Appendix  A.  AFRL  SUPERTRAPP  JP-8  Surrogate  Thermodynamic  Data 
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AFRL  SUPERTRAPP  JP-8  Surrogate 
Specific  Heat  [kJ/(kg-K)] 
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AFRL  SUPERTRAPP  JP-8  Surrogate 
Enthalpy  [kJ/kg] 
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1329 

1327 

1326 

1325 

1324 

1322 

1321 

1320 

1319 

1318 

1185 

1389 

1387 

1385 

1383 

1381 

1379 

1377 

1376 

1374 

1373 

1371 

1370 

1368 

1367 

1366 

1365 

1363 

1224 

1439 

1436 

1434 

1432 

1429 

1427 

1425 

1423 

1421 

1420 

1418 

1416 

1415 

1413 

1412 

1411 

1409 

1264 

1521 

1487 

1484 

1481 

1479 

1476 

1474 

1472 

1469 

1467 

1466 

1464 

1462 

1460 

1459 

1457 

1456 

1305 

1623 

1539 

1535 

1532 

1529 

1526 

1523 

1520 

1518 

1516 

1513 

1511 

1509 

1507 

1506 

1504 

1502 

1346 

1737 

1592 

1588 

1583 

1580 

1576 

1573 

1570 

1567 

1564 

1562 

1559 

1557 

1555 

1553 

1551 

1549 

1387 

1864 

1646 

1641 

1636 

1632 

1627 

1623 

1620 

1616 

1613 

1610 

1607 

1605 

1602 

1600 

1598 

1596 

1429 

1476 

1704 

1697 

1691 

1685 

1680 

1675 

1671 

1667 

1663 

1660 

1657 

1654 

1651 

1648 

1646 

1643 

1472 

1514 

1800 

1756 

1748 

1741 

1735 

1729 

1724 

1719 

1715 

1710 

1707 

1703 

1700 

1697 

1694 

1691 

1514 

1553 

1915 

1818 

1808 

1799 

1791 

1784 

1778 

1772 

1767 

1762 

1758 

1754 

1750 

1747 

1743 

1740 

1557 

1593 

2042 

1887 

1873 

1861 

1851 

1842 

1835 

1828 

1821 

1816 

1811 

1806 

1801 

1797 

1794 

1790 

1601 

1633 

1721 

1979 

1944 

1928 

1915 

1904 

1894 

1885 

1878 

1871 

1865 

1859 

1854 

1849 

1845 

1841 

1644 

1674 

1749 

2100 

2025 

2002 

1984 

1969 

1956 

1945 

1936 

1928 

1920 

1914 

1908 

1902 

1897 

1893 

1688 

1716 

1780 

2234 

2120 

2085 

2059 

2038 

2022 

2008 

1996 

1986 

1977 

1970 

1962 

1956 

1950 

1945 

1733 

1758 

1815 

1969 

2243 

2185 

2144 

2115 

2092 

2074 

2059 

2047 

2036 

2026 

2018 

2011 

2004 

1998 

1777 

1801 

1851 

1968 

2381 

2318 

2246 

2201 

2169 

2144 

2125 

2109 

2096 

2085 

2075 

2066 

2059 

2052 

1822 

1845 

1889 

1983 

2355 

2608 

2382 

2303 

2255 

2221 

2195 

2175 

2159 

2145 

2133 

2123 

2114 

2106 

1867 

1888 

1929 

2006 

2224 

3369 

2602 

2436 

2355 

2305 

2270 

2244 

2223 

2206 

2192 

2180 

2170 

2161 

1913 

1932 

1969 

2034 

2185 

3061 

3169 

2632 

2479 

2400 

2351 

2316 

2290 

2270 

2253 

2239 

2227 

2217 

1958 

1977 

2010 

2067 

2180 

2533 

5098 

2985 

2642 

2512 

2440 

2393 

2360 

2334 

2314 

2298 

2284 

2272 

2004 

2021 

2052 

2102 

2192 

2403 

3295 

3519 

2864 

2644 

2538 

2475 

2432 

2401 

2377 

2358 

2342 

2328 

2050 

2067 

2095 

2139 

2213 

2359 

2747 

3479 

3094 

2791 

2644 

2561 

2507 

2469 

2440 

2418 

2400 

2385 

2096 

2112 

2139 

2178 

2240 

2350 

2578 

3058 

3169 

2921 

2749 

2647 

2582 

2537 

2504 

2478 

2458 

2441 

2143 

2158 

2182 

2218 

2271 

2358 

2513 

2801 

3069 

2987 

2838 

2729 

2655 

2604 

2566 

2538 

2515 

2496 

2189 

2203 

2227 

2259 

2305 

2376 

2491 

2680 

2917 

2980 

2893 

2796 

2722 

2668 

2627 

2596 

2571 

2551 

2236 

2250 

2271 

2301 

2342 

2401 

2491 

2626 

2805 

2926 

2913 

2845 

2779 

2726 

2684 

2652 

2626 

2605 

2283 

2296 

2316 

2343 

2380 

2431 

2503 

2606 

2740 

2864 

2904 

2874 

2824 

2777 

2737 

2705 

2679 

2657 

2330 

2342 

2362 

2387 

2420 

2464 

2524 

2605 

2709 

2817 

2883 

2887 

2858 

2820 

2785 

2754 

2729 

2707 

2378 

2389 

2407 

2431 

2461 

2500 

2551 

2617 

2700 

2790 

2861 

2889 

2880 

2856 

2827 

2800 

2776 

2755 

2425 

2436 

2453 

2475 

2503 

2537 

2582 

2637 

2705 

2779 

2847 

2887 

2896 

2884 

2864 

2841 

2820 

2801 

2472 

2483 

2500 

2520 

2545 

2576 

2615 

2663 

2719 

2782 

2842 

2887 

2908 

2908 

2896 

2879 

2861 

2845 

2520 

2531 

2546 

2565 

2589 

2617 

2651 

2693 

2741 

2794 

2847 

2891 

2919 

2928 

2924 

2914 

2900 

2886 

2568 

2578 

2593 

2611 

2632 

2658 

2689 

2726 

2767 

2813 

2860 

2901 

2931 

2947 

2950 

2946 

2936 

2925 

2616 

2626 

2640 

2657 

2677 

2701 

2729 

2761 

2797 

2837 

2878 

2917 

2947 

2967 

2975 

2976 

2971 

2963 

2664 

2673 

2687 

2703 

2722 

2744 

2769 

2798 

2830 

2866 

2902 

2937 

2966 

2988 

3000 

3005 

3004 

2999 

2712 

2721 

2734 

2749 

2767 

2788 

2811 

2837 

2866 

2897 

2930 

2961 

2989 

3011 

3026 

3034 

3036 

3034 

2761 

2769 

2782 

2796 

2813 

2832 

2853 

2877 

2903 

2931 

2960 

2989 

3015 

3037 

3053 

3063 

3068 

3069 

2809 

2817 

2829 

2843 

2859 

2877 

2897 

2919 

2942 

2968 

2994 

3020 

3044 

3065 

3082 

3094 

3101 

3104 

2858 

2866 

2877 

2890 

2905 

2922 

2941 

2961 

2982 

3005 

3029 

3053 

3076 

3096 

3112 

3125 

3134 

3138 

2906 

2914 

2925 

2938 

2952 

2968 

2985 

3004 

3024 

3045 

3066 

3088 

3109 

3128 

3145 

3158 

3167 

3173 

2955 

2963 

2973 

2985 

2999 

3014 

3030 

3047 

3066 

3085 

3105 

3125 

3145 

3163 

3178 

3192 

3202 

3209 

3004 

3011 

3022 

3033 

3046 

3060 

3075 

3092 

3109 

3127 

3145 

3164 

3182 

3199 

3214 

3227 

3237 

3245 

3053 

3060 

3070 

3081 

3094 

3107 

3121 

3137 

3153 

3169 

3186 

3203 

3220 

3236 

3251 

3263 

3274 

3282 

3102 

3109 

3119 

3130 

3141 

3154 

3168 

3182 

3197 

3213 

3228 

3244 

3260 

3275 

3289 

3301 

3312 

3320 

3151 

3158 

3168 

3178 

3189 

3201 

3214 

3228 

3242 

3257 

3271 

3286 

3301 

3315 

3328 

3340 

3350 

3359 

3201 

3207 

3217 

3226 

3237 

3249 

3261 

3274 

3287 

3301 

3315 

3329 

3343 

3356 

3368 

3380 

3390 

3399 

3250 

3257 

3266 

3275 

3286 

3297 

3309 

3321 

3333 

3346 

3359 

3372 

3385 

3398 

3410 

3420 

3430 

3439 

3300 

3306 

3315 

3324 

3334 

3345 

3356 

3368 

3380 

3392 

3404 

3417 

3429 

3440 

3452 

3462 

3472 

3480 

3350 

3356 

3364 

3373 

3383 

3393 

3404 

3415 

3426 

3438 

3450 

3461 

3473 

3484 

3495 

3505 

3514 

3522 

3400 

3406 

3414 

3423 

3432 

3442 

3452 

3463 

3474 

3485 

3496 

3507 

3518 

3528 

3538 

3548 

3557 

3565 
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1 
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15 

20 

25 
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35 

Pressure  (atm) 
40  45 

50 

55 

60 

65 

70 

75 

80 

85 

0.137 

0.138 

0.138 

0.138 

0.138 

0.138 

0.139 

0.139 

0.139 

0.139 

0.139 

0.139 

0.140 

0.140 

0.140 

0.140 

0.140 

0.141 

0.137 

0.137 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.138 

0.138 

0.139 

0.139 

0.139 

0.139 

0.139 

0.140 

0.140 

0.140 

0.136 

0.136 

0.136 

0.136 

0.136 

0.137 

0.137 

0.137 

0.137 

0.137 

0.138 

0.138 

0.138 

0.138 

0.138 

0.139 

0.139 

0.139 

0.134 

0.135 

0.135 

0.135 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.137 

0.137 

0.137 

0.137 

0.137 

0.138 

0.138 

0.133 

0.133 

0.134 

0.134 

0.134 

0.134 

0.134 

0.135 

0.135 

0.135 

0.135 

0.135 

0.136 

0.136 

0.136 

0.136 

0.136 

0.137 

0.132 

0.132 

0.132 

0.132 

0.133 

0.133 

0.133 

0.133 

0.134 

0.134 

0.134 

0.134 

0.134 

0.135 

0.135 

0.135 

0.135 

0.135 

0.130 

0.131 

0.131 

0.131 

0.131 

0.132 

0.132 

0.132 

0.132 

0.132 

0.133 

0.133 

0.133 

0.133 

0.134 

0.134 

0.134 

0.134 

0.129 

0.129 

0.129 

0.130 

0.130 

0.130 

0.130 

0.131 

0.131 

0.131 

0.131 

0.132 

0.132 

0.132 

0.132 

0.132 

0.133 

0.133 

0.127 

0.128 

0.128 

0.128 

0.128 

0.129 

0.129 

0.129 

0.129 

0.130 

0.130 

0.130 

0.130 

0.131 

0.131 

0.131 

0.131 

0.132 

0.126 

0.126 

0.126 

0.127 

0.127 

0.127 

0.127 

0.128 

0.128 

0.128 

0.128 

0.129 

0.129 

0.129 

0.129 

0.130 

0.130 

0.130 

0.124 

0.125 

0.125 

0.125 

0.125 

0.126 

0.126 

0.126 

0.126 

0.127 

0.127 

0.127 

0.127 

0.128 

0.128 

0.128 

0.129 

0.129 

0.123 

0.123 

0.123 

0.123 

0.124 

0.124 

0.124 

0.125 

0.125 

0.125 

0.125 

0.126 

0.126 

0.126 

0.126 

0.127 

0.127 

0.127 

0.121 

0.121 

0.122 

0.122 

0.122 

0.122 

0.123 

0.123 

0.123 

0.124 

0.124 

0.124 

0.124 

0.125 

0.125 

0.125 

0.125 

0.126 

0.119 

0.119 

0.120 

0.120 

0.120 

0.121 

0.121 

0.121 

0.122 

0.122 

0.122 

0.122 

0.123 

0.123 

0.123 

0.124 

0.124 

0.124 

0.118 

0.118 

0.118 

0.118 

0.119 

0.119 

0.119 

0.120 

0.120 

0.120 

0.121 

0.121 

0.121 

0.122 

0.122 

0.122 

0.122 

0.123 

0.116 

0.116 

0.116 

0.117 

0.117 

0.117 

0.118 

0.118 

0.118 

0.119 

0.119 

0.119 

0.120 

0.120 

0.120 

0.120 

0.121 

0.121 

0.114 

0.114 

0.114 

0.115 

0.115 

0.115 

0.116 

0.116 

0.116 

0.117 

0.117 

0.117 

0.118 

0.118 

0.118 

0.119 

0.119 

0.119 

0.112 

0.112 

0.113 

0.113 

0.113 

0.114 

0.114 

0.114 

0.115 

0.115 

0.115 

0.116 

0.116 

0.116 

0.117 

0.117 

0.117 

0.118 

0.112 

0.110 

0.110 

0.111 

0.111 

0.112 

0.112 

0.112 

0.113 

0.113 

0.113 

0.114 

0.114 

0.114 

0.115 

0.115 

0.115 

0.116 

0.111 

0.108 

0.108 

0.108 

0.109 

0.109 

0.110 

0.110 

0.110 

0.111 

0.111 

0.111 

0.112 

0.112 

0.112 

0.113 

0.113 

0.114 

0.111 

0.105 

0.106 

0.106 

0.106 

0.107 

0.107 

0.108 

0.108 

0.108 

0.109 

0.109 

0.110 

0.110 

0.110 

0.111 

0.111 

0.111 

0.110 

0.103 

0.103 

0.103 

0.104 

0.104 

0.105 

0.105 

0.105 

0.106 

0.106 

0.107 

0.107 

0.107 

0.108 

0.108 

0.109 

0.109 

0.108 

0.100 

0.100 

0.101 

0.101 

0.102 

0.102 

0.103 

0.103 

0.103 

0.104 

0.104 

0.105 

0.105 

0.105 

0.106 

0.106 

0.107 

0.025 

0.097 

0.098 

0.098 

0.099 

0.099 

0.100 

0.100 

0.101 

0.101 

0.101 

0.102 

0.102 

0.103 

0.103 

0.103 

0.104 

0.104 

0.026 

0.095 

0.095 

0.096 

0.096 

0.097 

0.097 

0.098 

0.098 

0.099 

0.099 

0.100 

0.100 

0.100 

0.101 

0.101 

0.102 

0.102 

0.027 

0.093 

0.093 

0.094 

0.094 

0.095 

0.095 

0.095 

0.096 

0.096 

0.097 

0.097 

0.098 

0.098 

0.099 

0.099 

0.099 

0.100 

0.028 

0.091 

0.091 

0.091 

0.092 

0.092 

0.093 

0.093 

0.094 

0.094 

0.095 

0.095 

0.095 

0.096 

0.096 

0.097 

0.097 

0.098 

0.029 

0.091 

0.088 

0.089 

0.089 

0.090 

0.090 

0.091 

0.091 

0.092 

0.092 

0.093 

0.093 

0.094 

0.094 

0.095 

0.095 

0.096 

0.030 

0.090 

0.086 

0.087 

0.087 

0.088 

0.088 

0.089 

0.089 

0.090 

0.090 

0.091 

0.091 

0.092 

0.092 

0.093 

0.093 

0.094 

0.031 

0.089 

0.084 

0.084 

0.085 

0.086 

0.086 

0.087 

0.087 

0.088 

0.088 

0.089 

0.089 

0.090 

0.090 

0.091 

0.091 

0.092 

0.032 

0.033 

0.081 

0.082 

0.083 

0.083 

0.084 

0.085 

0.085 

0.086 

0.086 

0.087 

0.087 

0.088 

0.088 

0.089 

0.089 

0.090 

0.033 

0.034 

0.081 

0.080 

0.081 

0.081 

0.082 

0.083 

0.083 

0.084 

0.084 

0.085 

0.085 

0.086 

0.086 

0.087 

0.087 

0.088 

0.034 

0.035 

0.080 

0.078 

0.078 

0.079 

0.080 

0.080 

0.081 

0.082 

0.082 

0.083 

0.084 

0.084 

0.085 

0.085 

0.086 

0.086 

0.035 

0.036 

0.079 

0.075 

0.076 

0.077 

0.078 

0.078 

0.079 

0.080 

0.080 

0.081 

0.082 

0.082 

0.083 

0.083 

0.084 

0.084 

0.036 

0.037 

0.038 

0.074 

0.074 

0.075 

0.076 

0.077 

0.077 

0.078 

0.079 

0.079 

0.080 

0.080 

0.081 

0.082 

0.082 

0.083 

0.037 

0.038 

0.039 

0.073 

0.072 

0.073 

0.074 

0.075 

0.075 

0.076 

0.077 

0.077 

0.078 

0.079 

0.079 

0.080 

0.081 

0.081 

0.038 

0.039 

0.040 

0.073 

0.070 

0.071 

0.072 

0.073 

0.073 

0.074 

0.075 

0.076 

0.076 

0.077 

0.078 

0.078 

0.079 

0.079 

0.039 

0.040 

0.041 

0.042 

0.069 

0.068 

0.070 

0.071 

0.072 

0.072 

0.073 

0.074 

0.075 

0.075 

0.076 

0.077 

0.077 

0.078 

0.040 

0.041 

0.042 

0.043 

0.068 

0.066 

0.067 

0.068 

0.069 

0.070 

0.071 

0.072 

0.073 

0.073 

0.074 

0.075 

0.075 

0.076 

0.041 

0.042 

0.043 

0.043 

0.045 

0.061 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.070 

0.071 

0.071 

0.072 

0.072 

0.073 

0.042 

0.043 

0.044 

0.044 

0.046 

0.058 

0.061 

0.062 

0.064 

0.065 

0.066 

0.066 

0.067 

0.068 

0.069 

0.069 

0.070 

0.070 

0.044 

0.044 

0.045 

0.045 

0.046 

0.049 

0.059 

0.061 

0.063 

0.064 

0.065 

0.066 

0.067 

0.067 

0.068 

0.069 

0.069 

0.070 

0.045 

0.045 

0.046 

0.046 

0.047 

0.049 

0.054 

0.059 

0.062 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.069 

0.070 

0.046 

0.046 

0.047 

0.047 

0.048 

0.050 

0.052 

0.057 

0.060 

0.062 

0.063 

0.065 

0.066 

0.066 

0.067 

0.068 

0.069 

0.069 

0.047 

0.047 

0.048 

0.048 

0.049 

0.050 

0.052 

0.055 

0.059 

0.061 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.069 

0.048 

0.048 

0.049 

0.049 

0.050 

0.051 

0.052 

0.055 

0.058 

0.060 

0.062 

0.063 

0.065 

0.066 

0.066 

0.067 

0.068 

0.069 

0.049 

0.049 

0.050 

0.050 

0.051 

0.052 

0.053 

0.055 

0.057 

0.059 

0.061 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.050 

0.050 

0.051 

0.051 

0.052 

0.053 

0.054 

0.055 

0.057 

0.059 

0.061 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.051 

0.051 

0.052 

0.052 

0.053 

0.054 

0.055 

0.056 

0.057 

0.059 

0.061 

0.062 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.052 

0.052 

0.053 

0.053 

0.054 

0.055 

0.056 

0.057 

0.058 

0.059 

0.061 

0.062 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.053 

0.053 

0.054 

0.054 

0.055 

0.056 

0.057 

0.057 

0.059 

0.060 

0.061 

0.062 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.054 

0.054 

0.055 

0.055 

0.056 

0.057 

0.057 

0.058 

0.059 

0.060 

0.062 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.055 

0.055 

0.056 

0.057 

0.057 

0.058 

0.058 

0.059 

0.060 

0.061 

0.062 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.056 

0.057 

0.057 

0.058 

0.058 

0.059 

0.059 

0.060 

0.061 

0.062 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.069 

0.057 

0.058 

0.058 

0.059 

0.059 

0.060 

0.060 

0.061 

0.062 

0.063 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.069 

0.069 

0.058 

0.059 

0.059 

0.060 

0.060 

0.061 

0.061 

0.062 

0.063 

0.063 

0.064 

0.065 

0.066 

0.067 

0.068 

0.068 

0.069 

0.070 

0.059 

0.060 

0.060 

0.061 

0.061 

0.062 

0.062 

0.063 

0.063 

0.064 

0.065 

0.066 

0.067 

0.067 

0.068 

0.069 

0.070 

0.070 

0.060 

0.061 

0.061 

0.062 

0.062 

0.063 

0.063 

0.064 

0.064 

0.065 

0.066 

0.066 

0.067 

0.068 

0.069 

0.069 

0.070 

0.071 

0.061 

0.062 

0.062 

0.063 

0.063 

0.064 

0.064 

0.065 

0.065 

0.066 

0.067 

0.067 

0.068 

0.069 

0.069 

0.070 

0.071 

0.071 

0.062 

0.063 

0.063 

0.064 

0.064 

0.065 

0.065 

0.066 

0.066 

0.067 

0.067 

0.068 

0.069 

0.069 

0.070 

0.071 

0.071 

0.072 

0.063 

0.064 

0.064 

0.065 

0.065 

0.066 

0.066 

0.067 

0.067 

0.068 

0.068 

0.069 

0.070 

0.070 

0.071 

0.071 

0.072 

0.073 

0.065 

0.065 

0.065 

0.066 

0.066 

0.067 

0.067 

0.068 

0.068 

0.069 

0.069 

0.070 

0.070 

0.071 

0.072 

0.072 

0.073 

0.073 

0.066 

0.066 

0.066 

0.067 

0.067 

0.068 

0.068 

0.069 

0.069 

0.070 

0.070 

0.071 

0.071 

0.072 

0.072 

0.073 

0.073 

0.074 

0.067 

0.067 

0.067 

0.068 

0.068 

0.069 

0.069 

0.070 

0.070 

0.070 

0.071 

0.072 

0.072 

0.073 

0.073 

0.074 

0.074 

0.075 

0.068 

0.068 

0.068 

0.069 

0.069 

0.070 

0.070 

0.070 

0.071 

0.071 

0.072 

0.072 

0.073 

0.073 

0.074 

0.075 

0.075 

0.076 

0.069 

0.069 

0.069 

0.070 

0.070 

0.071 

0.071 

0.071 

0.072 

0.072 

0.073 

0.073 

0.074 

0.074 

0.075 

0.075 

0.076 

0.076 

0.070 

0.070 

0.070 

0.071 

0.071 

0.072 

0.072 

0.072 

0.073 

0.073 

0.074 

0.074 

0.075 

0.075 

0.076 

0.076 

0.077 

0.077 

0.071 

0.071 

0.071 

0.072 

0.072 

0.073 

0.073 

0.073 

0.074 

0.074 

0.075 

0.075 

0.076 

0.076 

0.077 

0.077 

0.077 

0.078 

0.072 

0.072 

0.072 

0.073 

0.073 

0.074 

0.074 

0.074 

0.075 

0.075 

0.076 

0.076 

0.077 

0.077 

0.077 

0.078 

0.078 

0.079 

0.073 

0.073 

0.073 

0.074 

0.074 

0.075 

0.075 

0.075 

0.076 

0.076 

0.077 

0.077 

0.077 

0.078 

0.078 

0.079 

0.079 

0.080 

0.074 

0.074 

0.074 

0.075 

0.075 

0.076 

0.076 

0.076 

0.077 

0.077 

0.078 

0.078 

0.078 

0.079 

0.079 

0.080 

0.080 

0.080 

0.075 

0.075 

0.075 

0.076 

0.076 

0.077 

0.077 

0.077 

0.078 

0.078 

0.078 

0.079 

0.079 

0.080 

0.080 

0.080 

0.081 

0.081 

0.076 

0.076 

0.076 

0.077 

0.077 

0.078 

0.078 

0.078 

0.079 

0.079 

0.079 

0.080 

0.080 

0.081 

0.081 

0.081 

0.082 

0.082 
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AFRL  SUPERTRAPP  JP-8  Surrogate 
Density  [kg/m3] 
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AFRL  SUPERTRAPP  JP-8  Surrogate 
Dynamic  Viscosity  [(N-s)/(m2)]*106 
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Appendix  B.  Heat  Exchanger  Design  Calculations 
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Constants  and  Unit  Conversions 


bar  :=  14.503773773021  psi 


R  :=  8.314 


kJ 

mol-K 


kJ  :=  1000J  kmol  :=  mol 


a  :=  5.67-10' 


■  8  W 


m2-K4 


Assumed  Heat  Exchanger  Inlet  and  Outlet  Temperatures 

Tfluidin  -=  290K  Tfluidout  =  530K 

Ambient  Temperature 

Tamb:=  290K 

Inner  Tube  Dimensions 

djj  :=  2.067in  dio  :=  2.375in  1  :=  30in 

Outer  Tube  Dimensions 

dQj  :=  2.635in  dQ0  :=  2.875in 

Fuel  Properties 

Choose  your  detonation  fuel 
For  Hydrogen  FUEL=1 

For  JP8  FUEL=2 

FUEL  :=  2 

PDE  Operating  Parameters 

frequency  :=  15Hz  Voltube  :=  245 in3  FF  :=  1 
Tubes  :=  2  Patm  :=  latm  Tmjx  :=  394K  <|)  :=  1.05 


Heat  Transfer  Parameters 
Turbulent  Reynolds  Number 

Returb :=  2300 
Laminar  Nusselt  Number 

Nulam  :=  3  66 


MJ  :=  106J 
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Programming  Parameters 


dl  :=  ,001m  qtQ,  :=  .1 W  dT  :=  .IK 

For  Turbulent  Flow  Effects  when  Returb  =  2300  Set  TURB  to  1 


For  Purely  Laminar  Flow  Effects  Set  TURB  to  0 
TURB  :=  1 

For  a  Non-Adiabatic  Outer  Wall  (includes  Free  Convection  and  Radiation)  Set  FCR  to  1. 
For  an  Adiabatic  Outer  Wall  Set  FCR  to  0. 

FCR:=  1 

Choose  Your  Cooling  Fluid 
For  Water  FLUID=1 

For  Nitrogen  FLUID=2 

For  JP8  Surrogate  FLU1D=3 

For JP8  FLUID=4 

FLUID  :=  3 
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Stainless  Steel  Material  Properties 

From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt, 
John  Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.1 


From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt, 
John  Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.8  (Lightly  Oxidized  Stainless  Steel) 
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Water  Properties 

Heat  of  Formation  for  Water.  From  Combution  (Third  Edition),  Irvin  Glassman,  Academic 
Press,  Copyright  1997,  Chap  1,  Table  1,  pg  6. 


H0Fwater:=  13.44 


kj 

gm 


From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt,  John 
Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.6 


water  • 


'273.15' 

f  1 

'4.217' 

\ 

275 

1 

4.211 

280 

1 

4.198 

285 

1 

4.189 

290 

1.001 

4.184 

295 

1.002 

4.181 

300 

1.003 

4.179 

305 

1.005 

4.178 

310 

1.007 

4.178 

315 

1.009 

4.179 

320 

1.011 

4.180 

325 

1.013 

4.182 

330 

1.016 

4.184 

335 

1.018 

4.186 

340 

1.021 

4.188 

345 

1.024 

4.191 

350 

1.027 

4.195 

355 

1.030 

4.199 

360 

1.034 

4.203 

365 

1.038 

4.209 

370 

1.041 

4.214 

373.15 

1.044 

4.217 

375 

1.045 

4.220 

380 

1.049 

4.229 

385 

1.053 

4.232 

390 

1.058 

4.239 

400 

1.067 

3 

4.256 

410 

K  u  x  ■= 

1  077 

io- 3 

m  1 

A  OTO 

kj 

uwater  ■ 

1  .V  /  / 

.  r  water  • 

k§  ^watpr 

cpwater : 

4.Z/0 

kg’K 

105 


420 

l.USS 

”  U*,VI 

4.3U2 

430 

1.099 

4.331 

440 

1.110 

4.36 

450 

1.123 

4.40 

460 

1.137 

4.44 

470 

1.152 

4.48 

480 

1.167 

4.53 

490 

1.184 

4.59 

500 

1.203 

4.66 

510 

1.222 

4.74 

520 

1.244 

4.84 

530 

1.268 

4.95 

540 

1.294 

5.08 

550 

1.323 

5.24 

560 

1.355 

5.43 

570 

1.392 

5.68 

580 

1.433 

6.00 

590 

1.482 

6.41 

600 

1.541 

7.00 

610 

1.612 

7.85 

620 

1.705 

9.35 

625 

1.778 

10.6 

630 

1.856 

12.6 

635 

1.935 

16.4 

640 

2.075 

26 

645 

2.351 

90 

,  647.3  j 

i  3.170  . 

v  200  ; 

106 


107 


108 

.84 

621 

104 

.85 

608 

101 

.86 

594 

97 

.87 

580 

94 

.90 

563 

91 

.94 

548 

88 

.99 

528 

84 

1.05 

513 

81 

1.14 

497 

77 

1.30 

467 

72 

1.52 

444 

70 

1.65 

67 

2.0 

412 

64 

2.7 

392 

59 

4.2 

367 

54 

12 

331 

45  j 

v  ioo  ; 

,238; 
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Air  Properties 


From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt  John 
Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.4 


'  100 

\ 

'  2.00 

\ 

'9.34 

150 

4.426 

13.8 

200 

7.590 

18.1 

250 

11.44 

22.3 

300 

15.89 

26.3 

350 

20.92 

30.0 

400 

26.41 

33.8 

450 

32.9 

37.3 

500 

38.79 

40.7 

550 

45.57 

43.9 

600 

52.69 

46.9 

650 

60.21 

,  2 

49.7 

700 

•K  vair:= 

68.10 

•!0-  SL  ^ 

52.4 

750 

76.37 

S 

54.9 

800 

84.93 

57.3 

850 

93.8 

59.6 

900 

102.9 

62 

950 

112.2 

64.3 

1000 

121.9 

66.7 

1100 

141.8 

71.5 

1200 

162.9 

76.3 

1300 

185.1 

82 

1400 

213 

91 

1500 

240 

100 

1600, 

V 

268  j 

V 

106, 

109 


^2.54^ 

'.786n 

5.84 

.758 

10.3 

.737 

15.9 

.720 

22.5 

.707 

29.9 

.7 

38.3 

.69 

47.2 

.686 

56.7 

.684 

66.7 

.683 

76.9 

.685 

87.3 

2 

.690 
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‘10_6"T  Prair:" 

.695 

109 

.702 

120 

.709 

131 

.716 

143 

.720 

155 

.723 

168 

.726 

195 

.728 

224 

.728 

238 

.719 

303 

.703 

350 

.685 

<390  j 

^.688  , 

Air  density  and  expansion  coefficent  for  Ideal  Gas 


„  „„„  J  ,  ^air 

R  :=  287 — —  P  :r  :=  latm  pair  := - 

air  kg-K  air  air  R.;^T 


air 
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Tables  of  Thermodynamic  and  Transport  Properties  of  Air,  Argon,  Carbon  Dioxide,  Carbon 
Monoxide,  Hydrogen,  Nitrogen,  Oxygen,  and  Steam;  Hilsenrath,  Hoge,  Beckett,  Masi, 
Benedict,  Nuttall,  Fano,  Touloukian,  Woolley;  Pergamon  Press  Inc.;  1960 


CP' 


^250 

\ 

^79.528 

f  4.124 

:=  5.985-10  2 — 

260 

75.689 

4.056 

hr-m 

270 

72.254 

4.000 

280 

69.172 

3.953 

:=  1 .25046- 10_  3.— 

290 

66.376 

3.913 

cm 

300 

63.810 

3.878 

I 

310 

61.475 

3.848 

:=  .296774 - 

gm-K 

320 

59.311 

3.822 

330 

57.320 

3.799 

_  -4  W 

340 

55.470 

3.779 

=  2.41-10  — — 

350 

cm-K 

53.747 

3.7619 

360 

52.141 

3.7466 

370 

50.632 

3.7331 

380 

49.213 

3.7215 

390 

47.881 

3.7113 

400 

46.625 

3.7023 

410 

45.434 

■ 

3.6948 

TN2  := 

420 

•K  Pn2:= 

44.309 

•convp  cpN2  := 

3.6882 

430 

43.215 

3.6829 

440 

42.229 

3.6784 

450 

41.267 

3.6746 

460 

40.346 

3.6720 

470 

39.473 

3.6698 

480 

38.631 

• 

3.6685 

490 

37.832 

3.6679 

500 

37.065 

3.6680 

510 

36.327 

3.6685 

520 

35.625 

3.6695 

530 

34.950 

3.6712 

540 

34.296 

3.6733 

550 

33.669 

3.6759 

560 

33.068 

3.6788 

570 

32.488 

3.6823 

580 

31.927 

3.6859 

V 

590; 

V 

31.386  , 

V 

3.6900 ; 

conv, 


cp 
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^N60  := 


^N80  :~ 


(iN2 ' 


( 1.127^ 
1.248 
1.362 
1.472 
1.578 
1.680 
1.778 
1.868 
^1.955 j 

f  1.154^ 
1.269 
1.380 
1.487 
1.591 
1.692 
1.789 
1.879 
^  1.965 j 


( 300^1 
350 
400 
450 
500 
550 
600 
650 
1^700  J 


•K 


*N2  ■ 


^N60  +  'UN80 

^N2  := - - - convri 

length(TN2)  =  35 

i  :=  0..34 

^intN2.:=  linterP^T(jN2’^N2,TN2 


(  .924  'N 
.957 
.990 
1.021 
1.051 
1.081 
1.111 
1.141 
1.172 
1.202 
1.232 
1.262 
1.292 
1.321 
1.349 
1.377 
1.405 
1.433 
1.460 
1.487 
1.513 
1.540 
1.566 
1.592 
1.619 
1.645 
1.671 
1.697 
1.722 
1.747 
1.771 
1.795 
1.819 
1.843 
1 .867 J 


•convi. 
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AFRL  SUPERTRAPP  JP-8  Surrogate  Properties 
Temperature  Range  for  data 

^suriow  :=  ^73 K  ^surhigh :=  998K 


i  :=  0.. 


^surhigh  ^surlow 


K 


Tsur;  ’  TsurIow  + 


Fuel  Pressure  for  Test 


Pressure  :=  round  f 700psi  >  q 
V  atm 


Assinging  Surrogate  Thermodynamic  Properties 
Psurjp8 


...\mujp8.xls  ■usur.  •  Psurjp8 


N-s 


i,  Pressure- 1  2 

m 


Psurjp8  * 


...\rojp8.xls  p  sur{  •  Psurjp8j 


kg 

.Pressure-1  3 
m 


*%urjp8 


W 


...\kjp8.xls  ksuri:  ksurjp8i>pressure_1-m.K 


Cpsurjp8  .  kJ 

..Acpjp8.XlS  V»V-  cp-r|p8,iPressure_ 
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From  U.S.  Oil  and  Refining  Co.,  https://www.usor.com/pdfs/specs/lpd/finished/JP-8.pdf 

HOCing  :=  42.8  — 

Jp8  k§ 

AFRL  Provided  PPDS  JP-8  Surrogate  Data  @  1000  psia 


*jp8  := 


,t  > 

^51.963^ 

^.436^ 

100 

49.438 

.499 

200 

46.798 

.563 

300 

43.996 

.626 

400 

40.953 

.689 

500 

5 

37.521 

lb 

.752 

600 

+  459.67 

7K  PjP»:= 

33.352 

'~T  cpjp8 :" 
ft 

.815 

700 

27.108 

.876 

800 

27.811 

’  .821 

900 

22.672 

.873 

1000 

17.557 

.891 

1 100  J 

^  13.995  j 

<•89, 

kjP8:= 


.U/Ul 

.0706 

2.303 

.0651 

1.278 

.0596 

.845 

0.0541 

.619 

.0486 

BTU 

.479 

.0431 

5  •” 

.378 

ft-hr— K 

.0376 

9 

.282 

.041 

.13 

.0408 

.0995 

.0412 

.0766 

„  .0429  , 

v0644 ) 

lb 

ft-hr 


BTU 

lb— -K 
9 
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Assign  Cooling  Fluid  Properties 


Z  := 


if  FLUID  =  I 

q-l 

^(o)  .  1  water 

A 

K 

I1  water 

N-s 

2 
m 

Kvater 
W 


,<1> 


,<2) 


Z'^  4- 


m-K 


pwater 


<4> 

Z  4- 


kJ 

kg-K 

Pwater 


kg 


m 


if  FLUID  =  2 


Z®  4— 


XN2 


Z'l)  4- 


K 

PintN2 


,<2) 


,<3> 


,<4> 


N-s 


m 


kN2 


4- 

w 

m-K 

cpN2 

kJ 

kg-K 

PN2 

kg 

m 


if  FLUID  =  3 


z<». 


sur 


,<i> 


K 

Psur 

M.c 
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,<2> 


,<3> 


,<4> 


*  2 
m 

<— 

ksur 

W 

m-K 

cpsur 

^ — 

kJ 

kg-K 

Psur 

kg 

m 


if  FLUID  =  4 


<0>  ,  Tjp8 


<— 


K 


N-s 


7<2>  .  kjp8 

L  v” 


,<3> 


,<4> 


W 

m-K 

CpjP« 

kJ 


kg-K 

pjp8 

Jcg 

3 

m 


Tfluid:=Z<0>.K 

7<4>  kg 
P  fluid  •=  z  ■“ 
m 

7<1>  N-s 
Pfluid  :=  z  — ~ 

m 

.  „<2>  W 

kfluid-Z  m-K 

„  7<3>  kJ 

pfIuid '  'kg-K 


j:=0..last(pfluid)-2  dPj  := 


pj:=— dpj  Tp:=Tfluid 
PP.  J 


Pfluidj+1  Pfluidj 


1  fluid: 


j+1 


‘fluidj 


PPj 1  Pfluidj 
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Mass  Flow  Calculator 


C  := 


H:= 


0  if  FUEL  =  1 

10.9  if  FUEL  =  2 

2  if  FUEL  =  1 

20.9  if  FUEL  *  2 


Fuel, 


moles 


:=  lmol  if  FUEL  =  1 
lmol  if  FUEL  =  2 


Stoichiometric  Calculations 


CO?  :=  C  Water  :=  — 
Z  2 


Water  =  10.45 


Water  ^ 


Molesajr  :=  |  C02  +  — - — Jmol 


mwfuel (c>  12-01'~ ]  +  1  '01  mfuel :=  FueImoles-mwfuel 


32--^-  +  —.28.02--^ 


mwair := 


mol  21 


mol 


l+2 

21 


ka  79  ks 

m,ir  :=  Moles,;,.-!  32 — —  + - 28.02 — - 

mol  21  mol 


m, 


air 


=  14.575 


mfuel 


mfuel 


=  0.068609 


m 


air 


Mixture  Calculations  for  Operating  Paramaters 

ci-  ,  mfuel  1 

fueltoair  :=  (p - airtofuel  := 


air 


fueltoair 
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Volumetric  and  Mass  Flow  Rate  Calculations 

V'dot  :=  frequency' Voltu[,e  FF-Tubcs 


lb 


min 


madot '  ^  ’ 

Given 

Ydot  “  Yadot  +  Yfdot 

V  P 
vadot  ratm 


m. 


adot 


R 


mix 


lb 

mfdot ' 

min 


Vadot  ■  Ydot 


V, 


mfdot  “ 


Vfdot'Patm 


R 


mw, 


fuel 


-T  • 

mix 


fd°t '  1Q0  'Vdot 


m 


fdot 


=  fueltoair 


madot 

Y  :=  Find(Vadot ,  Vfyot,  macjot ,  mfy0t) 

\r  y 

adot  •  o 
V»*:=YI 
madot ,=  Y~ 


mfdot Y~ 


mdot 1  mfdot 


mdottot :=  madot  +  mfdot 
79 

Nitrogen  :=  Molesa.-r — 


21 


Water-mol- 18- 


mdotwater '  mdottot 


kg 

mol 


f  Water-mol- 1 8-—^  j  +  f  Nitrogen- 28. 02-—^- 
\  mol  J  ^  mol 


mdotnitrogen  •  mdottot' 


kg 

Nitrogen- 28.02 - 

mol 


Water-mol- 18-—^-  )  +  Nitrogen- 28.02--^- 


V 


mol 


mol 
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Calculation  of  specific  energy  and  power  required  for  inlet  and  outlet  temperatures 


Area  Under  Cp  vs.  T  curve  is  the  required  specific  energy  required  to  raise  the  cooling  fluid 
temperature  across  the  required  range. 


n  :=  10000 


dt:= 


Tfluidout  ^fluidin 


n-  1 


dt  =  0.024  K 


i  :=  0..  n  -  1 


T  := 
er.  ’ 


j Tfluidip  +  i-dt  cPj linterP^Tjp8’cpjp8’Ter^ 
n-1 

SpecificEnergy  :=  ^  ^Cp-dt  EnergyRequired  :=  m^-SpecificEnergy 


i  =  0 


EnergyReleased  := 


PercentRequired  := 


mdotwater ' ^^^water  FUEL  1 

mdot'HOCjp8  ^  FUEL  =  2 

EnergyRequired 

EnergyReleased 
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Calculating  average  heat  transfer  coefficient  and  temperature  inside  the  PDE  tube  during  test 
From  Steady  State  Tube  Free  Convection  and  Radiation  Test 

Tsurf:=  1005K 

Tinf:=  290K 

Film  Temperature 
T  "I-  T* 

T  surf - mf  eq  7.2,  pg  326,  Introduction  to  Heat  Transfer  (Third  Edition) 

1  2 

Tube  Dimensions 
D0  :=  2.375in 
Dj  :=  2.157in 

Free  Convection  Heat  T ransfer  Coefficient  Calculation 


g-linterf 

’(Tair’Pair>Tf) 

(Tsurf  Tinf) 

)*o 

3 

EQ  9.24,  pg  456,  Introduction  to 

linterp^ 

Tair’vair’Tf)-HnterPl 

■  oc  • 
(airman 

r>Tf] 

l~  Heat  Transfer  (Third  Edition) 

Ra  =  4.58  x  IQ5 


Nu  := 


.6  + 


,387-Ra 


1  + 


.559 


9 

16 


8 

n27 


linterp(Tajr,Prair,T^ 


EQ  9.27,  pg  457,  Introduction  to 
Heat  Transfer  (Third  Edition) 


Nu  =  11.658 


h  := 


linterp(Tair,kair>Tf)'Nu 


h  =  9.577 


kg_ 

s3K 


EQ  6.58,  pg  307,  Introduction  to 
Heat  Transfer  (Third  Edition) 
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Heat  Transfer  from  Surface  of  Tube  Due  to  Free  Convection  and  Radiation 
General  Convective  Heat  Transfer  ,fc  =  [hfc ■«  <Baouter'(Tsurf  -  T„mb)] 

General  Radiative  Heat  Transfer  qrad  =  sodiWr<.(Tsurf4  -Tamb4) 
q  :=  [l>  >  D0  (Tsurf  -  Tinf)]  +  lmte^,(Te,sss,Tf).rt.D0  e(Tsurf4  -  Tinf4) 


3  W 

q  =  4.143  x  10  - 
rn 

Solving  For  Inner  Tube  Wall  Temperature 

Guess  Value 
Tjs  :=  1200K 

Given 


2'7r'linterp(Tss,kss,  Tsurf)-(Tsurf  Tjs) 


^  2  ; 

^iwfcrad  ,=  ^*nc^(^is) 

Inner  Wall  Temperature 

TiWfcrad  =  1-007x  1q3k 


EQ  3.27,  pg  91 ,  Introduction  to  Heat 
Transfer  (Third  Edition) 


Heat  Flux  Calculation  for  Free  Convection  and  Radiation 
fluxfcrad  :==  [h-(Tsurf  ~  Tinf)]  +  linterP(Te  >ess’Tf)'CT'(Tsurf  “  Tinf  ) 
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From  Water  Jacket  Tests 


Operating  Parameters 

Vdotwater:=  6-l~  AT:=  24’f K  AT=  13’333K 
min  y 


Tube  Dimensions 

Lengthtube  :=  36in  Dj  :=  2.067in  DQ  :=  2.375in 

Heat  Transfer  to  Water 


q  =  mdot-cp-(T  -  T  )  EQ  8.37,  pg  399,  Introduction  to  Heat  Transfer  (Third 

'  x+  x  x '  Edition) 


L  295'K+ (295-K+ AT)’ 

^  Motwater*  linterP  Twater?P  water  2 


linterp 


T  c 
water  vpwaters 


295-K+  (295- 


Heat  Flux  to  Water 

q 


fluxwj  := 


Lengthtube- Dj-n 


kW 

flux  :  =  141.651 - 

wj  2 

m 


V- 


f  100  Y  (65  s] 

200 

400 

v600y 


•K 


kal := 


163 
186 
U«6 ; 


W 

m-K 


Solving  for  Inner  Wall  Temperuture 

n  -=  q 

qprime- 


Tjwwj  ==  500K 


Given 


Qprime  “ 


_  2-7t-linterp(Tal,kal,295-K+  AT)-[Tiwj  -  (29.5 -K+  At)] 


f  D,  ^ 


Ini 


vDiy 


Ywwj  •  Find(Tj%VVf,j)  Tj^^j  324,653  K 
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Two  Equations  Two  Unknowns.  Solve  for  PDE  Tube  Avg.  Heat  Transfer  Coefficient  and  Avg.  Gas 
Temperature 

Solving  for  Heat  Transfer  Coefficient  and  Avg.  Gas  Temperature 

W 

hflame:=10°—  Tflame  :=  2°0°K 
m  K 

Given 

^uxfcrad  ”  Aflame' (Aflame  “  ^iwfcrad) 

^uxwj  “  Aflame  (Tflame  ~  Tiwwj) 

Z:= 

Aflame  :=  Zo  Tflame  :=  z\ 

W 

hflame  =  175'428  —  Tflame  =  H32.111K 
m  -K 
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Finite  Difference  Method  Heat  Exchanger  Design  Calculation 
Concentric  Tube  Calculations 

W=^oi2-<io2) 

Pconc  :=  n'  (^oi  +  ^io) 

Hydraulic  Diameter  of  Concentric  Tube  Heat  Exchanger 


4- A, 


cone 


1hconc ' 


cone 


EQ  8.67,  pg  417,  Introduction  to  Heat  Transfer  (Third  Edition) 


Cooling  Fluid  Velocity  in  Concentric  Tube 
mdot 


cone ' 


(linterp(Tfiuid,  Pfluid>Tamb)- Aconc) 

Finite  Difference  Method  Heat  Transfer  Calculations 


uconc  • 


in 


X  := 


i  <-  0 

Ti  Tfluidin 
while  T.  <  f  £jujcjout 

linterP(Tfluid ,  pfluid,Ti)uconc-dhconc 

Re  ^ — - - - 

linterP(Tfluid^fluid>Ti 

if  TURB  =  0 
Nu  <-  Nulam 

Nu-  linterp^T^uid ,  kf|uid ,  T.j 

h°  * - H 

ahconc 

if  TURB  =  1 
if  Re  <  Returb 

Nu  <"  Nulam 

Nu-linterp(Tfluid)kfluid,T.) 

ho  “ 

uhconc 


otherwise 

I  f  <—  (,790-ln(Re)  -  1.64) 


f  m  linterp(Tf]uid ,  cpfiuid  >  Ti)  ‘  linterP(Tfluid  ’  ^fluid  >  Ti) 

i  ivv  ivvvi . 

8 

linterpfTfl^.kfl^T.) 

1 

,  .  0 
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.  .  ^  ....  T\.lm+o«JT..  .  .. 

2 

.  tW3 

Nu  <- 


1  +  12.7- 


fV 
8 


LV 


utllc‘H^fluid^pfluid> 1  jj’lllucIF^ 1  fluid 1  H-flujd . 1  j J 
linterp(Tfluid,kfluid)T.J 


ho<- 


if  FCR  =  0 


Too<-°-K 


Nu-linterp(T{1uid,kfluid,T,) 


3hconc 


hfc  °“ 


W 


m2-K 

Toi<-O.K 
qwi  <-0  W 
qw2<-0-W 
if  FCR  =  1 

Too  Tamb  +  dT 
m  Too  +  Tamb 


oo 


ri  .  g* linterp^Tp , P , Tf] •  (Tq0  -  Tamb)' dc 
linterp(Tair,vair,Tf).linterp(Tair,aair>Tf) 


Nu 


fc ' 


.60  + 


.387-Ra 


1  + 


.559 


9 

16 


8 

i27 


‘fc 


linterp(Tajr ,  pra{r » Tf) 
NUfc'litlterp(Tajr,  kajr , 


■*00 


qwl  hfcdoort-dl-(Too  -  Tamb)  +  linterP(T£>ess’Too)-CT-(doo-ndl)(To04  "  Tmb) 


T  ■ 
oi 


‘Jwr'n 

ro 

v  doi  J 

2-u-dl-linterp( 

.^SS’^sS’^Oo) 

+  T 


00 


^w2  ^o,Jt  doi'd^(^i  Toi) 
while  |qw2  -  qwl|  >  qtoI 

|  Too  Too  +  dT 
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TV  <— 


Ra  4 — 


T  |  T 

oo  amb 


g-linterp(Tp,P,Tf).(T00  -  Tamb)-dop3 
IinterP(Tair  ■  vair  ’  Tf) ' H 11  terp(Tair  > a  air  ’  Tf) 


Nufc  <- 


.60  + 


,387-Ra 


1  + 


.559 


_9_ 

16 


^air’^rair’Tf) 


_8_ 

27 


hfc<" 


Nufclinterp(Tair,kair,Tf) 


OO 


^wl  ^fc‘^oo-7l‘^^(^oo  ^amb)  +  *'nterP(Te ’sss’^oo)'CT'(^oo'7I'^^)'^oo  ^aml 


fA  ^ 


qwrln 


v  d»y 


2.jr-dMinterp(Tss,kss,TO0) 

<lw2  <"  h0^-doi-dl-(Ti  -  Toj) 

Tio  Tfluidin  +  dT 

qini  <-  h0-n-d£0-dl-(Ti0  -  T.) 

qinrii 


+  T, 


oo 


f  d-  > 

uio 


T;:  <r-  ' 


Vduy 


+  T: 


2  •  linterp^Tss ,  kgs,  Tj0^  ■  dl-  jc 
^in2  ^flame,Tt'dii‘d*'(^flame  -  "^ii) 

while  |qinl  -  qin2[  >  qt0l 


10 


Tio^Tio  +  dT 


qini  <"  ho^  dio  dl  (Tio  -  Ti) 

1  fdi0' 

qinrH 


T;;  <— 


v  dn  y 


+  ^ 


2-linterp(Tss,kss,Tj0^-dl'7r 

qin2  hfIame‘n‘dudl'(Tflarae  “  Tii) 
qtrans  ^irS  “  ^w2 
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m 

W 
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Fluid  Temperature 


■"p  . _  Y 

1coolingfluid  •“  A  -JV 

Inner  Tube  Inner  Surface  Temperature 

T1  . _  Y 

pdeinner  •  A 

Inner  Tube  Outer  Surface  Temperature 

t  -  y<2>  k 

ipdeouter  A 

Outer  Tube  Inner  Surface  Temperature 

T.  -=x(3)k 
dinner-  yv  ^ 


Outer  Tube  Outer  Surface  Temperature 
T„ 


<4> 

.  :=  X  -K 


oouter  • 

Heat  Transfer  to  the  cooling  fluid 

_  y<5>  W 

‘Trans  •“  A  ' 

m 

Heat  Transfer  into  the  System 

-  v(6>  Z 

Qin  '  ^ 


m 


Heat  Transfer  rejected  by  Free  Convection  and  Radiation 

•=  x<7>  — 

^  waste  ■  A  ' 

m 

_  v<8>  W 
9waste2  x  ~ 

Heat  Losses  due  to  Radiation 

-  v<9>  Z 
^rad  ‘  A 

m 

Heat  Losses  due  to  Free  Convection 

Y<10>  w 

‘tfc 1  ^  ' 

m 

Tube  Length 
Length  :=  X^^-mm 
Reynolds  Number 


Re  :=  X 


.<12> 
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Nusselt  Number 


Nu : 


„<13> 


Fluid  Heat  Transfer  Coefficient 


W 

m2K 


Free  Convection  Heat  Transfer  Coefficient  on  Outside  of  Outer  Tube 


W 

m2-K 


x  :=  length(Length)  x  =  742 

LengthRequired  :=  Length^^ 

(LengthRequired  ■  Aconc) '  **nterp(Tfiu{d .  P  fluid  >  ^amb) 

ResidenceTime  :=  - - - - 1 - - 


mdot 


FluidVelocity  := 


LengthRequired 

ResidenceTime 


Theatexchangerexit -=  ^coolingfluidxl 
9  transtotal ,=  ^.^trans'^ 

^intotal  ~ 


Percentloss  := 


^transtotal 
^  intotal 


•100 
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Heat  Transfer  (kW/m) 


Temperature  vs.  Length 


Length 
Length  (m) 


Vans 

1000 

V 

1000 

Vaste 

1000 

Vaste2 

1000 

Vd 

1000 

Vc 

1000 


Heat  Transfer  vs.  Length 
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Reynolds  Number  Heat  Transfer  Coefficient  (W/(mA2*K)) 


Heat  Transfer  Coefficient  vs.  Length 


Length  (m) 


Reynolds  Number  vs.  Length 
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Final  Results 


Stochiometric  Values 
Moles„ir  =  16.125  mol 

Stochiometric  Fuel  and  Air  Molecular  Weight 
mwfci.152.018-^  mwair  =  28.856-^ 

Mixture  Values  for  Operating  Parameters 

fueltoair  =  0.07204  airtofuel  =  13.881 

Volumetric  and  Mass  Flow  Rate  Values 

3  3  _  3 

Vdot  =  (m~  Vadot  =  0.119—  Vfdot=  1.625  x  10“  — 
s  s  s 


madot  =  °-1067 

Energy  Required  to  Raise  Cooling  Fluid  Temperature 
From  Tfluidin  =  290  K  To  Tfluidout  =  530K 
kJ 

SpecificEnergy  =  615.163  —  EnergyRequired  =  4.7kW 
kg 

EnergyReleased  =  327  kW  PercentRequired  =  1.437 
Heat  Transfer  Coefficient  and  Inner  Tube  Temperature  Values 

Aflame  =  175.428  ~ T~  Tflame  =  1132.11 1 K 

m  K 

Final  Values 
^transtotal  =  4-226  kW 
^intotal  =  4.462  kW 
Percentloss  =  5.29 

^heatexchangerexit  =  529.928  K 

LengthRequired  =  29  in 
ResidenceTime  =  52.901  s 


132 


Appendix  C.  Plain  Tube  Heat  Transfer  Calculations 
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Constants  and  Unit  Conversions 


bar  :=.  14.503773773021psi 


kJ :=  1000J 


MJ  :=  10  J 


MPa  :=  10  Pa 


kmol  :=  mol 


R  :=  8.314 — — —  c»:=5.67'10  8  W 
mol-K  2  ..4 

m  -K 


Material  and  Air  Properties 

From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt,  John 
Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.1 


From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt,  John 
Wiley  &  Sons,  Inc.  Copyright  1996 

Appendix  A  Table  A.8  (Lightly  Oxidized  Stainless  Steel) 
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From  Introduction  to  Heat  Transfer  (Third  Edition),  Frank  P.  Incropera,  David  P.  DeWitt,  John 
Wiley  &  Sons,  Inc.  Copyright  1996 


Appendix  A  Table  A.4 


f  100  N 

f  2.00  N 

r9.34N 

150 

4.426 

13.8 

200 

7.590 

18.1 

250 

11.44 

22.3 

300 

15.89 

26.3 

350 

20.92 

30.0 

400 

26.41 

33.8 

450 

32.9 

37.3 

500 

38.79 

40.7 

550 

45.57 

43.9 

600 

52.69 

46.9 

650 

60.21 

2 

49.7 

700 

•K  vair := 

68.10 

•10"6i7  “air:= 

52.4 

750 

76.37 

54.9 

800 

84.93 

57.3 

850 

93.8 

59.6 

900 

102.9 

62 

950 

112.2 

64.3 

1000 

121.9 

66.7 

1100 

141.8 

71.5 

1200 

162.9 

76.3 

1300 

185.1 

82 

1400 

213 

91 

1500 

240 

100 

^1600 ) 

,  268  , 

^  106  y 

-3  W 
m-K 
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5-1-20  Avgas  Clean  Tube  Heat  Transfer  Calcs 
Steady-state  data  obtained  from  testing 


§  := 


r.95^ 
1.02 
1.09 
1.12 
1.22 
VI -26  J 


loc  := 


f  19.685  ^ 
46.99 
73.8188 
105.886 
121.761 
137.001 
154.623 
v  169.863  j 


•cm 


Temp  := 


^596.51 

605.4541667 

601.3773889 

595.3083889 

575.4475 

573.50811  lP 

692.76 

705.947 

669.0364444 

644.0752778 

613.7165 

618.1916111 

763.06 

781.1624444 

699.8215 

691.7939444 

648.9514444 

641.4213333 

802.92 

816.6038889 

800.0493889 

796.4642222 

746.0194444 

733.0732222 

835.58 

849.4266667 

853.8216667 

846.1422222 

798.5135556 

789.2047222 

808.82 

832.4044444 

851.9672222 

846.1422222 

788.8633889 

780.444 

842.04 

869.0783333 

878.4661111 

876.4894444 

803.4879444 

789.4714444 

^814.11 

835.9805556 

854.2244444 

854.985 

756.0436667 

738.1221667  , 

Outer  Tube  Diameter 
Dq  :=  2.375in 

Ambient  Air  Temperature 
Tinf:=  290K 


•K 
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Heat  transfer  calculation  based  on  radiation  and  free  convection  from  a  long  horizontal  cyclinder 


X  := 


for  ye  0..5 
for  xe  0..  7 


Tsurf  TemP: 


x,y 

Tsurf  +  Tinf 


Ra  4 — 


g4inter{ 

’(*air’  Pair’^f)' 

(Psurf  Pinf) 

|-Do3 

linterp^ 

Tair>vair-Tf)‘linterPl 

f-r 

T  •  ct  • 
{  air’  an 

r.Tf) 

Nu<- 


.6  + 


,387-Ra 


~\2 


1  H~ 


.559 


linterP(Tair’Prair’Tf) 
linterp(Tair,kair,Tf)-Nu 


9 

ye 

)  J 


_8_ 

27 


qx,y  [h'^'Do'(Tsurf  "  Tinf)]  +  linterp(Te,ess,Tf)^-D0^{TSurf  "  Tinf 


326) 


456) 


465) 


465) 


-10) 


Heat  Transfer  Values 


0.803 

0.839 

0.823 

0.798 

0.721 

0.714^ 

1.251 

1.323 

1.129 

1.008 

0.874 

0.893 

1.667 

1.789 

1.289 

1.246 

1.031 

0.996 

1.944 

2.047 

1.923 

1.897 

1.558 

1.479 

kW 

2.196 

2.31 

2.348 

2.283 

1.912 

1.845 

m 

1.988 

2.17 

2.332 

2.283 

1.842 

1.784 

2.249 

2.481 

2.566 

2.548 

1.948 

1.847 

2.028 

2.199 

2.351 

2.357 

1.622 

1.51  , 

„  v<0>  n  Jl) 

^primel  ■“  A  "prime2  A 


clprimc3  • 


:=  X 


-<2> 


x<3>  y<4>  y<5> 

1prime4  ■  A  ^prime5  •  A  ^primefi  ■  A 
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5-1-21  JP-8  Clean  Tube  Heat  Transfer  Calcs 
Steady-state  data  obtained  from  testing 


f  19.685  A 

'l.05N 

46.99 

73.8188 

1.07 

105.886 

1.08 

loc  := 

121.761 

1.19 

,1.26y 

137.001 

154.623 

J  69.863  ; 

Temp  := 


618.198 

610.3761667 

610.3451667 

607.8732778 

602.6850556^ 

764.2152222 

701.2448333 

751.4777222 

676.3736111 

746.2182222 

838.235 

809.0104444 

819.2511111 

821.0644444 

777.8316111 

885.5205556 

866.5633333 

869.2855556 

874.6472222 

864.8983333 

903.9816667 

878.685 

887.7288889 

896.0127778 

894.4655556 

882.1638889 

857.4505556 

868.9688889 

883.9744444 

883.1605556 

918.24 

894.0527778 

902.6311111 

911.6044444 

907.5194444 

858.1288889 

821.6572222 

838.5005556 

850.1838889 

842.4988889  ; 

Outer  Tube  Diameter 
D0  :=  2.375in 

Ambient  Air  Temperature 

Tinf  :=  290K 
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From  Steady  State  Tube  Free  Convection  and  Radiation  Test 

Y:=  |  for  ye  0..4 
for  x  e  0 ..  7 


Tsurf  TemP. 


Tf  <— 


Ra 


x,y 

Tsurf  +  Tinf 


{Incropera  and  Dewitt,  1996:326) 


gTinterp(Tair,Pair,Tf)'(Tsurf  -  T{nf)-D03 

Iinterp(Tair,vair,Tf)  linterp(Tair,aair)Tf)  (Incropera  and  Dewitt,  1 996:456) 

2 


Nu 


.6  + 


,387-Ra 


1  + 


.559 


9 

16 


_8_ 
n  27 


linterP(Tair-Prair>Tf)> 
lintcrp(Tajr,kair,Tf)-Nu 


(Incropera  and  Dewitt,  1996:465) 


(Incropera  and  Dewitt,  1996:465) 


D„ 


qx,y  [h  "  D0-(Tsurf  ~  Tinf)]  +  linterp(Te  ^ss’^^-Do'^  ^surf4  "  Tinf4) 

(Incropera  and  Dewitt,  1996: 8-10) 


Heat  Transfer  Values 


'0.893  0.86  0.86  0.849 

1.675  1.297  1.592  1.166 
2.217  1.989  2.067  2.081 
2.631  2.458  2.483  2.531 
2.809  2.568  2.652  2.731 
2.6  2.379  2.48  2.617 

2.963  2.712  2.795  2.887 
^2.385  2.085  2.22  2.317 


0.828N 

1.559 

1.766 

2.444 

2.716 

2.609 

2.844 

2.252 , 


kW 

m 


^primel  ■ 


9prime2 


9prime3 


9pnme4  ’ 


9prime5  * 
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Appendix  D.  Sample  Wavespeed  Calculation 
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Constants  and  Unit  Conversions 

bar  :=  14.503773773021psi  kj  :=  1000J  101101  :=  mo1  MJ  :=  106J 

kJ  -  8  W  1  1 

R:=  8.314 -  a  :=  5.67-10  -  ms:= - s  ps:= - ms 

mol-K  2  4  1000  1000 

m  K 

Assigning  Data 
WSIG  := 

..\8.xls 

datapoints  :=  9804 
rows  :=  9804 

shift  :=  11 
colshifts :=  0 

ion  probe  voltage  tolerance  for  wave  time 
iontol :=  4-45 

Spark  derivative  tolerance  for  spark  time 

sParkdertol  :=  -J5 

High  and  Low  Tolerances  for  spark  time 
sParkt0ll0W  := 
sParktolhigh  :=  3-5 
skip  :=  200 

res  :=  [(wSIG'°^)1  -  (wSIG^)0  -s 

Data  Resolution 

res  =  51  ps 

Data  Column  Numbers 

timecol  :=  0  sparkcol  :=  1  ionlcol:=5  ion2col  :=  6  ion3col  :=  7  ion4coI 
ion5coI  :=  9  ion6col:=  10  ion7col:=  11  ion8col:=  12 


Ion  Probe  Locations 

f  7.75  ^ 
18.5 
29.0625 
41.6875 

47.9375 

53.9375 

60.875 

66.875  ) 


ionloc :: 


in 
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Assigning  Time  Data 


timevector:= 


z  <-  0 

while  z  <  datapoints 
for  xe  O..colshifts 

for  row  e  0 ..  rows  -  I 

col  <-  x- shift 

time  4-  WSIG 

z  row,  col 


z  +  1 


time-s 


i  :=  0.. datapoints  -  1 


time.  :=  t«nevecto 


Ion  Probe  Distance  Calculation 


i  :=  0..  last/ionloc)  -  1  dist.  :=  ionloc.  ,  -  ionloc. 

l  l+l  l 


dist  = 


Assigning  Spark 
sparkyec  := 


Data 

z  +-  0 

channel  <—  sparkcol 
while  z  <  datapoints 
for  xe0..colshifts 

for  row  e  0 ..  rows  -  1 

|  col  +-  x- shift  +  channel 

var  <-  WSIG 

z  row,  col 

z  <-  z  +  1 


var 


i  :=  0..  datapoints  -  1 
spark.  :=  sparky^ 

i  :=  1 ..  last(time) 


(  10.75  N 
10.563 
12.625 
6.25 
6 

6.938 

V  6  j 
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Spark  Derivative 

sParkder. := 


spark.  -  spark._j 


time.  -  time.  . 
1  i-l 


sparkder  := 


sparkd( 


er 


maxi 


Assigning  Ion  Probe  Data 


ionl 


vec  ■ 


z  <-  0 

channel  <—  ionl  col 
while  z  <  datapoints 
for  xe0  ..colshifts 

for  row  e  0 ..  rows  -  1 
col  <-  x-shift  +  channel 


var  <—  WSIG 


row,  col 


z  <—  z  +  1 


var 


i  :=  0..  datapoints  -  1 

ionl.  :=  ionl™,. 
i  vec; 


ion3. 


vec  - 


z  <r-  0 

channel  ion3col 
while  z  <  datapoints 
for  xe  O..colshifts 

for  row  e  0 ..  rows  -  1 
I  col  -t—  x-shift  +  channel 


var  <-WSIG 


row,  col 


z  <—  z  +  1 


|  var 

i  :=  0 ..  datapoints  -  1 
ion3  :=  ion3yec 


(sparkder) 


ion^vec  := 


ion4. 


vec ' 


z  <r-  0 

channel  <—  ion2col 
while  z  <  datapoints 
for  x  e  0 ..  cols|^ds 

for  row  e  0 ..  rows  -  1 

col  <-  x-shift  +  channel 

var  <—  WSIG 

z  row,  col 

_  z  <—  z  +  1 

var 

i  :=  0 ..  datapoints  -  1 
ion2.  :=  ion2yeCi 
z  <-  0 

channel  <—  ion4col 
while  z  <  datapoints 
for  xe0..colshjfts 

for  row  e  0 ..  rows  -  1 
I  col  •<—  x-shift  +  channel 


var  <-  WSIG 


row,  col 


z  +  1 


var 


i  :=  0 ..  datapoints  -  1 
ion4.  :=  ion4. 


'vec. 
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ion5, 


vec  ‘ 


z  <—  0 

channel  ion5col 
while  z  <  datapoints 
for  xe  O..colshifts 

for  row e  O..rows -  I 

|  col  <-  x- shift  +  channel 

var  4-WSIG 

z  row,  col 

Z  <r-  Z  +  1 


var 


i  :=  0..  datapoints  -  1 
ionSj  :=  ion5yec 


ion6vec  := 


z  0 

channel  <-  ion6col 
while  z  <  datapoints 
for  xe  O..colshifts 

for  rowe  O..rows  -  1 

|  col  <—  x-  shift  +  channel 

var  <-  WSIG  . 
z  row,  col 


z  +  1 


var 


i  :=  0..  datapoints  -  1 
ion6.  :=  ion6yec 


ion7vec := 


z,<—  0 

channel  <—  ion7col 
while  z  <  datapoints 
for  xe0..colshifts 

for  row  e  0 ..  rows  -  1 

|  col  <—  x- shift  +  channel 

var  WSIG 

z  row,  col 

z  4-  z  +  1 


|  var 

i  :=  0 datapoints  —  1 
i°n7j  :=  ion7yec 


ion8. 


vec  • 


z  <—  0 

channel  ion8col 
while  z  <  datapoints 
for  x  e  0 ..  col 


shifts 

for  row  e  0 ..  rows  -  1 
|  col  4—  x-shift  +  channel 


var 


WSIG 


z  row,  col 

z  <-  z  +  1 


var 


i  :=  0..  datapoints  —  1 
ion8.  :=  ion8. 


vec- 
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Determining  Spark  Time 

tspark :=  x<_  0 

y  0 

while  x  <  last(time) 

if  sparky  <  sparkdertol 

if  sparky  <  sparkto|high 

^sparky  timex 

y  <-  y  +  1 
x<-  x+  1 

x  <—  x  +  1  otherwise 
x  <-  x  +  1  otherwise 

lspark 
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Determining  wave  times  at  each  Ion  Probe 


W 


ion  «-  ionl 
x  <-  0 
y  <-  0 

while  x  <  last(time) 
if  ionx  <  ionto] 

if  y  <  las<tspark) 

I  if  time  > 
x 


ti°ny 


timev-rionto!  +  timev‘ionx_i  _  timex-iontoj  -  ion^-time 


x-1 


ion  ,  -  ion 

x-1  x 


Ix  <-  x  +  skip 
y  <-y+  1 

x  x  +  1  otherwise 
x  4-  x  +  1  otherwise 
x  <—  x  +  1  otherwise 
a  0 
b  <—  0 

while  a  <  Iast(tspark) 

|lf  ti°nb<Vrka+1 


timel  t:__ 
a  ion. 


I  a  a  +  1 
b  < —  b  +  1 

otherwise 

timela<-  .OOlps 

a  <-  a  +  1 
timel 
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*10112 ' 


ion  <r-  ion2 
x  <-  0 
y  <-  0 

while  x  <  last(time) 
I  if  ionx  <  iontol 


if  y  <  last(tspark) 
if  time  >  t 


sparky 


+  5ms 


time  ,  iont.,i  +  time  -ion  .  -  time  -ionf„i  -  ion  -time  . 
X— 1  tol  x  X— 1  X  tol  x  x-1 

tiony *  ion  ,-ion 

x- 1  x 

x  < —  x  4-  skip 
y  «-  y  +  1 

x  <—  x  4-  1  otherwise 
x  <-  x  +  1  otherwise 
x  <-  x  4-  1  otherwise 
a  <-  0 
b  <-  0 

while  a  <  last(tspark) 

|if,io„b«spa*wl 


timel  <-  t:ntl 
a 


a  <—  a  +  1 


|b<-b+  1 

otherwise 

timela<—  .002jas 

a  <-  a  4-  1 
time! 


150 


lion3  : 


ion  <-  ion3 
x  <-  0 
y  4-  0 

while  x  <  last(time) 
if  ionx  <  ionto[ 

if  y  <  'as^spark) 

if  timex>tsparky 

timex-f  iontol  +  timex'ionx-l  _  timex'iontol  -  ionx'timex-l 

tion  <  :  : - 

y  ion  ion 

x-1  x 

x  4-  x  +  skip 
y  4-  y  +  1 

x  4-  x  +  1  otherwise 
x  4-  x  +  I  otherwise 
x  4-  x  +  1  otherwise 
a  4-  0 
b  -<-0 

while  a  <  last(tspark) 


if  t;„„  <t, 
b 


ionu  "  ‘■spark. 


a+1 


timel 


ion. 


a  4-  a  +  1 
b  4-  b  +  1 


otherwise 

timela  4-  0.003 ps 

a  4-  a  +  1 
timel 
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lion4  ’ 


ion  <-  ion4 
x  <-  0 
y  0 

while  x  <  last(time) 


if  ionx  <  ionto, 


if  y  <  last(^tSpar^ 
if  time  >  t, 


sparky, 

timex_j-iontoj  +  time^ion^  -  timex- ionto[  -  ion  -time  j 


ion. 


ion  ,  -  ion 
x-1  x 


X  <r-  X  4-  skip 
y  <—  y  +  1 

x  <r-  x  +  1  otherwise 
x  x  +  1  otherwise 
x  <—  x  +  1  otherwise 
a  0 
b  <—  0 

while  a  <  last(tspark) 


lf  tionb<tsParka+1 


timel  <—  t:-.n 
a  1(5nb 

a  <r-  a  +  1 

b  b  +  1 


otherwise 

timely  <—  0.004ps 

a  a  +  1 
timel 


152 


*ion5  :_ 


ion  4—  ion5 
x  4-  0 
y  4-  0 

while  x  <  last(time) 
if  ionx  <  iontd 

if  y  <  last^spark) 

if  time  > 
x 

tirnvriontoi +  timex-ionx-i -  timex'iontoi  -  ionx'timVi 

tjon  <  :  :  ~ 

y  ion  .  -  ion 

X— 1  X 

x  4-  x  +  skip 
y  4-  y  +  1 

x  4—  x  +  1  otherwise 
x  4-  x  +  1  otherwise 
x  4-  x  +  1  otherwise 
a  4-  0 
b  4-0 

while  a  <  last(tspark) 

,f  tionb  <  tsParka+l 

timela  4-  tiollb 

a  4-  a  +  1 
b  4-  b  +  1 

otherwise 

time  1  a  4-  0.005ps 

a  4-  a  +  1 
timel 
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lion6  1 


ion  <—  ion6 
x  <-  0 
y  <-  0 

while  x  <  last(time) 
if  ionx  <  iontol 

I  if  y  <  last(tspark) 


if  timex>tsparky 


time  ,-iont  i  +  time  ’ion  ,  -  time  -iontril  -  ion  -time  , 
X— 1  tol  x  X-l  X  tol  X  X—  1 


ion  ,  -  ion 

X-l  X 


Ix  <—  x  4-  skip 
y  4-  y  +  1 

x  <—  x  +  1  otherwise 
x  <r-  x  +  1  otherwise 
x  x  +  1  otherwise 
a  <-  0 
b  0 

while  a  <  last(tspark) 

|'f  tionb<tsParka+l 


timel 


a  ti°n 


b 


a  <-  a  +  1 
b  <-  b  +  1 

otherwise 

timely  <—  0.006ps 

a  < —  a  +  1 
timel 
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*10117 


ion  4-  ion7 
x  <-  0 
y  4-  0 

while  x  <  last(time) 
if  ionx  £  iontol 

if  y<iast(tspark) 

if  time  > 
x 

timex-l'iontol  +  timVonx_i  -  timex,iontol  “  i°nx  timVl 

*ion  * - : - : - - 

y  ion  ion 

X-l  X 

x  4-  x  +  skip 
y  <-  y  +  1 

x  4—  x  4-  1  otherwise 
x  4-  x  4-  1  otherwise 
x  4—  x  +  1  otherwise 
a  4-  0 
b  4-0 

while  a  <  last(tspark) 

,f  *ionb  <  tsparka+1 

timela^tionb 

a  4-  a  +  1 
b  4-  b  +  1 

otherwise 

timela  4-  0.007ps 

a  4—  a  +  I 
timel 
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tion8  :“ 


ion  <—  ion8 
x  <-  0 
y  <r-  0 

while  x  <  last(time) 
I  if  ionx  <  ionto] 


if  y  <  last(tspark) 
if  time  >  t, 


spark. 


time  , -iont„i  +  time__-ion  ,  -  time  -ionto|  -  ion  -time 


lion  . 


x-l'wntol 


x-1 


x-1 


ion  ,  -  ion 
x-I  x 


x  «-  x  +  skip 
y<-y+  1 

x  4—  x  +  1  otherwise 
x  4-  x  +  1  otherwise 
x  4—  x  +  1  otherwise 
a  4-  0 
b  <-  0 

while  a  <  last(tspark) 

lf  tionb  <  tsParka+  l 
timela  <“  tionb 

a  4-  a  +  1 
b  4-  b  +  1 

otherwise 

timel  4-  0.006ps 

£L 

a  4-  a  +  1 
timel 
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Determining  Velocities  from  time  and  distance 


dist, 


0 


dist 


dist. 


dist. 


Vj2 V23  ■— -  V34  := -  V45  :=  — — — — 

*ion2  “  ^ionl  *ion3  -  lion4  ”  ^ion3  ^ionS  ~  tion4 


dist. 


disU 


dist^ 


v56  := 


v67 := 


tion6  *ion5  ^ion7  ^onfi 

Assign  Data  to  Matrix  for  Exporting 


v78  := 


{ion8  lion7 


(o>  tspark  <l)  tionl  <2>  tion2  (3)  lion3 


ans  := 


ans 


ans  := 


ans  := 


{4)  tion4 


ans  := 


ans 


(5)  ^onS  (6>  lion6 


ans  := 


<7)  tion?  (s)  tion8 

ans  := -  ans  := - 

s  s 


ans 


(9)  v12  <10>  v23 


m 

s 


ans 


m 

s 


(ll)  v34  <12>  v45 

ans  := -  ans  := - 

m  m 


ans 


(l3>  v56 


m 

s 


ans 


<14> 


v67 

m 

s 


ans 


<15)  v78 


m 

s 
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Appendix  E.  Flash  Vaporization  System  Heat  Transfer  Calculations 
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Constants  and  Unit  Conversions 


bar  :=  14.503773773021psi  kJ  :=  1000J 
R:=  8.314— cr  :=  5.67T0”8  W 


MJ  :=  106J  MPa  :=  106Pa 


kmol  := 


mol-K 


2„4 
m  K 


Importing  surrogate  data 

^surjp8  :=  ■  ^surjp8  ,= 

...\mujp8.xls  ...\kjp8.xls 


Psurjp8  • 


...\rojp8.xls 


cPsurjp8 


Acpjp8.xls 


6  Dec  04  FN  3.2  H2-AIR  EOT 
Water  volumetric  flow 


hottest :  7 


4L 


15min 
Importing  test  data 
data  := 

4-12-6.xls 


Tamb:=  290K 


Tin  :=  data^-K 
Tout:=  data<6>'K 


Max  outlet  temperature  and  inlet  temperature  at  same  time 


^fluidout :  max(Tout) 


Tfluidout  =  512-073K 


Tfluidin :  Tin|match(Tfluidout,Tout)| 
^fluidin  =  301.986  K 


Time  :=  data^-s 


mol 
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600 

- 1 - f - 1 - 1 - 

Tin  500 

1 

V 

V 

\ 

\ 

X 

X 

X 

l 

1 

1 

1 

1 

1 

1 

t 

1 

1 

1 

1 

i 

t 

1 

I 

1 

1 

« 

1 

f 

_ 1 _ 

K 

✓ 

/ 

/ 

^  400 

—  /  — 

Tout 

1 

1 

K  300 

/ 

200 

_ 1 _ 1 _ 1 _ 1 

200  1 - 1 - 1 - 1 - 1 - 1 

0  100  200  300  400  500 

Time 


s 

Heat  Exchanger  Dimensions 

di5  :=  2.067in  dio  :=  2.375in  I  :=  30in  doi  :=  2.635in  dQ0  :=  2.875in 

Operating  Parameters 

frequency  :=  10Hz  Vol^g  :=  245in3  FF  :=  1  Tubes  :=  2 
Patm:=latm  W=  290K  *  -  1-0 

0  Properties  Information - - - 

0  Mass  Flow  Calculator———— - 

0  Energy  Required - 
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10  Dec  04  FN  2.4  H2-AIR  BOT 


hottest  • 


4L 

12min 


data  := 

l-12-10.xls 


^amb  ■“  290K 

<9> 

Tin  :=  data  -K 

Tout  -  K 

Tfluidout :=  max(Tout) 

Tfluidin-  Tin|match(Tfluidout,Tout; 

Tfluidout  ~  53 1.207  K 

Tfluidin  =  291-961 K 

(6> 

Time  :=  data  -s 

djj  -  2.067in  dlo  -  2.375in  1  s.  30in  dol  :=  2.635in  d00  :=  2.875in 

3 

frequency  :=  10Hz  Vol^g  :=  245 in  FF  :=  1 

Tubes  :=  2  Patm  :=  latm  T^  :=  290K  *  :=  1.0 

0  Properties - — — - — - ■ - 

0  Mass  Flow  Calculator  — - — - 

0  Energy  Required — — - — - 
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P  test '  ^nterp(^water » P  water » ^fluidin) 
hottest :=  P  test' hottest 


HeatTransfer  := 


mdottest4*nterFj  TWater,cpwater 


(' ^in  +  ^out)  1 

to 

'((Tout  ~  ^in)) 

hottest  ”  0333  mfn 

mfdot  =  0-401” 
min 


EnergyRequired  =  3.693  kW 


HeatTransfer 

kW 
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5-1-12  Avgas-H20  phi  1 .04-1 .1  15hz  .837kgmin 
data  := 

5-1-12.xls 


Tamb:=  290K 

Tjn  :=  data^-K 

<5> 

Tout:= data  ‘K 

Tfluidout :=  max(Tout) 

Tfluidin :  Tin|match(Tfluidout)Tout)0| 

Tfluidout  =  546-773K 

^fiuidin  ~  296.349  K 

Time  :=  data  s 


mdottest -837'  ■ 


du  :=  2.067in  dio  :=  2.375in  1  :=  30in 


Volt,,^  :=  245in3  FF:=  1 


frequency  :=  15Hz  '  ”  tu^e 

[f -459.67^  (  0 


Theatedair:=  linH 


1340.33 )  ^1000 


K,  120 


dQj  :=  2.635in  d0Q  :=  2.875in 

Tubes  :=  2  Patm  1  atm  <J)  :=  mean 

^heatedair'^  +  ^fluidin'^ 

Tmix :=  72 


[3  Properties - 

0  Mass  Flow  Calculator 
[3  Energy  Required - 
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5-1-26  Avgas-H20  phi  1.05-1.1  15hz  1.1541bmin 


data  := 

5-1-26  1.xls 


amb :=  290K 

,  .  <8>  v 
jn  :=  data  K 

,  (9>„ 

out :=  data  -K 

fluidout max(Tout) 

■^fluidin :  Tin|match(Tfluidout,Tout 

fluidout  =  479.581  K 

Tfluidin  =  289086  K 

<0> 

Time  :=  data  -s 


mdottest 


:=  1.154- 


lb 


min 


djj  :=  2.067in  dio  :=  2.375in  1  :=  15in 


d0j  ;=  2.635in  dQ0  :=  2.875in 


frequency  :=  15Hz  Vol^g  :=  245in' 


3‘ 


FF  :=  1  Tubes  ;=  2  Patm  :=  latm  <}»  :=  1.01 


Tair  :=  linH  I 


~( -459.67 (  0 


LV  1340.33^  UOOO 


K,  120 


Tair-15  +  T, 


fluidin' 


atm 
1 


mix  • 
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Tmix  =  319-979  K 


0  Properties - 

□'Mass  Flow  Calculator 
□  Energy  Required - 
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HeatTransfer  := 


mdottest'^nlerP]  ^water,cp water’ 


(Tin+Tout)~ 

2 

■((^out  ”  ^in)) 

mdottest  0.523 


”fiot=0'572'5; 

min 


EnergyRequired  =  3.998  kW 


HeatTransfer 

kW 


HeatTransfer 

EnergyRequired 
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5-1-26  Avgas-H20  phi  1 .05-1.1  15hz  .802  Ibmin 
data  := 

5-1-26  2.xls 


Tamb:=  290K 

(o) 

Tin  :=  data  -K 

(g) 

Tout:=  dataW-K 

^fluidout max(T0Ut) 

Tfluidin  := 

Tfluidout  ~  538.763  K 

I'fluidin  = 

Time data^-s 

„„„  lb 

hottest  :=  -802-— 
min 


600 


da  :=  2.067in  dio  :=  2.375in  1  :=  15in 


..  3 


d0j  :=  2.635in  dQ0  :=  2.875in 


frequency  :=  15Hz  Vo,tube  245  in  FF:- 1  Tubes  2  Patm:-latm  <j>  :=  1.05 
^-459.67  (O'] 


Taa  :=  lintei 


1340.33 )  UOOOy 


K,  120 


^air"  ^  +  Guidin'  1 

Tmix  := - - -  Tmix  =  320.084  K 


0  Properties - - 

0  Mass  Flow  Calculator 
0  Energy  Required - 
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5-2-1  JP8-JP8 


data  := 

5-2-1  .xls 


Tamb  :=  290K 

~  .  ,  <28>  „ 
Tjn  :=  data  K 

-  .  ,  <29>  „ 

Tout  :=  data  *K 


Tfluidout  ’  max(2'out) 

Tfluidin  •  Tin|match(Tfluidout!Tout)0| 

Tfluidout  =  632.054  K 

Tfluidin  =  291.548  K 

a  .  ^33>  lb 

“flowmeter  data  — 
min 

Tfuelinlet :=  ^ 

AP  :=  (data*'35'*  -  data'34^)-psi 

(o> 

Time  :=  data  -s 


djj  :=  2.067in  djQ  :=  2.375in  1  :=  15in 


(  9.65  N 

'l0> 

frequency  :=  linterp 

13.17 

15 

,  10.65 

v  16.7  , 

,20; 

dQj  :=  2.635in  dQ0  :=  2.875in 


.  3 


Voltube  :=  245in 


FF  :=  1 


Patm latm  <t):= 


mix  ■ 


^air +  ^fluidin'^ 
16 


mix 


Tubes  :=  2 
1.05 


=  387.842  K 


|T]  Properties - 

00  Mass  Flow  Calculator 
[3  Energy  Required - 
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Time 


172 


Psurr  ■  linterP(Tsur’  Psur’  ^fuelinlet) 
hi  :=  20 


lo:=  10 

avg  :=  hi  -  lo  +  1 


CRC  density  data  used  for  FN  calibration 

( 293  ^ 


^805  —  ^ 
3 
m 


PCRC  := 


760 


kg 


lCRC 


1,358 J 


K 


linterp 


hi  ^ 

I  T«“i 


TCRC’PCRC> 


i  =  lo 


avg  ) 


=  710.775-^- 


m 


hi 


Y.  H^owmeter. 


i  =  lo 


FN  := 


avg 


hi 

1 - I 

z  ^ 

i  =lo 

linterp 

V 

>1  avg  A 

hi 

Z  T°uti 


TCRC>PCRC’ 


i  -  lo 


avg  ) 


Given 


2 

mdotxx 


dPxx’^xx 


FN  =  0.0000001478  m2 


Peak  : 


- > 

2 

mflowmeter 

AP-FN2 


173 


5-2-3  JP8-JP8 
data  := 

5-2-3.xls 

Tamb:=  290K 

Tin  :=  data<n)-K 
Tout:=  data(l4>.K 
Tfluidout :=  max(Tout) 


Tfluidout  =  7G5-424K 


mftowmeter :=  data 


AP  :=  (data*'23''  -data^J-psi 


<17>  lb 

min 

S2V\ 


Tfluidin  ’  Tin|match(Tfluidout,Tout 
Tfluidin  =  285.972  K 

(13) 

^fuelinlet 


Time  :=  data^-s 


djj  :=  2.067in  djQ  :=  2.375in  1  :=  15in 


doi  :=  2.635in  dQ0  :=  2.875in 


frequency  :=  13.88Hz  ’  U,tube 

/ -459.67^  f  0  > 

1 1340.33  1000, 


Tajr  :=  limed 


Volfllhe  :=  245in  FF  :=  1  Tubes  :=  2  P, 
K,  250 


Tajr'15  +  Tfi„;,j;„-1 


atm-1  atin  $  :=  1.00 

fluidin- 


mix  • 


16 


Tmix  =  387.493  K 


0  Properties - 

0  Mass  Flow  Calculator 

{►]  Energy  Required - 
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Psurr '  ^nterp(Tsur’  Psur’  ^fuelinlet) 
hi  :=  20 
lo  :=  10 

avg  :=  hi  -  lo  +  1 


CRC  density  data  used  for  FN  calibration 


f  v„  \ 

805 


PCRC  :~ 


kg 


m 


760^- 

V  m3; 


lCRC  : 


(  293 
358 


•K 


linterp 


hi  ^ 

Z  T°uti 


tcrc>Pcrc> 


i  =  lo 


=  684.799 


kg 


avg  ) 


m 


hi 


/  ,  “flowmeter. 


i  =  lo 


FN:= 


avg 


hi 

- 7 

i  ^ 

i  =  lo 

linterp 

\ 

^  avg  > 

TCRC>PCRC 


hi 


i  =  lo 


lout. 


avg  ) 


Given 


FN  = 
xx 


mdotxx 


Find(pxx) 


m, 


dotxx 


dPxx^xx 


Pcalc  • 


“flowmeter 


AP- 


’■FN2 


FN  =  0.0000001512m2 
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5-2-4  JP8-JP8 


data  := 


5-2-4.xls 


Tamb:.290K 


Tin:=  data(l5)-K 
Tout :=  data<16>.K 


^fluidout max(Tout) 


Tfluidout  “  769.571  K 


mflowmeter  •  ^ata 


<12)  kg 


min 


AP  :=  (data^  -  data^j-psi 


<0> 

Time data  -s 


Tfluidin :  Tin|match(Tf]ujdout,Tout)0| 
Tfluidin  =  293244  K 
^fuelinlet ^ata 


dj'  :=  2.067in  dj0  :=  2.375in  1  :=  30in 


frequency  :=  12Hz  ^°*tube  •  245in 
Patm  :=  latm  <t>  :=  100 
Tair :=:  linterP| 


doi  :=  2.635in  dQ0  :=  2.875in 


FF:=  1  Tubes  :=  2 


7-459.67 "I 

f  0  ^  ] 

J 

K.250 

LU340.33;  ' 

- 1 

0s 

0 

0 

r-H 

Tair' 15  +  Tfluidin' 1 


‘mix  • 


16 


T  ■  = 
Xmix 


387.948  K 


0  Properties - 

[T]  Mass  Flow  Calculator 
[►]  Energy  Required - 
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HeatTransfer  := 


mfdofimted  Tsur>cpsur 


(' ^in  +  ^out)  1 


((^out  ^in)) 


lb 

"^dot  =  a783~ 
min 


HeatTransfer 

EnergyRequired 
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pSUrr  •  linterp^Tsur,  psur,  Tfuejjnjetj 

hi  :=  20 


lo  :=  10 

avg  :=  hi  -  lo  +  1 


CRC  density  data  used  for  FN  calibration 


PCRC  := 


f  kg  ^ 
805  -a. 

3 

m 


760 


kg 


V  to.  ) 


lCRC 


293 

358 


■K 


hi  A 


out- 


i  =  lo 


linterp^TCRC,pCRC)  ^ 


=  699.549-^- 


m 


hi 


71  m  flowmeter. 


i  =  lo 


Given 


FN 


FNxx  = 


mdotxx 


Find(Pxx) 


m. 


2 


'dotxx 


dpxx-FNxx 


Pcalc  ■ 


'’’flowmeter 


AP' 


•FN2 


0.0000002021  m2 
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